We consider the tilting instability of a magnetically confined spheromak using 3D MHD and relativistic PIC calculations with an application to astrophysical plasmas, specifically those occurring in magnetar magnetospheres. The instability is driven by the counter alignment of the spheromak's intrinsic magnetic dipole with the external magnetic field. Initially the spheromak rotates - tilts - trying to lower its magnetic potential energy. As a result a current sheet forms between the internal magnetic field of a spheromak and the confining field. Magnetic reconnection sets in; this leads to the annihilation of the newly counter-aligned magnetic flux of the spheromak. This occurs on few Alfv\'en time scales. In the case of higher order (second order) spheromak, the internal core is first pushed out of the envelope, resulting in formation of two nearly independent tilting spheromaks. Thus, the magnetically twisted outer shell cannot stabilize the inner core. During dissipation, helicity of the initial spheromak is carried away by torsional Alfv\'en waves, violating the assumptions of the Taylor relaxation theorem. In applications to magnetars' giant flares, fast development of tilting instabilities, and no stabilization of the higher order spheromaks, make it unlikely that trapped spheromaks are responsible for the tail emission lasting hundreds of seconds.

We investigate with relativistic MHD simulations the dissipation physics of BL Lac jets, by studying the synchrotron polarization signatures of particles accelerated by the kink instability in a magnetically-dominated plasma column. The nonlinear stage of the kink instability generates current sheets, where particles can be efficiently accelerated via magnetic reconnection. We identify current sheets as regions where s = J d/B is above some predefined threshold (where B is the field strength, J the current density and d the grid scale), and assume that the particle injection efficiency scales as proportional to the square of the current. X-ray emitting particles have short cooling times, so they only probe the field geometry of their injection sites. In contrast, particles emitting in the optical band, which we follow self-consistently as they propagate away from their injection sites while cooling, sample a larger volume, and so they may be expected to produce different polarimetric signatures. We find that the degree of polarization is roughly the same between X-ray and optical bands, because even the optical-emitting particles do not travel far from the current sheet where they were injected, due to lack of sufficient kink-generated turbulence. The polarization angle shows a different temporal evolution between the two bands, due to the different regions probed by X-ray and optical emitting particles. In view of the upcoming IXPE satellite, our results can help constrain whether kink-induced reconnection (as opposed to shocks) can be the source of multi-wavelength emission from BL Lacs.

Plasmoids -- magnetized quasi-circular structures formed self-consistently in reconnecting current sheets -- were previously considered to be the graveyards of energetic particles. In this paper, we demonstrate the important role of plasmoids in shaping the particle energy spectrum in relativistic reconnection (i.e., with upstream magnetization $\sigma_{\rm up} \gg 1$). Using two dimensional particle-in-cell simulations in pair plasmas with $\sigma_{\rm up}=10$ and $100$, we study a secondary particle energization process that takes place inside compressing plasmoids. We demonstrate that plasmoids grow in time, while their interiors compress, amplifying the internal magnetic field. The magnetic field felt by particles injected in an isolated plasmoid increases linearly with time, which leads to particle energization as a result of magnetic moment conservation. For particles injected with a power-law distribution function, this energization process acts in such a way that the shape of the injected power law is conserved, while producing an additional non-thermal tail $f(E)\propto E^{-3}$ at higher energies followed by an exponential cutoff. The cutoff energy, which increases with time as $E_{\rm cut}\propto\sqrt{t}$, can greatly exceed $\sigma_{\rm up} m_e c^2$. We analytically predict the secondary acceleration timescale and the shape of the emerging particle energy spectrum, which can be of major importance in certain astrophysical systems, such as blazar jets.

Electromagnetic precursor waves generated by the synchrotron maser instability at relativistic magnetized shocks have been recently invoked to explain the coherent radio emission of Fast Radio Bursts. By means of two-dimensional particle-in-cell simulations, we explore the properties of the precursor waves in relativistic electron-positron perpendicular shocks as a function of the pre-shock magnetization $\sigma\gtrsim 1$ (i.e., the ratio of incoming Poynting flux to particle energy flux) and thermal spread $\Delta \gamma \equiv kT/mc^2=10^{-5}-10^{-1}$. We measure the fraction $f_\xi$ of total incoming energy that is converted into precursor waves, as computed in the post-shock frame. At fixed magnetization, we find that $f_{\xi}$ is nearly independent of temperature as long as $\Delta \gamma \lesssim 10^{-1.5}$ (with only a modest decrease of a factor of three from $\Delta\gamma =10^{-5}$ to $\Delta\gamma =10^{-1.5}$), but it drops by nearly two orders of magnitude for $\Delta \gamma \gtrsim 10^{-1}$. At fixed temperature, the scaling with magnetization $f_\xi\sim 10^{-3}\,\sigma^{-1}$ is consistent with our earlier one-dimensional results. For our reference $\sigma=1$, the power spectrum of precursor waves is relatively broad (fractional width $\sim 1-3$) for cold temperatures, whereas it shows pronounced line-like features with fractional width $\sim 0.2$ for $10^{-3} \lesssim \Delta \gamma \lesssim 10^{-1.5} $. For $\sigma\gtrsim 1$, the precursor waves are beamed within an angle $\simeq \sigma^{-1/2}$ from the shock normal (as measured in the post-shock frame), as required so they can outrun the shock. Our results can provide physically-grounded inputs for FRB emission models based on maser emission from relativistic shocks.

Synchrotron emission from astrophysical nonthermal sources usually assumes that the emitting particles are isotropic. By means of large-scale two- and three-dimensional particle-in-cell simulations, we demonstrate that the dissipation of magnetically-dominated ($\sigma_0\gg1$) turbulence in pair plasmas leads to strongly anisotropic particle distributions. At Lorentz factors $\sim \sigma_0 \gamma_{th0}$ (here, $\gamma_{th0}$ is the initial Lorentz factor), the particle velocity is preferentially aligned with the local magnetic field; instead, the highest energy particles are roughly isotropic. This energy-dependent anisotropy leads to a synchrotron spectral flux $\nu F_\nu\propto \nu^s$ that is much harder than for isotropic particles. Remarkably, for $\sigma_0\gg1$ we find that the solid-angle-averaged spectral slope in the slow cooling regime is $s\sim 0.5-0.7$ for a wide range of turbulence fluctuations, $0.25\lesssim \delta B_{\rm rms0}^2/B_0^2\lesssim 4$, despite significant variations in the power-law energy spectrum of nonthermal particles. This is because weaker turbulence levels imprint a stronger degree of anisotropy, thereby counteracting the effect of the steeper particle spectrum. The synchrotron spectral slope may be even harder, $s\gtrsim 0.7$, if the observer is in the plane perpendicular to the mean magnetic field. Our results are independent of domain size and dimensionality. Our findings may help explaining the origin of hard synchrotron spectra of astrophysical nonthermal sources, most notably the radio spectrum of Pulsar Wind Nebulae.

We model the sample of fast radio bursts (FRB), including the newly discovered CHIME repeaters, using the synchrotron blast wave model of Metzger, Margalit & Sironi (2019). This model postulates that FRBs are precursor radiation from ultra-relativistic magnetized shocks generated as flare ejecta from a central engine collide with an effectively stationary external medium. Downward drifting of the burst frequency structure naturally arises from deceleration of the blast-wave. The data are consistent with FRBs being produced by flares of energy $E_{\rm flare} \sim 10^{43}-10^{46}(f_{\xi}/10^{-3})^{-4/5}$ erg, where $f_{\xi}$ is the maser efficiency, and minimum bulk Lorentz factors $\Gamma \approx 10^2-10^3$, which generate the observed FRBs at shock radii $r_{\rm sh} \sim 10^{12}-10^{13}$ cm. We infer upstream densities $n_{\rm ext}(r_{\rm sh}) \sim 10^{2}-10^{4}$ cm$^{-3}$ and radial profiles $n_{\rm ext} \propto r^{-k}$ showing a range of slopes $k \approx [-2,1]$ (which are seen to evolve between bursts), broadly consistent with the upstream medium being the inner edge of an ion-loaded shell released by a recent energetic flare. The burst timescales, energetics, rates, and external medium properties are consistent with repeating FRBs arising from young, hyper-active flaring magnetars, but the methodology presented is generally applicable to any central engine which injects energy impulsively into a dense magnetized medium. Uncertainties and variations of the model are discussed, including the effects of the strong electric field of the FRB wave (strength parameter $a \gg 1$) on the upstream medium. One-dimensional particle-in-cell simulations of magnetized shocks into a pair plasma are presented which demonstrate that high maser efficiency can be preserved, even in the limit $a \gg 1$ in which the FRB wave accelerates the upstream electrons to ultra-relativistic speeds.

Collisionless shocks heat electrons in the solar wind, interstellar blast waves, and hot gas permeating galaxy clusters. How much shock heating goes to electrons instead of ions, and what plasma physics controls electron heating? We simulate 2-D perpendicular shocks with a fully kinetic particle-in-cell code. For magnetosonic Mach number $\mathcal{M}_\mathrm{ms} \sim 1$-$10$ and plasma beta $\beta_\mathrm{p} \lesssim 4$, the post-shock electron/ion temperature ratio $T_\mathrm{e}/T_\mathrm{i}$ decreases from $1$ to $0.1$ with increasing $\mathcal{M}_\mathrm{ms}$. In a representative $\mathcal{M}_\mathrm{ms}=3.1$, $\beta_\mathrm{p}=0.25$ shock, electrons heat above adiabatic compression in two steps: ion-scale $E_\parallel = \vec{E} \cdot \hat{b}$ accelerates electrons into streams along $\vec{B}$, which then relax via two-stream-like instability. Shock rippling also allows quasi-static shock-normal electric fields to heat electrons; we find that quasi-static fields generally contribute half of the electron heating beyond adiabatic compression.

We present Chandra and VLA observations of GW170817 at ~521-743 days post merger, and a homogeneous analysis of the entire Chandra data set. We find that the late-time non-thermal emission follows the expected evolution from an off-axis relativistic jet, with a steep temporal decay $F_{\nu}\propto t^{-1.95\pm0.15}$ and a simple power-law spectrum $F_{\nu}\propto \nu^{-0.575\pm0.007}$. We present a new method to constrain the merger environment density based on diffuse X-ray emission from hot plasma in the host galaxy and we find $n\le 9.6 \times 10^{-3}\,\rm{cm^{-3}}$. This measurement is independent from inferences based on the jet afterglow modeling and allows us to partially solve for model degeneracies. The updated best-fitting model parameters with this density constraint are a fireball kinetic energy $E_0 = 1.5_{-1.1}^{+3.6}\times 10^{49}\,\rm{erg}$ ($E_{iso}= 2.1_{-1.5}^{+6.4}\times10^{52}\, \rm{erg}$), jet opening angle $\theta_{0}= 5.9^{+1.0}_{-0.7}\,\rm{deg}$ with characteristic Lorentz factor $\Gamma_j = 163_{-43}^{+23}$, expanding in a low-density medium with $n_0 = 2.5_{-1.9}^{+4.1} \times 10^{-3}\, \rm{cm^{-3}}$ and viewed $\theta_{obs} = 30.4^{+4.0}_{-3.4}\, \rm{deg}$ off-axis. The synchrotron emission originates from a power-law distribution of electrons with $p=2.15^{+0.01}_{-0.02}$. The shock microphysics parameters are constrained to $\epsilon_{\rm{e}} = 0.18_{-0.13}^{+0.30}$ and $\epsilon_{\rm{B}}=2.3_{-2.2}^{+16.0} \times 10^{-3}$. We investigate the presence of X-ray flares and find no statistically significant evidence of $\ge2.5\sigma$ of temporal variability at any time. Finally, we use our observations to constrain the properties of synchrotron emission from the deceleration of the fastest kilonova ejecta with energy $E_k^{KN}\propto (\Gamma\beta)^{-\alpha}$ into the environment, finding that shallow stratification indexes $\alpha\le6$ are disfavored.

Blazar emission models based on magnetic reconnection succeed in reproducing many observed spectral and temporal features, including the short-duration luminous flaring events. Plasmoids, a self-consistent by-product of the tearing instability in the reconnection layer, can be the main source of blazar emission. Kinetic simulations of relativistic reconnection have demonstrated that plasmoids are characterized by rough energy equipartition between their radiating particles and magnetic fields. This is the main reason behind the apparent shortcoming of plasmoid-dominated emission models to explain the observed Compton ratios of BL Lac objects. Here, we demonstrate that the radiative interactions among plasmoids, which have been neglected so far, can assist in alleviating this contradiction. We show that photons emitted by large, slow-moving plasmoids can be a potentially important source of soft photons to be then up-scattered, via inverse Compton, by small fast-moving, neighboring plasmoids. This inter-plasmoid Compton scattering process can naturally occur throughout the reconnection layer, imprinting itself as an increase in the observed Compton ratios from those short and luminous plasmoid-powered flares within BL Lac sources, while maintaining energy equipartition between radiating particles and magnetic fields.

Magnetized turbulence and magnetic reconnection are often invoked to explain the nonthermal emission observed from a wide variety of astrophysical sources. By means of fully-kinetic 2D and 3D PIC simulations, we investigate the interplay between turbulence and reconnection in generating nonthermal particles in magnetically-dominated pair plasmas. A generic by-product of the turbulence evolution is the generation of a nonthermal particle spectrum with a power-law energy range. The power-law slope $p$ is harder for larger magnetizations and stronger turbulence fluctuations, and it can be as hard as $p < 2$. The Larmor radius of particles at the high-energy cutoff is comparable to the size $l$ of the largest eddies. Plasmoid-mediated reconnection, which self-consistently occurs in the turbulent plasma, controls the physics of particle injection. Then, particles are further accelerated by stochastic scattering off turbulent fluctuations. The work done by parallel electric fields - naturally expected in reconnection layers - is responsible for most of the initial energy increase, and is proportional to the magnetization $\sigma$ of the system, while the subsequent energy gain, which dominates the overall energization of high-energy particles, is powered by the perpendicular electric fields of turbulent fluctuations. The two-stage acceleration process leaves an imprint in the particle pitch-angle distribution: low-energy particles are aligned with the field, while the highest energy particles move preferentially orthogonal to it. The energy diffusion coefficient of stochastic acceleration scales as $D_\gamma\sim 0.1\sigma(c/l)\gamma^2$, where $\gamma$ is the particle Lorentz factor. This results in fast acceleration timescales $t_{acc}\sim (3/\sigma)\,l/c$. Our findings have important implications for understanding the generation of nonthermal particles in high-energy astrophysical sources.

We perform two-dimensional particle-in-cell simulations of reconnection in magnetically dominated electron-positron plasmas subject to strong Compton cooling. Magnetic reconnection under such conditions can operate in accretion disk coronae around black holes, which produce hard X-rays through Comptonization. Our simulations show that most of the plasma in the reconnection layer is kept cold by Compton losses and locked in magnetically dominated plasmoids with a small thermal pressure. Compton drag clears cavities inside plasmoids and also affects their bulk motions. These effects, however, weakly change the reconnection rate and the plasmoid size distribution from those in non-radiative reconnection. This demonstrates that the reconnection dynamics is governed by similar magnetic stresses in both cases and weakly affected by thermal pressure. We examine the energy distribution of particles energized by radiative reconnection and observe two distinct components. (1) A mildly relativistic peak, which results from bulk motions of cooled plasmoids. This component receives most of the reconnection power and dominates the output X-ray emission. The peak has a quasi-Maxwellian shape with an effective temperature of ~100 keV. Thus, it mimics thermal Comptonization used previously to fit hard-state spectra of accreting black holes. (2) A high-energy tail, which receives ~20% of the reconnection power. It is populated by particles accelerated impulsively at X-points or "picked up" by fast outflows from X-points. The high-energy particles immediately cool, and their inverse Compton emission explains the MeV spectral tail detected in the hard state of Cyg X-1. Our first-principle simulations support magnetic reconnection as a mechanism powering hard X-ray emission from accreting black holes.

We present a revised and complete optical afterglow light curve of the binary neutron star merger GW170817, enabled by deep Hubble Space Telescope (HST) F606W observations at $\approx\!584$ days post-merger, which provide a robust optical template. The light curve spans $\approx 110-362$ days, and is fully consistent with emission from a relativistic structured jet viewed off-axis, as previously indicated by radio and X-ray data. Combined with contemporaneous radio and X-ray observations, we find no spectral evolution, with a weighted average spectral index of $\langle \beta \rangle = -0.583 \pm 0.013$, demonstrating that no synchrotron break frequencies evolve between the radio and X-ray bands over these timescales. We find that an extrapolation of the post-peak temporal slope of GW170817 to the luminosities of cosmological short GRBs matches their observed jet break times, suggesting that their explosion properties are similar, and that the primary difference in GW170817 is viewing angle. Additionally, we place a deep limit on the luminosity and mass of an underlying globular cluster of $L \lesssim 6.7 \times 10^{3}\,L_{\odot}$, or $M \lesssim 1.3 \times 10^{4}\,M_{\odot}$, at least 4 standard deviations below the peak of the globular cluster mass function of the host galaxy, NGC4993. This limit provides a direct and strong constraint that GW170817 did not form and merge in a globular cluster. As highlighted here, HST (and soon JWST) enables critical observations of the optical emission from neutron star merger jets and outflows.

Electron acceleration during magnetic reconnection is thought to play a key role in time-variable high-energy emission from astrophysical systems. By means of particle-in-cell simulations of trans-relativistic reconnection, we investigate electron injection and acceleration mechanisms in low-$\beta$ electron-proton plasmas. We set up a diversity of density and field structures (e.g., X-points and plasmoids) by varying the guide field strength and choosing whether to trigger reconnection or let it spontaneously evolve. We show that the number of X-points and plasmoids controls the efficiency of electron acceleration, with more X-points leading to a higher efficiency. Using on-the-fly acceleration diagnostics, we also show that the non-ideal electric fields associated with X-points play a critical role in the first stages of electron acceleration. As a further diagnostic, we include two populations of test particles that selectively experience only certain components of electric fields. We find that the out-of-plane component of the parallel electric field determines the hardness of the high-energy tail of the electron energy distribution. These results further our understanding of electron acceleration in this regime of magnetic reconnection and have implications for realistic models of black hole accretion flows.

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 outline a model of the Crab Pulsar Wind Nebula with two different populations of synchrotron emitting particles, arising from two different acceleration mechanisms: (i) Component-I due to Fermi-I acceleration at the equatorial portion of the termination shock, with particle spectral index $p_I \approx 2.2$ above the injection break corresponding to $\gamma_{wind} \sigma_{wind} \sim 10^5$, peaking in the UV ($\gamma_{wind} \sim 10^2$ is the bulk Lorentz factor of the wind, $\sigma _{wind} \sim 10^3$ is wind magnetization); (ii) Component-II due to acceleration at reconnection layers in the bulk of the turbulent Nebula, with particle index $p_{II} \approx 1.6$. The model requires relatively slow but highly magnetized wind. For both components the overall cooling break is in the infra-red at $\sim 0.01$ eV, so that the Component-I is in the fast cooling regime (cooling frequency below the peak frequency). In the optical band Component-I produces emission with the cooling spectral index of $\alpha_o \approx 0.5$, softening towards the edges due to radiative losses. Above the cooling break, in the optical, UV and X-rays, Component-I mostly overwhelms Component-II. We hypothesize that acceleration at large-scale current sheets in the turbulent nebula (Component-II) extends to the synchrotron burn-off limit of $\epsilon_s \approx 100$ MeV. Thus in our model acceleration in turbulent reconnection (Component-II) can produce both hard radio spectra and occasional gamma-ray flares. This model may be applicable to a broader class of high energy astrophysical objects, like AGNe and GRB jets, where often radio electrons form a different population from the high energy electrons.

Magnetic reconnection is often invoked to explain the non-thermal radiation of relativistic outflows, including jets of active galactic nuclei (AGN). Motivated by the largely unknown plasma composition of AGN jets, we study reconnection in the unexplored regime of electron-positron-proton (pair-proton) plasmas with large-scale two-dimensional particle-in-cell simulations. We cover a wide range of pair multiplicities (lepton-to-proton number ratio $\kappa=1-199$) for different values of the all-species plasma magnetization ($\sigma=1,3$ and 10) and electron temperature ($\Theta_e\equiv kT_e/m_ec^2=0.1-100$). We focus on the dependence of the post-reconnection energy partition and lepton energy spectra on the hot pair plasma magnetization $\sigma_{e,h}$ (i.e., the ratio of magnetic to pair enthalpy densities). We find that the post-reconnection energy is shared roughly equally between magnetic fields, pairs, and protons for $\sigma_{e,h}\gtrsim 3$. We empirically find that the mean lepton Lorentz factor in the post-reconnection region depends on $\sigma, \Theta_e$, and $\sigma_{e,h}$ as $\langle \gamma_e-1\rangle \approx \sqrt{\sigma}(1+4\Theta_e)\left(1+\sigma_{e,h}/30\right)$, for $\sigma\ge1$. The high-energy part of the post-reconnection lepton energy distributions can be described by a power law, whose slope is mainly controlled by $\sigma_{e,h}$ for $\kappa \gtrsim 3-6$, with harder power laws obtained for higher magnetizations. We finally show that reconnection in pair-proton plasmas with multiplicities $\kappa \sim 1-20$, magnetizations $\sigma \sim 1-10$, and temperatures $\Theta_e \sim 1-10$ results in particle power law slopes and average electron Lorentz factors that are consistent with those inferred in leptonic models of AGN jet emission.

Using 1D and 2D particle-in-cell (PIC) simulations of a plasma with a growing magnetic field $\vec{B}$, we show that ions can be stochastically accelerated by the ion-cyclotron (IC) instability. As $\vec{B}$ grows, an ion pressure anisotropy $p_{\perp,i} > p_{||,i}$ arises, due to the adiabatic invariance of the ion magnetic moment ($p_{||,i}$ and $p_{\perp,i}$ are the ion pressures parallel and perpendicular to $\vec{B}$). When initially $\beta_i = 0.5$ ($\beta_i \equiv 8\pi p_i/|\vec{B}|^2$, where $p_i$ is the ion isotropic pressure), the pressure anisotropy is limited mainly by inelastic pitch-angle scattering provided by the IC instability, which in turn produces a non-thermal tail in the ion energy spectrum. After $\vec{B}$ is amplified by a factor $\sim 2.7$, this tail can be approximated as a power-law of index $\sim 3.4$ plus two non-thermal bumps, and accounts for $2-3\%$ of the ions and $\sim 18\%$ of their kinetic energy. On the contrary, when initially $\beta_i =2$, the ion scattering is dominated by the mirror instability and the acceleration is suppressed. This implies that efficient ion acceleration requires that initially $\beta_i \lesssim 1$. Although we focus on cases where $\vec{B}$ is amplified by plasma shear, we check that the acceleration occurs similarly if $\vec{B}$ grows due to plasma compression. Our results are valid in a sub-relativistic regime where the ion thermal energy is $\sim 10\%$ of the ion rest mass energy. This acceleration process can thus be relevant in the inner region of low-luminosity accretion flows around black holes.

Direct observations of compact objects, in the form of radiation spectra, gravitational waves from VIRGO/LIGO, and forthcoming direct imaging, are currently one of the primary source of information on the physics of plasmas in extreme astrophysical environments. The modeling of such physical phenomena requires numerical methods that allow for the simulation of microscopic plasma dynamics in presence of both strong gravity and electromagnetic fields. In Bacchini et al. (2018) we presented a detailed study on numerical techniques for the integration of free geodesic motion. Here we extend the study by introducing electromagnetic forces in the simulation of charged particles in curved spacetimes. We extend the Hamiltonian energy-conserving method presented in Bacchini et al. (2018) to include the Lorentz force and we test its performance compared to that of standard explicit Runge-Kutta and implicit midpoint rule schemes against analytic solutions. Then, we show the application of the numerical schemes to the integration of test particle trajectories in general relativistic magnetohydrodynamic (GRMHD) simulations, by modifying the algorithms to handle grid-based electromagnetic fields. We test this approach by simulating ensembles of charged particles in a static GRMHD configuration obtained with the Black Hole Accretion Code (BHAC).

Fast radio bursts (FRB) can arise from synchrotron maser emission at ultra-relativistic magnetized shocks, such as produced by flare ejecta from young magnetars. We combine PIC simulation results for the synchrotron maser with the dynamics of self-similar shock deceleration, as commonly applied to GRBs, to explore the implications for FRB emission. We assume the upstream environment into which the ultra-relativistic ejecta collides is a mildly relativistic baryon-loaded shell released following a previous flare, motivated by the high electron-ion injection rate Mdot ~ 1e19-1e21 g/s needed on larger scales to power the persistent radio nebula coincident with the repeating burster FRB 121102 and its high inferred rotation measure. The observed radio fluence peaks once the optical depth ahead of the shock to induced Compton scattering decreases to <~ few, a condition which places a GHz observer on the high frequency tail of the maser SED. Given intervals between ion shell ejection events ~1e5 s similar to the occurrence rate of the most powerful bursts from FRB 121102, we demonstrate the production of FRBs of frequency ~ 0.1-10 GHz, isotropic radiated energies ~1e37-1e40 erg and durations ~0.1-10 ms for flares of energy ~1e43-1e45 erg. Deceleration of the blast wave, and increasing transparency of the upstream medium, generates a temporal decay of the peak frequency, similar to the downward drift seen in the sub-bursts of FRB 121102 and FRB 180814.J0422+73. The delay >~ 1e5 s between major ion-injection events needed to clear sufficiently low densities around the engine for FRB emission could explain prolonged "dark" periods and clustered burst arrival times, and lead to stochastic variation in the dispersion measure. Thermal electrons heated at the shock generate a short-lived <~ 1 ms (1 s) synchrotron transient at gamma-ray (X-ray) energies, analogous to a GRB afterglow.

Active galactic nuclei (AGN) with relativistic jets, powered by gas accretion onto their central supermassive black hole (SMBH), are unique laboratories for studying the physics of matter and elementary particles in extreme conditions that cannot be realized on Earth. For a long time since the discovery of AGN, photons were the only way to probe the underlying physical processes. The recent discovery of a very high energy neutrino, IceCube-170922A, coincident with a flaring blazar, TXS 0506+056, provides the first evidence that AGN jets are multi-messenger sources; they are capable of accelerating hadrons to very high energies, while producing non-thermal EM radiation and high-energy neutrinos. This new era of multi-messenger astronomy, which will mature in the next decade, offers us the unprecedented opportunity to combine more than one messenger to solve some long-standing puzzles of AGN jet physics: How do jets dissipate their energy to accelerate particles? What is the jet total kinetic power? Where and how do jets produce high-energy emission and neutrinos? What physical mechanisms drive the particle acceleration?

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

The last decade has seen the rapid, concurrent development of new classes of energetic astrophysical transients, including Fast Radio Bursts, Superluminous Supernovae, and ultralong gamma-ray bursts. We consider the implications of an emerging consensus on the origin of these diverse phenomena as different signatures of magnetar birth. We describe how magnetar-powered radio transients from milliseconds to gigaseconds will be used to address diverse physical questions related to the stellar evolution and death, the circum- and intergalactic medium, and fundamental physics on scales ranging from the quantum to cosmological. We argue that high-sensitivity radio interferometer concepts in development will be critical to accessing this science. There is also a growing need for scientists with expertise in new statistical tools and software techniques.

This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood: we lack a predictive theory for how it responds to perturbations, how it transports momentum and energy, and how it generates and amplifies magnetic fields. Such topics are foundational to the scientific study of plasmas, and are of intrinsic interest to those who regard plasma physics as a fundamental physics discipline. But these topics are also of extrinsic interest: addressing them directly informs upon our understanding of a wide variety of space and astrophysical systems, including accretion flows around supermassive black holes, the intracluster medium (ICM) between galaxies in clusters, and regions of the near-Earth solar wind. Specific recommendations to advance this field of study are discussed.

The plasma in low-luminosity accretion flows, such as the one around the black hole at the center of M87 or Sgr A* at our Galactic Center, is expected to be collisioness and two-temperature, with protons hotter than electrons. Here, particle heating is expected to be controlled by magnetic reconnection in the transrelativistic regime $\sigma_{w}\sim 0.1$-$1$, where the magnetization $\sigma_{w}$ is the ratio of magnetic energy density to plasma enthalpy density. By means of large-scale 2D particle-in-cell simulations, we explore for a fiducial $\sigma_w=0.1$ how the dissipated magnetic energy gets partitioned between electrons and protons, as a function of $\beta_{\rm i}$ (the ratio of proton thermal pressure to magnetic pressure) and of the strength of a guide field $B_{\rm g}$ perpendicular to the reversing field $B_0$. At low $\beta_{\rm i}\;(\lesssim 0.1)$, we find that the fraction of initial magnetic energy per particle converted into electron irreversible heat is nearly independent of $B_{\rm g}/B_0$, whereas protons get heated much less with increasing $B_{\rm g}/B_0$. As a result, for large $B_{\rm g} /B_{0}$, electrons receive the overwhelming majority of irreversible particle heating (${\sim}93\%$ for $B_{\rm g} /B_{0}=6$). This is significantly different than the antiparallel case $B_{\rm g}/B_0=0$, in which electron irreversible heating accounts for only ${\sim}18\%$ of the total particle heating. At $\beta_{\rm i} \sim 2$, when both species start already relativistically hot (for our fiducial $\sigma_w=0.1$), electrons and protons each receive ${\sim}50\%$ of the irreversible particle heating, regardless of the guide field strength. Our results provide important insights into the plasma physics of electron and proton heating in hot accretion flows around supermassive black holes.

The emission process of Fast Radio Bursts (FRBs) remains unknown. We investigate whether the synchrotron maser emission from relativistic shocks in a magnetar wind can explain the observed FRB properties. We perform particle-in-cell (PIC) simulations of perpendicular shocks in cold pair plasmas, checking our results for consistency among three PIC codes. We confirm that a linearly polarized X-mode wave is self-consistently generated by the shock and propagates back upstream as a precursor wave. We find that at magnetizations $\sigma\gtrsim 1$ (i.e., ratio of Poynting flux to particle energy flux of the pre-shock flow) the shock converts a fraction $f_\xi' \approx 7 \times 10^{-4}/\sigma^2$ of the total incoming energy into the precursor wave, as measured in the shock frame. The wave spectrum is narrow-band (fractional width $\lesssim 1-3$), with apparent but not dominant line-like features as many resonances concurrently contribute. The peak frequency in the pre-shock (observer) frame is $\omega^{\prime \prime}_{\rm peak} \approx 3 \gamma_{\rm s | u} \omega_{\rm p}$, where $\gamma_{\rm s|u}$ is the shock Lorentz factor in the upstream frame and $\omega_{\rm p}$ the plasma frequency. At $\sigma\gtrsim1$, where our estimated $\omega''_{\rm peak}$ differs from previous works, the shock structure presents two solitons separated by a cavity, and the peak frequency corresponds to an eigenmode of the cavity. Our results provide physically-grounded inputs for FRB emission models within the magnetar scenario.

We compare relativistic particle integrators commonly used in plasma physics showing several test cases relevant for astrophysics. Three explicit particle pushers are considered, namely the Boris, Vay, and Higuera-Cary schemes. We also present a new relativistic fully implicit particle integrator that is energy conserving. Furthermore, a method based on the relativistic guiding center approximation is included. The algorithms are described such that they can be readily implemented in magnetohydrodynamics codes or Particle-in-Cell codes. Our comparison focuses on the strengths and key features of the particle integrators. We test the conservation of invariants of motion, and the accuracy of particle drift dynamics in highly relativistic, mildly relativistic, and non-relativistic settings. The methods are compared in idealized test cases, i.e., without considering feedback on the electrodynamic fields, collisions, pair creation, or radiation. The test cases include uniform electric and magnetic fields, $\mathbf{E}\times\mathbf{B}$-fields, force-free fields, and setups relevant for high-energy astrophysics, e.g., a magnetic mirror, a magnetic dipole, and a magnetic null. These tests have direct relevance for particle acceleration in shocks and in magnetic reconnection.

We apply the general relativistic resistive magnetohydrodynamics code {\tt BHAC} to perform a 2D study of the formation and evolution of a reconnection layer in between two merging magnetic flux tubes in Minkowski spacetime. Small-scale effects in the regime of low resistivity most relevant for dilute astrophysical plasmas are resolved with very high accuracy due to the extreme resolutions obtained with adaptive mesh refinement. Numerical convergence in the highly nonlinear plasmoid-dominated regime is confirmed for a sweep of resolutions. We employ both uniform resistivity and non-uniform resistivity based on the local, instantaneous current density. For uniform resistivity we find Sweet-Parker reconnection, from $\eta = 10^{-2}$ down to $\eta = 10^{-4}$, for a reference case of magnetisation $\sigma = 3.33$ and plasma-$\beta = 0.1$. {For uniform resistivity $\eta=5\times10^{-5}$ the tearing mode is recovered, resulting in the formation of secondary plasmoids. The plasmoid instability enhances the reconnection rate to $v_{\rm rec} \sim 0.03c$ compared to $v_{\rm rec} \sim 0.01c$ for $\eta=10^{-4}$.} For non-uniform resistivity with a base level $\eta_0 = 10^{-4}$ and an enhanced current-dependent resistivity in the current sheet, we find an increased reconnection rate of $v_{\rm rec} \sim 0.1c$. The influence of the magnetisation $\sigma$ and the plasma-$\beta$ is analysed for cases with uniform resistivity $\eta=5\times10^{-5}$ and $\eta=10^{-4}$ in a range $0.5 \leq \sigma \leq 10$ and $0.01 \leq \beta \leq 1$ in regimes that are applicable for black hole accretion disks and jets. The plasmoid instability is triggered for Lundquist numbers larger than a critical value of $S_{\rm c} \approx 8000$.

The multi-wavelength spectral and temporal variability observed in blazars set tight constraints on current theoretical emission models. Here, we investigate the relativistic magnetic reconnection process as a source of blazar emission in which quasi-spherical plasmoids, containing relativistic particles and magnetic fields, are associated with the emission sites in blazar jets. By coupling recent two-dimensional particle-in-cell simulations of relativistic reconnection with a time-dependent radiative transfer code, we compute the non-thermal emission from a chain of plasmoids formed during a reconnection event. The derived photon spectra display characteristic features observed in both BL Lac sources and flat spectrum radio quasars, with the distinction made by varying the strength of the external photon fields, the jet magnetization, and the number of pairs per proton contained within. Light curves produced from reconnection events are composed of many fast and powerful flares that appear on excess of a slower evolving envelope produced by the cumulative emission of medium-sized plasmoids. The observed variability is highly dependent upon the orientation of the reconnection layer with respect to the blazar jet axis and to the observer. Our model provides a physically motivated framework for explaining the multi-timescale blazar variability across the entire electromagnetic spectrum.

Magnetic reconnection is invoked as an efficient particle accelerator in a variety of astrophysical sources of non-thermal high-energy radiation. With large-scale two-dimensional particle-in-cell simulations of relativistic reconnection (i.e., with magnetization $\sigma\gg1$) in pair plasmas, we study the long-term evolution of the power-law slope and high-energy cutoff of the spectrum of accelerated particles. We find that the high-energy spectral cutoff does not saturate at $\gamma_{\rm cut}\sim 4\sigma$, as claimed by earlier studies, but it steadily grows with time as long as the reconnection process stays active. At late times, the cutoff scales approximately as $\gamma_{\rm cut}\propto \sqrt{t}$, regardless of the flow magnetization and initial temperature. We show that the particles dominating the high-energy spectral cutoff reside in plasmoids, and in particular in a strongly magnetised ring around the plasmoid core. The growth of their energy is driven by the increase in the local field strength, coupled with the conservation of the first adiabatic invariant. We also find that the power-law slope of the spectrum ($p=-{\rm d}\log N/{\rm d}\log \gamma$) evolves with time. For $\sigma\gtrsim10$, the spectrum is hard at early times ($p\lesssim 2$), but it tends to asymptote to $p\sim 2$; the steepening of the power-law slope allows the spectral cutoff to extend to higher and higher energies, without violating the fixed energy budget of the system. Our results demonstrate that relativistic reconnection is a viable candidate for accelerating the high-energy particles emitting in relativistic astrophysical sources.

We present recent developments on numerical algorithms for computing photon and particle trajectories in the surrounding of compact objects. Strong gravity around neutron stars or black holes causes relativistic effects on the motion of massive particles and distorts light rays due to gravitational lensing. Efficient numerical methods are required for solving the equations of motion and compute i) the black hole shadow obtained by tracing light rays from the object to a distant observer, and ii) obtain information on the dynamics of the plasma at the microscopic scale. Here, we present generalized algorithms capable of simulating ensembles of photons or massive particles in any spacetime, with the option of including external forces. The coupling of these tools with GRMHD simulations is the key point for obtaining insight on the complex dynamics of accretion disks and jets and for comparing simulations with upcoming observational results from the Event Horizon Telescope.

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 non-thermal particles in high-energy astrophysical sources.

We present the discovery of a slowly-evolving, extragalactic radio transient, FIRST J141918.9+394036, identified by comparing a catalog of radio sources in nearby galaxies against new observations from the Very Large Array Sky Survey. Analysis of other archival data shows that FIRST J141918.9+394036 faded by a factor of ~50 over 23 years, from a flux of ~26 mJy at 1.4 GHz in 1993 to an upper limit of 0.4 mJy at 3 GHz in 2017. FIRST J141918.9+394036 is likely associated with the small star-forming galaxy SDSS J141918.81+394035.8 at a redshift z=0.01957 (d=87 Mpc), which implies a peak luminosity $\nu L_\nu \gtrsim 3\times10^{38}$ erg s$^{-1}$. If interpreted as an isotropic synchrotron blast wave, the source requires an explosion of kinetic energy ~10^{51} erg some time prior to our first detection in late 1993. This explosion could plausibly be associated with a long gamma-ray burst (GRB) or the merger of two neutron stars. Alternatively, FIRST J141918.9+394036 could be the nebula of a newly-born magnetar. The radio discovery of any of these phenomena would be unprecedented. Joint consideration of the event light curve, host galaxy, lack of a counterpart gamma-ray burst, and volumetric rate suggests that FIRST J141918.9+394036 is the afterglow of an off-axis (`orphan') long GRB. The long time baseline of this event offers the best available constraint in afterglow evolution as the bulk of shock-accelerated electrons become non-relativistic. The proximity, age, and precise localization of FIRST J141918.9+394036 make it a key object for understanding the aftermath of rare classes of stellar explosion.

Radio images of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), are dominated by interstellar scattering. Previous studies of Sgr A* have adopted an anisotropic Gaussian model for both the intrinsic source and the scattering, and they have extrapolated the scattering using a purely $\lambda^2$ scaling to estimate intrinsic properties. However, physically motivated source and scattering models break all three of these assumptions. They also predict that refractive scattering effects will be significant, which have been ignored in standard model fitting procedures. We analyze radio observations of Sgr A* using a physically motivated scattering model, and we develop a prescription to incorporate refractive scattering uncertainties when model fitting. We show that an anisotropic Gaussian scattering kernel is an excellent approximation for Sgr A* at wavelengths longer than 1cm, with an angular size of $(1.380 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the major axis, $(0.703 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the minor axis, and a position angle of $81.9^\circ \pm 0.2^\circ$. We estimate that the turbulent dissipation scale is at least $600\,{\rm km}$, with tentative support for $r_{\rm in} = 800 \pm 200\,{\rm km}$, suggesting that the ion Larmor radius defines the dissipation scale. We find that the power-law index for density fluctuations in the scattering material is $\beta < 3.47$, shallower than expected for a Kolmogorov spectrum ($\beta=11/3$), and we estimate $\beta = 3.38^{+0.08}_{-0.04}$ in the case of $r_{\rm in} = 800\,{\rm km}$. We find that the intrinsic structure of Sgr A* is nearly isotropic over wavelengths from 1.3mm to 1.3cm, with a size that is roughly proportional to wavelength. We discuss implications for models of Sgr A*, for theories of interstellar turbulence, and for imaging Sgr A* with the Event Horizon Telescope.

We present new observations of the binary neutron star merger GW170817 at $\Delta t\approx 220-290$ days post-merger, at radio (Karl G. Jansky Very Large Array; VLA), X-ray (Chandra X-ray Observatory) and optical (Hubble Space Telescope; HST) wavelengths. These observations provide the first evidence for a turnover in the X-ray light curve, mirroring a decline in the radio emission at $\gtrsim5\sigma$ significance. The radio-to-X-ray spectral energy distribution exhibits no evolution into the declining phase. Our full multi-wavelength dataset is consistent with the predicted behavior of our previously published models of a successful structured jet expanding into a low-density circumbinary medium, but pure cocoon models with a choked jet cannot be ruled out. If future observations continue to track our predictions, we expect that the radio and X-ray emission will remain detectable until $\sim 1000$ days post-merger.

The dissipation of energy flux in blazar jets plays a key role in the acceleration of relativistic particles. Two possibilities are commonly considered for the dissipation processes, magnetic reconnection -- possibly triggered by instabilities in magnetically-dominated jets -- , or shocks -- for weakly magnetized flows. We consider the polarimetric features expected for the two scenarios analyzing the results of state-of-the-art simulations. For the magnetic reconnection scenario we conclude, using results from global relativistic MHD simulations, that the emission likely occurs in turbulent regions with unstructured magnetic fields, although the simulations do not allow us to draw firm conclusions. On the other hand, with local particle-in-cell simulations we show that, for shocks with a magnetic field geometry suitable for particle acceleration, the self-generated magnetic field at the shock front is predominantly orthogonal to the shock normal and becomes quasi-parallel downstream. Based on this result we develop a simplified model to calculate the frequency-dependent degree of polarization, assuming that high-energy particles are injected at the shock and cool downstream. We apply our results to HBLs, blazars with the maximum of their synchrotron output at UV-soft X-ray energies. While in the optical band the predicted degree of polarization is low, in the X-ray emission it can ideally reach 50\%, especially during active/flaring states. The comparison between measurements in the optical and in the X-ray band made during active states (feasible with the planned {\it IXPE} satellite) are expected to provide valuable constraints on the dissipation and acceleration processes.

The numerical integration of particle trajectories in curved spacetimes is fundamental for obtaining realistic models of the particle dynamics around massive compact objects such as black holes and neutron stars. Generalized algorithms capable of handling generic metrics are required for studies of both standard spacetimes (Schwarzschild and Kerr metrics) and non-standard spacetimes (e.g. Schwarzschild metric plus non-classical perturbations or multiple black hole metrics). The most commonly employed explicit numerical schemes (e.g. Runge-Kutta) are incapable of producing highly accurate results at critical points, e.g. in the regions close to the event horizon where gravity causes extreme curvature of the spacetime, at an acceptable computational cost. Here, we describe a generalized algorithm for the numerical integration of time-like (massive particles) and null (photons) geodesics in any given 3+1 split spacetime. We introduce a new, exactly energy-conserving implicit integration scheme based on the preservation of the underlying Hamiltonian, and we compare its properties with a standard fourth-order Runge-Kutta explicit scheme and an implicit midpoint scheme. We test the numerical performance of the three schemes against analytic solutions of particle and photon orbits in Schwarzschild and Kerr spacetimes. We also prove the versatility of our framework in handling more exotic metrics such as Morris-Thorne wormholes and quantum-perturbed Schwarzschild black holes. The generalized approach is also discussed in the perspective of future extensions to more complex particle dynamics, e.g. the addition of the Lorentz force acting on charged particles, which allows for test particle diagnostics in GRMHD simulations.

The accretion flow around the Galactic Center black hole Sagittarius A* (Sgr A*) is expected to have an electron temperature that is distinct from the ion temperature, due to weak Coulomb coupling in the low-density plasma. We present four two-temperature general relativistic radiative magnetohydrodynamic (GRRMHD) simulations of Sgr A* performed with the code KORAL. These simulations use different electron heating prescriptions, motivated by different models of the underlying plasma microphysics. We compare the Landau-damped turbulent cascade model used in previous work with a new prescription we introduce based on the results of particle-in-cell simulations of magnetic reconnection. With the turbulent heating model, electrons are preferentially heated in the polar outflow, whereas with the reconnection model electrons are heated by nearly the same fraction everywhere in the accretion flow. The spectra of the two models are similar around the submillimetre synchrotron peak, but the models heated by magnetic reconnection produce variability more consistent with the level observed from Sgr A*. All models produce 230~GHz images with distinct black hole shadows which are consistent with the image size measured by the Event Horizon Telescope, but only the turbulent heating produces an anisotropic `disc-jet' structure where the image is dominated by a polar outflow or jet at frequencies below the synchrotron peak. None of our models can reproduce the observed radio spectral slope, the large near-infrared and X-ray flares, or the near-infrared spectral index, all of which suggest non-thermal electrons are needed to fully explain the emission from Sgr A*.

We develop a model of gamma-ray flares of the Crab Nebula resulting from the magnetic reconnection events in highly-magnetized relativistic plasma. We first discuss physical parameters of the Crab nebula and review the theory of pulsar winds and termination shocks. We also review the principle points of particle acceleration in explosive reconnection events (Lyutikov et al. 2017a,b). It is required that particles producing flares are accelerated in highly magnetized regions of the nebula. Flares originate from the poleward regions at the base of Crab's polar outflow, where both the magnetization and the magnetic field strength are sufficiently high. The post-termination shock flow develops macroscopic (not related to the plasma properties on the skin-depth scale) kink-type instabilities. The resulting large-scales magnetic stresses drive explosive reconnection events on the light-crossing time of the reconnection region. Flares are produced at the initial stage of the current sheet development, during the X-point collapse. The model has all the ingredients needed for Crab flares: natural formation of highly magnetized regions, explosive dynamics on light travel time, development of high electric fields on macroscopic scales and acceleration of particles to energies well exceeding the average magnetic energy per particle.

Using analytical and numerical methods (fluid and particle-in-cell simulations) we study a number of model problems involving merger of magnetic flux tubes in relativistic magnetically-dominated plasma. Mergers of current-carrying flux tubes (exemplified by the two dimensional `ABC' structures) and zero total current magnetic flux tubes are considered. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration during flux mergers: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. During the X-point collapse particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $\sigma$. For plasma magnetization $\sigma \leq 10^2$ the spectrum power law index is $p> 2$; in this case the maximal energy depends linearly on the size of the reconnecting islands. For higher magnetization, $\sigma \geq 10^2$, the spectra are hard, $p< 2$, yet the maximal energy $\gamma_{max}$ can still exceed the average magnetic energy per particle, $ \sim \sigma$, by orders of magnitude (if $p$ is not too close to unity). The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced magnetic reconnection, further accelerating the particles, but proceeds at a slower reconnection rate.

The extreme properties of the gamma ray flares in the Crab Nebula present a clear challenge to our ideas on the nature of particle acceleration in relativistic astrophysical plasma. It seems highly unlikely that standard mechanisms of stochastic type are at work here and hence the attention of theorists has switched to linear acceleration in magnetic reconnection events. In this series of papers, we attempt to develop a theory of explosive magnetic reconnection in highly-magnetized relativistic plasma which can explain the extreme parameters of the Crab flares. In the first paper, we focus on the properties of the X-point collapse. Using analytical and numerical methods (fluid and particle-in-cell simulations) we extend Syrovatsky's classical model of such collapse to the relativistic regime. We find that the collapse can lead to the reconnection rate approaching the speed of light on macroscopic scales. During the collapse, the plasma particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $\sigma$. For sufficiently high magnetizations and vanishing guide field, the non-thermal particle spectrum consists of two components: a low-energy population with soft spectrum, that dominates the number census; and a high-energy population with hard spectrum, that possesses all the properties needed to explain the Crab flares.

Non-thermal electron acceleration via magnetic reconnection is thought to play an important role in powering the variable X-ray emission from radiatively inefficient accretion flows around black holes. The trans-relativistic regime of magnetic reconnection, where the magnetization $\sigma$, defined as the ratio of magnetic energy density to enthalpy density, is $\sim 1$, is frequently encountered in such flows. By means of a large suite of two-dimensional particle-in-cell simulations, we investigate electron and proton acceleration in the trans-relativistic regime. We focus on the dependence of the electron energy spectrum on $\sigma$ and the proton $\beta$ (i.e., the ratio of proton thermal pressure to magnetic pressure). We find that the electron spectrum in the reconnection region is non-thermal and can be generally modeled as a power law. At $\beta \lesssim 3 \times 10^{-3}$, the slope, $p$, is independent of $\beta$ and it hardens with increasing $\sigma$ as $p\simeq 1.8 +0.7/\sqrt{\sigma}$. Electrons are primarily accelerated by the non-ideal electric field at X-points, either in the initial current layer or in current sheets generated in between merging magnetic islands. At higher values of $\beta$, the electron power law steepens for all values of $\sigma$. At values of $\beta$ near $\beta_{\rm max}\approx1/4\sigma$, when both electrons and protons are relativistically hot prior to reconnection, the spectra of both species display an additional component at high energies, containing a few percent of particles. These particles are accelerated via a Fermi-like process by bouncing in between the reconnection outflow and a stationary magnetic island. We provide an empirical prescription for the dependence of the power-law slope and the acceleration efficiency on $\beta$ and $\sigma$, which can be used in global simulations of collisionless accretion disks.

We report deep Chandra, HST and VLA observations of the binary neutron star event GW170817 at $t<160$ d after merger. These observations show that GW170817 has been steadily brightening with time and might have now reached its peak, and constrain the emission process as non-thermal synchrotron emission where the cooling frequency $\nu_c$ is above the X-ray band and the synchrotron frequency $\nu_m$ is below the radio band. The very simple power-law spectrum extending for eight orders of magnitude in frequency enables the most precise measurement of the index $p$ of the distribution of non-thermal relativistic electrons $N(\gamma)\propto \gamma^{-p}$ accelerated by a shock launched by a NS-NS merger to date. We find $p=2.17\pm0.01$, which indicates that radiation from ejecta with $\Gamma\sim3-10$ dominates the observed emission. While constraining the nature of the emission process, these observations do \emph{not} constrain the nature of the relativistic ejecta. We employ simulations of explosive outflows launched in NS ejecta clouds to show that the spectral and temporal evolution of the non-thermal emission from GW170817 is consistent with both emission from radially stratified quasi-spherical ejecta traveling at mildly relativistic speeds, \emph{and} emission from off-axis collimated ejecta characterized by a narrow cone of ultra-relativistic material with slower wings extending to larger angles. In the latter scenario, GW170817 harbored a normal SGRB directed away from our line of sight. Observations at $t\le 200$ days are unlikely to settle the debate as in both scenarios the observed emission is effectively dominated by radiation from mildly relativistic material.

Recent X-ray observations of merger shocks in galaxy clusters have shown that the post-shock plasma is two-temperature, with the protons being hotter than the electrons. In this work, the second of a series, we investigate by means of two-dimensional particle-in-cell simulations the efficiency of electron irreversible heating in perpendicular low Mach number shocks. We consider values of plasma beta (ratio of thermal and magnetic pressures) in the range $4\lesssim \beta_{p0}\lesssim 32$ and sonic Mach number (ratio of shock speed to pre-shock sound speed) in the range $2\lesssim M_{s}\lesssim 5$, as appropriate for galaxy cluster shocks. As shown in Paper I, magnetic field amplification - induced by shock compression of the pre-shock field, or by strong proton cyclotron and mirror modes accompanying the relaxation of proton temperature anisotropy - can drive the electron temperature anisotropy beyond the threshold of the electron whistler instability. The growth of whistler waves breaks the electron adiabatic invariance, and allows for efficient entropy production. We find that the post-shock electron temperature $T_{e2}$ exceeds the adiabatic expectation $T_{e2,\rm ad}$ by an amount $(T_{e2}-T_{e2,\rm ad})/T_{e0}\simeq 0.044 \,M_s (M_s-1)$ (here, $T_{e0}$ is the pre-shock temperature), which depends only weakly on the plasma beta, over the range $4\lesssim \beta_{p0}\lesssim 32$ which we have explored, and on the proton-to-electron mass ratio (the coefficient of $\simeq 0.044$ is measured for our fiducial $m_i/m_e=49$, and we estimate that it will decrease to $\simeq 0.03$ for the realistic mass ratio). Our results have important implications for current and future observations of galaxy cluster shocks in the radio band (synchrotron emission and Sunyaev-Zel'dovich effect) and at X-ray frequencies.

Recent X-ray observations of merger shocks in galaxy clusters have shown that the post-shock plasma is two-temperature, with the protons hotter than the electrons. By means of two-dimensional particle-in-cell simulations, we study the physics of electron irreversible heating in perpendicular low Mach number shocks, for a representative case with sonic Mach number of 3 and plasma beta of 16. We find that two basic ingredients are needed for electron entropy production: (i) an electron temperature anisotropy, induced by field amplification coupled to adiabatic invariance; and (ii) a mechanism to break the electron adiabatic invariance itself. In shocks, field amplification occurs at two major sites: at the shock ramp, where density compression leads to an increase of the frozen-in field; and farther downstream, where the shock-driven proton temperature anisotropy generates strong proton cyclotron and mirror modes. The electron temperature anisotropy induced by field amplification exceeds the threshold of the electron whistler instability. The growth of whistler waves breaks the electron adiabatic invariance, and allows for efficient entropy production. We find that the electron heating efficiency displays only a weak dependence on mass ratio (less than 30 percent drop, as we increase the mass ratio from 49 up to 1600). We develop an analytical model of electron irreversible heating and show that it is in excellent agreement with our simulation results.

Hot collisionless accretion flows, such as the one in Sgr A$^{*}$ at our Galactic center, provide a unique setting for the investigation of magnetic reconnection. Here, protons are non-relativistic while electrons can be ultra-relativistic. By means of two-dimensional particle-in-cell simulations, we investigate electron and proton heating in the outflows of trans-relativistic reconnection (i.e., $\sigma_w\sim 0.1-1$, where the magnetization $\sigma_w$ is the ratio of magnetic energy density to enthalpy density). For both electrons and protons, we find that heating at high $\beta_{\rm i}$ (here, $\beta_{\rm i}$ is the ratio of proton thermal pressure to magnetic pressure) is dominated by adiabatic compression ('adiabatic heating'), while at low $\beta_{\rm i}$ it is accompanied by a genuine increase in entropy ('irreversible heating'). For our fiducial $\sigma_w=0.1$, the irreversible heating efficiency at $\beta_{\rm i}\lesssim 1$ is nearly independent of the electron-to-proton temperature ratio $T_{\rm e}/T_{\rm i}$ (which we vary from $0.1$ up to $1$), and it asymptotes to $\sim 2\%$ of the inflowing magnetic energy in the low-$\beta_{\rm i}$ limit. Protons are heated more efficiently than electrons at low and moderate $\beta_{\rm i}$ (by a factor of $\sim7$), whereas the electron and proton heating efficiencies become comparable at $\beta_{\rm i}\sim 2$ if $T_{\rm e}/T_{\rm i}=1$, when both species start already relativistically hot. We find comparable heating efficiencies between the two species also in the limit of relativistic reconnection ($\sigma_w\gtrsim 1$). Our results have important implications for the two-temperature nature of collisionless accretion flows, and may provide the sub-grid physics needed in general relativistic MHD simulations.

Plasmoids, quasi-spherical regions of plasma containing magnetic fields and high-energy particles, are a self-consistent by-product of the reconnection process in the relativistic regime. Recent two-dimensional particle-in-cell (PIC) simulations have shown that plasmoids can undergo a variety of processes (e.g. mergers, bulk acceleration, growth, and advection) within the reconnection layer. We developed a Monte Carlo (MC) code, benchmarked with the recent PIC simulations, to examine the effects of these processes on the steady-state size and momentum distributions of the plasmoid chain. The differential plasmoid size distribution is shown to be a power law, $N(w)\propto w^{-\chi}$, ranging from a few plasma skin depths to $\sim 0.1$ of the reconnection layer's length. We demonstrate numerically and analytically that the power law slope $\chi$ is linearly dependent upon the ratio of the plasmoid acceleration and growth rates and that it slightly decreases (from $\sim 2$ to $\sim 1.3$) with increasing plasma magnetization (from 3 to 50). We perform a detailed comparison of our results with those of recent PIC simulations and briefly discuss the astrophysical implications of our findings through the representative case of flaring events from blazar jets.

We discuss the role of particle-in-cell (PIC) simulations in unveiling the origin of the emitting particles in PWNe. After describing the basics of the PIC technique, we summarize its implications for the quiescent and the flaring emission of the Crab Nebula, as a prototype of PWNe. A consensus seems to be emerging that, in addition to the standard scenario of particle acceleration via the Fermi process at the termination shock of the pulsar wind, magnetic reconnection in the wind, at the termination shock and in the Nebula plays a major role in powering the multi-wavelength signatures of PWNe.

Type IIn supernovae (SNe), a rare subclass of core collapse SNe, explode in dense circumstellar media that have been modified by the SNe progenitors at their last evolutionary stages. The interaction of the freely expanding SN ejecta with the circumstellar medium gives rise to a shock wave propagating in the dense SN environment, which may accelerate protons to multi-PeV energies. Inelastic proton-proton collisions between the shock-accelerated protons and those of the circumstellar medium lead to multi-messenger signatures. Here, we evaluate the possible neutrino signal of type IIn SNe and compare with IceCube observations. We employ a Monte Carlo method for the calculation of the diffuse neutrino emission from the SN IIn class to account for the spread in their properties. The cumulative neutrino emission is found to be ~ 10 per cent of the observed IceCube neutrino flux above 60 TeV. Type IIn SNe would be the dominant component of the diffuse astrophysical flux, only if 4 per cent of all core collapse SNe were of this type and 20 to 30 per cent of the shock energy was channeled to accelerated protons. Lower values of the acceleration efficiency are accessible by the observation of a single type IIn SN as a neutrino point source with IceCube using up-going muon neutrinos. Such an identification is possible in the first year following the SN shock breakout for sources within 20 Mpc.

Non-ideal MHD effects may play a significant role in determining the dynamics, thermal properties, and observational signatures of radiatively inefficient accretion flows onto black holes. In particular, particle acceleration during magnetic reconnection events may influence black hole spectra and flaring properties. We use representative GRMHD simulations of black hole accretion flows to identify and explore the structures and properties of current sheets as potential sites of magnetic reconnection. In the case of standard and normal (SANE) disks, we find that, in the reconnection sites, the plasma beta ranges from $0.1$ to $1000$, the magnetization ranges from $10^{-4}$ to $1$, and the guide fields are weak compared to the reconnecting fields. In magnetically arrested (MAD) disks, we find typical values for plasma beta from $10^{-2}$ to $10^3$, magnetizations from $10^{-3}$ to $10$, and typically stronger guide fields, with strengths comparable to or greater than the reconnecting fields. These are critical parameters that govern the electron energy distribution resulting from magnetic reconnection and can be used in the context of plasma simulations to provide microphysics inputs to global simulations. We also find that ample magnetic energy is available in the reconnection regions to power the fluence of bright X-ray flares observed from the black hole in the center of the Milky Way.

According to magnetohydrodynamics (MHD), the encounter of two collisional magnetized plasmas at high velocity gives rise to shock waves. Investigations conducted so far have found that the same conclusion still holds in the case of collisionless plasmas. For the case of a flow-aligned field, MHD stipulates that the field and the fluid are disconnected, so that the shock produced is independent of the field. We present a violation of this MHD prediction when considering the encounter of two cold pair plasmas along a flow-aligned magnetic field. As the guiding magnetic field grows, isotropization is progressively suppressed, resulting in a strong influence of the field on the resulting structure. A micro-physics analysis allows to understand the mechanisms at work. Particle-in-cell simulations also support our conclusions and show that the results are not restricted to a strictly parallel field.

Pulsar Wind Nebulae, Blazars, Gamma Ray Bursts and Magnetars all contain regions where the electromagnetic energy density greatly exceeds the plasma energy density. These sources exhibit dramatic flaring activity where the electromagnetic energy distributed over large volumes, appears to be converted efficiently into high energy particles and gamma-rays. We call this general process magnetoluminescence. Global requirements on the underlying, extreme particle acceleration processes are described and the likely importance of relativistic beaming in enhancing the observed radiation from a flare is emphasized. Recent research on fluid descriptions of unstable electromagnetic configurations are summarized and progress on the associated kinetic simulations that are needed to account for the acceleration and radiation is discussed. Future observational, simulation and experimental opportunities are briefly summarized.

Powerful flares from blazars with short ($\sim$ min) variability timescales are challenging for current models of blazar emission. Here, we present a physically motivated ab initio model for blazar flares based on the results of recent particle-in-cell (PIC) simulations of relativistic magnetic reconnection. PIC simulations demonstrate that quasi-spherical plasmoids filled with high-energy particles and magnetic fields are a self-consistent by-product of the reconnection process. By coupling our PIC-based results (i.e., plasmoid growth, acceleration profile, particle and magnetic content) with a kinetic equation for the evolution of the electron distribution function we demonstrate that relativistic reconnection in blazar jets can produce powerful flares whose temporal and spectral properties are consistent with the observations. In particular, our model predicts correlated synchrotron and synchrotron self-Compton flares of duration of several hours--days powered by the largest and slowest moving plasmoids that form in the reconnection layer. Smaller and faster plasmoids produce flares of sub-hour duration with higher peak luminosities than those powered by the largest plasmoids. Yet, the observed fluence in both types of flares is similar. Multiple flares with a range of flux-doubling timescales (minutes to several hours) observed over a longer period of flaring activity (days or longer) may be used as a probe of the reconnection layer's orientation and the jet's magnetization. Our model shows that blazar flares are naturally expected as a result of magnetic reconnection in a magnetically-dominated jet.

Blobs, or quasi-spherical emission regions containing relativistic particles and magnetic fields, are often assumed ad hoc in emission models of relativistic astrophysical jets, yet their physical origin is still not well understood. Here, we employ a suite of large-scale two-dimensional particle-in-cell simulations in electron-positron plasmas to demonstrate that relativistic magnetic reconnection can naturally account for the formation of quasi-spherical plasmoids filled with high-energy particles and magnetic fields. Our simulations extend to unprecedentedly long temporal and spatial scales, so we can capture the asymptotic physics independently of the initial setup. We characterize the properties of the plasmoids that are continuously generated as a self-consistent by-product of the reconnection process: they are in rough energy equipartition between particles and magnetic fields; the upper energy cutoff of the plasmoid particle spectrum is proportional to the plasmoid width w, corresponding to a Larmor radius ~0.2 w; the plasmoids grow in size at ~0.1 of the speed of light, with most of the growth happening while they are still non-relativistic (first they grow); their growth is suppressed once they get accelerated to relativistic speeds by the field line tension, up to the Alfven speed (then they go). The largest plasmoids, whose typical recurrence interval is ~2.5 L/c, reach a characteristic size w ~ 0.2 L independently of the system length L, they have nearly isotropic particle distributions and they contain the highest energy particles, whose Larmor radius is ~0.03 L. The latter can be regarded as the Hillas criterion for relativistic reconnection. We briefly discuss the implications of our results for the high-energy emission from relativistic jets and pulsar winds.

Collisionless shocks, that is shocks mediated by electromagnetic processes, are customary in space physics and in astrophysics. They are to be found in a great variety of objects and environments: magnetospheric and heliospheric shocks, supernova remnants, pulsar winds and their nebul\ae, active galactic nuclei, gamma-ray bursts and clusters of galaxies shock waves. Collisionless shock microphysics enters at different stages of shock formation, shock dynamics and particle energization and/or acceleration. It turns out that the shock phenomenon is a multi-scale non-linear problem in time and space. It is complexified by the impact due to high-energy cosmic rays in astrophysical environments. This review adresses the physics of shock formation, shock dynamics and particle acceleration based on a close examination of available multi-wavelength or in-situ observations, analytical and numerical developments. A particular emphasize is made on the different instabilities triggered during the shock formation and in association with particle acceleration processes with regards to the properties of the background upstream medium. It appears that among the most important parameters the background magnetic field through the magnetization and its obliquity is the dominant one. The shock velocity that can reach relativistic speeds has also a strong impact over the development of the micro-instabilities and the fate of particle acceleration. Recent developments of laboratory shock experiments has started to bring some new insights in the physics of space plasma and astrophysical shock waves. A special section is dedicated to new laser plasma experiments probing shock physics

We develop a model of particle acceleration in explosive reconnection events in relativistic magnetically-dominated plasmas and apply it to explain gamma-ray flares from the Crab Nebula. The model relies on development of current-driven instabilities on macroscopic scales (not related to plasma skin depths). Using analytical and numerical methods (fluid and particle-in-cell simulations), we study a number of model problems in relativistic magnetically-dominated plasma: (i) we extend Syrovatsky's classical model of explosive X-point collapse to magnetically-dominated plasmas; (ii) we consider instability of two-dimensional force-free system of magnetic flux tubes; (iii) we consider merger of two zero total poloidal current magnetic flux tubes. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization. The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced reconnection, further accelerating the particles, but proceeds at a slower reconnection rate. Crab flares result from the initial explosive stages of magnetic island mergers of magnetic flux tubes produced in the bulk of nebula at intermediate polar regions. The post-termination shock plasma flow in the wind sectors with mild magnetization naturally generates large-scale highly magnetized structures. Internal kink-like instabilities lead to the formation of macroscopic current-carrying magnetic flux tubes that merge explosively.

Several supernovae (SNe) with an unusually dense circumstellar medium (CSM) have been recently observed at radio frequencies. Their radio emission is powered by relativistic electrons that can be either accelerated at the SN shock (primaries) or injected as a by-product (secondaries) of inelastic proton-proton collisions. We investigate the radio signatures from secondary electrons, by detailing a semi-analytical model to calculate the temporal evolution of the distributions of protons, primary and secondary electrons. With our formalism, we track the cooling history of all the particles that have been injected into the emission region up to a given time, and calculate the resulting radio spectra and light curves. For a SN shock propagating through the progenitor wind, we find that secondary electrons control the early radio signatures, but their contribution decays faster than that of primary electrons. This results in a flattening of the light curve at a given radio frequency that depends only upon the radial profiles of the CSM density and of the shock velocity, $v_0$. The relevant transition time at the peak frequency is $\sim 190 \, K_{\rm ep,-3}^{-1} A_{\rm w, 16} \beta_{0, -1.5}^{-2}\,{\rm d}$, where $A_{\rm w}$ is the wind mass-loading parameter, $\beta_0=v_0/c$ and $K_{\rm ep}$ is the electron-to-proton ratio of accelerated particles. We explicitly show that late peak times at 5 GHz (i.e., $t_{\rm pk} \gtrsim 300-1000$ d) suggest a shock wave propagating in a dense wind ($A_{\rm w} \gtrsim 10^{16}-10^{17}$ gr cm$^{-1}$), where secondary electrons are likely to power the observed peak emission.

We formulate a magnetohydrodynamic-particle-in-cell (MHD-PIC) method for describing the interaction between collisionless cosmic ray (CR) particles and a thermal plasma. The thermal plasma is treated as a fluid, obeying equations of ideal MHD, while CRs are treated as relativistic Lagrangian particles subject to the Lorentz force. Backreaction from CRs to the gas is included in the form of momentum and energy feedback. In addition, we include the electromagnetic feedback due to CR-induced Hall effect that becomes important when the electron-ion drift velocity of the background plasma induced by CRs approaches the Alfv\'en velocity. Our method is applicable on scales much larger than the ion inertial length, bypassing the microscopic scales that must be resolved in conventional PIC methods, while retaining the full kinetic nature of the CRs. We have implemented and tested this method in the Athena MHD code, where the overall scheme is second-order accurate and fully conservative. As a first application, we describe a numerical experiment to study particle acceleration in non-relativistic shocks. Using a simplified prescription for ion injection, we reproduce the shock structure and the CR energy spectra obtained with more self-consistent hybrid-PIC simulations, but at substantially reduced computational cost. We also show that the CR-induced Hall effect reduces the growth rate of the Bell instability and affects the gas dynamics in the vicinity of the shock front. As a step forward, we are able to capture the transition of particle acceleration from non relativistic to relativistic regimes, with momentum spectrum $f(p)\sim p^{-4}$ connecting smoothly through the transition, as expected from the theory of Fermi acceleration.

We review the physics of relativistic shocks, which are often invoked as the sources of non-thermal particles in pulsar wind nebulae (PWNe), gamma-ray bursts (GRBs), and active galactic nuclei (AGN) jets, and as possible sources of ultra-high energy cosmic-rays. We focus on particle acceleration and magnetic field generation, and describe the recent progress in the field driven by theory advances and by the rapid development of particle-in-cell (PIC) simulations. In weakly magnetized or quasi parallel-shocks (where the magnetic field is nearly aligned with the flow), particle acceleration is efficient. The accelerated particles stream ahead of the shock, where they generate strong magnetic waves which in turn scatter the particles back and forth across the shock, mediating their acceleration. In contrast, in strongly magnetized quasi-perpendicular shocks, the efficiencies of both particle acceleration and magnetic field generation are suppressed. Particle acceleration, when efficient, modifies the turbulence around the shock on a long time scale, and the accelerated particles have a characteristic energy spectral index of ~ 2.2 in the ultra-relativistic limit. We discuss how this novel understanding of particle acceleration and magnetic field generation in relativistic shocks can be applied to high-energy astrophysical phenomena, with an emphasis on PWNe and GRB afterglows.

Observations of gamma-ray-bursts and jets from active galactic nuclei reveal that the jet flow is characterized by a high radiative efficiency and that the dissipative mechanism must be a powerful accelerator of non-thermal particles. Shocks and magnetic reconnection have long been considered as possible candidates for powering the jet emission. Recent progress via fully-kinetic particle-in-cell simulations allows us to revisit this issue on firm physical grounds. We show that shock models are unlikely to account for the jet emission. In fact, when shocks are efficient at dissipating energy, they typically do not accelerate particles far beyond the thermal energy, and vice versa. In contrast, we show that magnetic reconnection can deposit more than 50% of the dissipated energy into non-thermal leptons as long as the energy density of the magnetic field in the bulk flow is larger than the rest mass energy density. The emitting region, i.e., the reconnection downstream, is characterized by a rough energy equipartition between magnetic fields and radiating particles, which naturally accounts for a commonly observed property of blazar jets.

This review discusses the physics of magnetic reconnection, a process in which the magnetic field topology changes and magnetic energy is converted to kinetic energy, in pair plasmas in the relativistic regime. We focus on recent progress in the field driven by theory advances and the maturity of particle-in-cell codes. This work shows that fragmentation instabilities at the current sheet can play a critical role in setting the reconnection speed and affect the resulting particle acceleration, anisotropy, bulk flows, and radiation. Then, we discuss how this novel understanding of relativistic reconnection can be applied to high-energy astrophysical phenomena, with an emphasis on pulsars, pulsar wind nebulae, and active galactic nucleus jets.

In the innermost regions of low-luminosity accretion flows, including Sgr A* at the center of our Galaxy, the frequency of Coulomb collisions is so low that the plasma is two-temperature, with the ions substantially hotter than the electrons. This paradigm assumes that Coulomb collisions are the only channel for transferring the ion energy to the electrons. In this work, the second of a series, we assess the efficiency of electron heating by ion velocity-space instabilities in collisionless accretion flows. The instabilities are seeded by the pressure anisotropy induced by magnetic field amplification, coupled to the adiabatic invariance of the particle magnetic moments. Using two-dimensional (2D) particle-in-cell (PIC) simulations, we showed in Paper I that if the electron-to-ion temperature ratio is < 0.2, the ion cyclotron instability is the dominant mode for values of ion beta_i ~ 5-30 (here, beta_i is the ratio of ion thermal pressure to magnetic pressure), as appropriate for the midplane of low-luminosity accretion flows. In this work, we employ analytical theory and 1D PIC simulations (with the box aligned with the fastest growing wavevector of the ion cyclotron mode) to fully characterize how the electron heating efficiency during the growth of the ion cyclotron instability depends on the electron-to-proton temperature ratio, the plasma beta, the Alfven speed, the amplification rate of the mean field (in units of the ion Larmor frequency) and the proton-to-electron mass ratio. Our findings can be incorporated as a physically-grounded sub-grid model into global fluid simulations of low-luminosity accretion flows, thus helping to assess the validity of the two-temperature assumption.

In systems accreting well below the Eddington rate, the plasma in the innermost regions of the disk is collisionless and two-temperature, with the ions hotter than the electrons. Yet, whether a collisionless faster-than-Coulomb energy transfer mechanism exists in two-temperature accretion flows is still an open question. We study the physics of electron heating during the growth of ion velocity-space instabilities, by means of multi-dimensional particle-in-cell (PIC) simulations. A large-scale compression - embedded in a novel form of the PIC equations - continuously amplifies the field. This constantly drives a pressure anisotropy P_perp > P_parallel, due to the adiabatic invariance of the particle magnetic moments. We find that, for ion plasma beta values beta_i ~ 5-30 appropriate for the midplane of low-luminosity accretion flows, mirror modes dominate if the electron-to-proton temperature ratio is > 0.2, whereas if it is < 0.2 the ion cyclotron instability triggers the growth of strong Alfven-like waves, that pitch-angle scatter the ions to maintain marginal stability. We develop an analytical model of electron heating during the growth of the ion cyclotron instability, which we validate with PIC simulations. We find that for cold electrons (beta_e < m_e/m_i), the electron energy gain is controlled by the magnitude of the E-cross-B velocity induced by the ion cyclotron waves. This term is independent of the initial electron temperature, so it provides a solid energy floor even for electrons starting with extremely low temperatures. On the other hand, the electron energy gain for beta_e > m_e/m_i - governed by the conservation of the magnetic moment in the growing fields of the instability - is proportional to the initial electron temperature. Our results have implications for two-temperature accretion flows as well as the solar wind and intracluster plasmas. [abridged]

Electron acceleration to non-thermal energies is known to occur in low Mach number (M<5) shocks in galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Using two-dimensional (2D) particle-in-cell (PIC) plasma simulations, we showed in Paper I that electrons are efficiently accelerated in low Mach number (M=3) quasi-perpendicular shocks via a Fermi-like process. The electrons bounce between the upstream region and the shock front, with each reflection at the shock resulting in energy gain via shock drift acceleration. The upstream scattering is provided by oblique magnetic waves, that are self-generated by the electrons escaping ahead of the shock. In the present work, we employ additional 2D PIC simulations to address the nature of the upstream oblique waves. We find that the waves are generated by the shock-reflected electrons via the firehose instability, which is driven by an anisotropy in the electron velocity distribution. We systematically explore how the efficiency of wave generation and of electron acceleration depend on the magnetic field obliquity, the flow magnetization (or equivalently, the plasma beta), and the upstream electron temperature. We find that the mechanism works for shocks with high plasma beta (>20) at nearly all magnetic field obliquities, and for electron temperatures in the range relevant for galaxy clusters. Our findings offer a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

Electron acceleration to non-thermal energies in low Mach number (M<5) shocks is revealed by radio and X-ray observations of galaxy clusters and solar flares, but the electron acceleration mechanism remains poorly understood. Diffusive shock acceleration, also known as first-order Fermi acceleration, cannot be directly invoked to explain the acceleration of electrons. Rather, an additional mechanism is required to pre-accelerate the electrons from thermal to supra-thermal energies, so they can then participate in the Fermi process. In this work, we use two- and three-dimensional particle-in-cell plasma simulations to study electron acceleration in low Mach number shocks. We focus on the particle energy spectra and the acceleration mechanism in a reference run with M=3 and a quasi-perpendicular pre-shock magnetic field. We find that about 15 percent of the electrons can be efficiently accelerated, forming a non-thermal power-law tail in the energy spectrum with a slope of p~2.4. Initially, thermal electrons are energized at the shock front via shock drift acceleration. The accelerated electrons are then reflected back upstream, where their interaction with the incoming flow generates magnetic waves. In turn, the waves scatter the electrons propagating upstream back toward the shock, for further energization via shock drift acceleration. In summary, the self-generated waves allow for repeated cycles of shock drift acceleration, similarly to a sustained Fermi-like process. This mechanism offers a natural solution to the conflict between the bright radio synchrotron emission observed from the outskirts of galaxy clusters and the low electron acceleration efficiency usually expected in low Mach number shocks.

We study the dynamical interaction of the G2 cloud with the accretion flow around Sgr A* by means of three-dimensional, hydrodynamic simulations. We show the effects of the rotating accretion flow on the evolution of G2 by projecting the cloud density onto the plane of the sky, and extracting position-velocity diagrams. We study a number of possible orientations of the cloud orbit with respect to the disk. We find that once the center of mass of the cloud has crossed the pericenter, the differences between models becomes significant. Models with the cloud counter-rotating with respect to the disk are expected to reach higher blue-shifted line of sight velocities. The spatial extent of the emission depends strongly on the cloud-to-disk inclination angle. Future imaging and spectroscopy of G2 emission will shed light both on the structure of the Sgr A* disk and on the properties of the cloud.

In magnetized astrophysical outflows, the dissipation of field energy into particle energy via magnetic reconnection is often invoked to explain the observed non-thermal signatures. By means of two- and three-dimensional particle-in-cell simulations, we investigate anti-parallel reconnection in magnetically-dominated electron-positron plasmas. Our simulations extend to unprecedentedly long temporal and spatial scales, so we can capture the asymptotic state of the system beyond the initial transients, and without any artificial limitation by the boundary conditions. At late times, the reconnection layer is organized into a chain of large magnetic islands connected by thin X-lines. The plasmoid instability further fragments each X-line into a series of smaller islands, separated by X-points. At the X-points, the particles become unmagnetized and they get accelerated along the reconnection electric field. We provide definitive evidence that the late-time particle spectrum integrated over the whole reconnection region is a power-law, whose slope is harder than -2 for magnetizations sigma>10. Efficient particle acceleration to non-thermal energies is a generic by-product of the long-term evolution of relativistic reconnection in both two and three dimensions. In three dimensions, the drift-kink mode corrugates the reconnection layer at early times, but the long-term evolution is controlled by the plasmoid instability, that facilitates efficient particle acceleration, in analogy to the two-dimensional physics. Our findings have important implications for the generation of hard photon spectra in pulsar winds and relativistic astrophysical jets.

The interaction of TeV photons from blazars with the extragalactic background light produces a relativistic beam of electron-positron pairs streaming through the intergalactic medium (IGM). The fate of the beam energy is uncertain. By means of two- and three-dimensional particle-in-cell simulations, we study the non-linear evolution of dilute ultra-relativistic pair beams propagating through the IGM. We explore a wide range of beam Lorentz factors gamma_b>>1 and beam-to-plasma density ratios alpha<<1, so that our results can be extrapolated to the extreme parameters of blazar-induced beams (gamma_b~10^6 and alpha~10^(-15), for the most powerful blazars). For cold beams, we show that the oblique instability governs the early stages of evolution, but its exponential growth terminates - due to self-heating of the beam in the transverse direction - when only a negligible fraction ~(alpha/gamma_b)^(1/3)~10^(-7) of the beam energy has been transferred to the IGM plasma. Further relaxation of the beam proceeds through quasi-longitudinal modes, until the momentum dispersion in the direction of propagation saturates at DeltaP_{b,parallel}/gamma_b m_e c ~ 0.2. This corresponds to a fraction ~10% of the beam energy being ultimately transferred to the IGM plasma, irrespective of gamma_b or alpha. If the initial dispersion in beam momentum satisfies DeltaP_{b,parallel}/gamma_b m_e c > 0.2 (as typically expected for blazar-induced beams), the fraction of beam energy deposited into the IGM is much smaller than ~10%. It follows that at least ~90% of the beam energy is still available to power the GeV emission produced by inverse Compton up-scattering of the Cosmic Microwave Background by the beam pairs.

We present late-time radio and X-ray observations of the nearby sub-energetic Gamma-Ray Burst (GRB)100316D associated with supernova (SN) 2010bh. Our broad-band analysis constrains the explosion properties of GRB100316D to be intermediate between highly relativistic, collimated GRBs and the spherical, ordinary hydrogen-stripped SNe. We find that ~10^49 erg is coupled to mildly-relativistic (Gamma=1.5-2), quasi-spherical ejecta, expanding into a medium previously shaped by the progenitor mass-loss with rate Mdot ~10^-5 Msun yr^-1 (for wind velocity v_w = 1000 km s^-1). The kinetic energy profile of the ejecta argues for the presence of a central engine and identifies GRB100316D as one of the weakest central-engine driven explosions detected to date. Emission from the central engine is responsible for an excess of soft X-ray radiation which dominates over the standard afterglow at late times (t>10 days). We connect this phenomenology with the birth of the most rapidly rotating magnetars. Alternatively, accretion onto a newly formed black hole might explain the excess of radiation. However, significant departure from the standard fall-back scenario is required.

Radio observations of Gamma Ray Bursts afterglows are fundamental in providing insights into their physics and environment, and in constraining the true energetics of these sources. Nonetheless, radio observations of GRB afterglows are presently sparse in the time/frequency domain. Starting from a complete sample of 58 bright Swift long bursts (BAT6), we constructed a homogeneous sub-sample of 38 radio detections/upper limits which preserves all the properties of the parent sample. One half of the bursts have detections between 1 and 5 days after the explosion with typical fluxes F>100 muJy at 8.4 GHz. Through a Population SYnthesis Code coupled with the standard afterglow Hydrodynamical Emission model (PSYCHE) we reproduce the radio flux distribution of the radio sub-sample. Based on these results we study the detectability in the time/frequency domain of the entire long GRB population by present and future radio facilities. We find that the GRBs that typically trigger Swift can be detected at 8.4 GHz by JVLA within few days with modest exposures even at high redshifts. The final SKA can potentially observe the whole GRB population provided that there will be a dedicated GRB gamma-ray detector more sensitive than Swift. For a sizable fraction (50%) of these bursts, SKA will allow us to perform radio-calorimetry, after the trans-relativistic transition (occurring ~100 d), providing an estimate of the true (collimation corrected) energetics of GRBs.

The afterglow emission from Gamma-Ray Bursts (GRBs) is usually interpreted as synchrotron radiation from relativistic electrons accelerated at the GRB external shock, that decelerates from ultra-relativistic to non-relativistic speeds as it sweeps up the surrounding medium. We investigate the temporal decay of the emission from GRB afterglows at late times, when the bulk of the shock-accelerated electrons are non-relativistic. For a uniform circumburst medium, we show that such "deep Newtonian phase" begins at t_{DN} ~ 3 epsilon_{e,-1}^{5/6} t_{ST}, where t_{ST} marks the transition of the blast wave to the non-relativistic spherically-symmetric Sedov-Taylor solution, and epsilon_e = 0.1 epsilon_{e,-1} quantifies the amount of shock energy transferred to the post-shock electrons. For typical parameters, the deep Newtonian stage starts ~ 0.5-several years after the GRB. The radio flux in this phase decays as F_nu ~ t^{-3(p+1)/10}, with decay slope between -0.9 and -1.2, for a power-law distribution of shock-accelerated electrons with index 2 < p < 3. This is shallower than the commonly assumed scaling F_nu ~ t^{-3(5p-7)/10} (with decay slope between -0.9 and -2.4) derived by Frail et al. (2000), which only applies if the GRB shock is non-relativistic, but the electron distribution still peaks at ultra-relativistic energies (a regime that is relevant for a narrow time interval, and only if t_{DN} > t_{ST}, or epsilon_e > 0.03). We discuss how the deep Newtonian phase can be reliably used for GRB calorimetry, and we comment on the good detection prospects of trans-relativistic blast waves at 0.1-10 GHz with EVLA and LOFAR.

Forward shocks caused by the interaction between a relativistic blast wave and the circum-burst medium are thought to be responsible for the afterglow emission in Gamma-Ray Bursts (GRBs). We consider the hydrodynamics of a spherical relativistic blast wave expanding into the surrounding medium and we generalize the standard theory in order to account for several effects that are generally ignored. In particular, we consider the role of adiabatic and radiative losses on the hydrodynamical evolution of the shock, under the assumption that the cooling losses are fast. Our model can describe adiabatic, fully radiative and semi-radiative blast waves, and can describe the effects of a time-varying radiative efficiency. The equations we present are valid for arbitrary density profiles, and also for a circum-burst medium enriched with electron-positron pairs. The presence of pairs enhances the fraction of shock energy gained by the leptons, thus increasing the importance of radiative losses. Our model allows to study whether the high-energy (>0.1 GeV) emission in GRBs may originate from afterglow radiation. In particular, it is suitable to test whether the fast decay of the high-energy light curve observed in several Fermi LAT GRBs can be ascribed to an initial radiative phase, followed by the standard adiabatic evolution.

We perform detailed magnetohydrodynamic simulations of the gas cloud G2 interacting with the accretion flow around the Galactic Center black hole Sgr A*. We take as our initial conditions a steady-state, converged solution of the accretion flow obtained earlier using the general-relativistic magnetohydrodynamic code HARM. Using the observed parameters for the cloud's orbit, we compute the interaction of the cloud with the ambient gas and identify the shock structure that forms ahead of the cloud. We show that for many configurations, the cloud front crosses orbit pericenter 7 to 9 months earlier than the center-of-mass.

We calculate radio light curves produced by the bow shock that is likely to form in front of the G2 cloud when it penetrates the accretion disk of Sgr A*. The shock acceleration of the radio-emitting electrons is captured self-consistently by means of first-principles particle-in-cell simulations. We show that the radio luminosity is expected to reach maximum in early 2013, roughly a month after the bow shock crosses the orbit pericenter. We estimate the peak radio flux at 1.4 GHz to be 1.4 - 22 Jy depending on the assumed orbit orientation and parameters. We show that the most promising frequencies for radio observations are in the 0.1<nu<1 GHz range, for which the bow shock emission will be much stronger than the intrinsic radio flux for all the models considered.

How accretion proceeds around the massive black hole in the Galactic center and other highly sub-Eddington accretors remains poorly understood. The orbit of the S2 star in the Galactic center passes through the accretion disk of the massive black hole and any observational signature from such interaction may be used as an accretion probe. Because of its early stellar type, S2 is expected to possess a fairly powerful wind. We show here that the ram pressure of the accretion disk shocks the stellar wind fairly close to the star. The shocked fluid reaches a temperature of ~ 1 keV and cools efficiently through optically thin, thermal bremsstrahlung emission. The radiation from the shocked wind peaks around the epoch of the pericenter passage of the star at a luminosity potentially comparable to the quiescent emission detected from Sgr A*. Detection of shocked wind radiation can constrain the density of the accretion disk at a distance of several thousands of gravitational radii from the black hole.

Bright star-forming galaxies radiate well below their Eddington Limit. The value of the flux-mean opacity that mediates the radiation force onto matter is orders of magnitude smaller than the UV or optical dust opacity. On empirical grounds, it is shown that high-redshift ULIRGs radiate at two orders of magnitude below their Eddington Limit, while the local starbursters M82 and Arp 220 radiate at a few percent of their Eddington Limit. A simple model for the radiative transfer of UV and optical light in dust-rich environments is considered. Radiation pressure on dust does not greatly affect the large-scale gas dynamics of star-forming galaxies.

The afterglow emission from gamma-ray bursts (GRBs) is usually interpreted as synchrotron radiation from electrons accelerated at the GRB external shock, that propagates with relativistic velocities into the magnetized interstellar medium. By means of multi-dimensional particle-in-cell simulations, we investigate the acceleration performance of weakly magnetized relativistic shocks, in the magnetization range 0<sigma<1e-1. The pre-shock magnetic field is orthogonal to the flow, as generically expected for relativistic shocks. We find that relativistic perpendicular shocks propagating in electron-positron plasmas are efficient particle accelerators if the magnetization is sigma<1e-3. For electron-ion plasmas, the transition to efficient acceleration occurs for sigma<3e-5. Here, the acceleration process proceeds similarly for the two species, since the electrons enter the shock nearly in equipartition with the ions, as a result of strong pre-heating in the self-generated upstream turbulence. In both electron-positron and electron-ion shocks, we find that the maximum energy of the accelerated particles scales in time as t^(1/2). This scaling is shallower than the so-called (and commonly assumed) Bohm limit, and it naturally results from the small-scale nature of the Weibel turbulence generated in the shock layer. In magnetized plasmas, the energy of the accelerated particles increases until it reaches a saturation value that scales with the magnetization as sigma^(-1/4). Further energization is prevented by the fact that the self-generated turbulence is confined within a finite region of thickness proportional to sigma^(-1/2) around the shock. Our results can provide physically-grounded inputs for models of non-thermal emission from a variety of astrophysical sources, with particular relevance to GRB afterglows.

A dense ionized cloud of gas has been recently discovered to be moving directly toward the supermassive black hole, Sgr A*, at the Galactic Center. In June 2013, at the pericenter of its highly eccentric orbit, the cloud will be approximately 3100 Schwarzschild radii from the black hole and will move supersonically through the ambient hot gas with a velocity of v_p ~ 5400 km/s. A bow shock is likely to form in front of the cloud and could accelerate electrons to relativistic energies. We estimate via particle-in-cell simulations the energy distribution of the accelerated electrons and show that the non-thermal synchrotron emission from these electrons might exceed the quiescent radio emission from Sgr A* by a factor of several. The enhanced radio emission should be detectable at GHz and higher frequencies around the time of pericentric passage and in the following months. The bow shock emission is expected to be displaced from the quiescent radio emission of Sgr A* by ~33 mas. Interferometric observations could resolve potential changes in the radio image of Sgr A* at wavelengths < 6 cm.

By means of two- and three-dimensional particle-in-cell simulations, we investigate the process of driven magnetic reconnection at the termination shock of relativistic striped flows. In pulsar winds and in magnetar-powered relativistic jets, the flow consists of stripes of alternating magnetic field polarity, separated by current sheets of hot plasma. At the wind termination shock, the flow compresses and the alternating fields annihilate by driven magnetic reconnection. Irrespective of the stripe wavelength "lambda" or the wind magnetization "sigma" (in the regime sigma>>1 of magnetically-dominated flows), shock-driven reconnection transfers all the magnetic energy of alternating fields to the particles, whose average Lorentz factor increases by a factor of sigma with respect to the pre-shock value. In the limit lambda/(r_L*sigma)>>1, where r_L is the relativistic Larmor radius in the wind, the post-shock particle spectrum approaches a flat power-law tail with slope around -1.5, populated by particles accelerated by the reconnection electric field. The presence of a current-aligned "guide" magnetic field suppresses the acceleration of particles only when the guide field is stronger than the alternating component. Our findings place important constraints on the models of non-thermal radiation from Pulsar Wind Nebulae and relativistic jets.

The relativistic wind of obliquely-rotating pulsars consists of toroidal stripes of opposite magnetic field polarity, separated by current sheets of hot plasma. By means of two- and three-dimensional particle-in-cell simulations, we investigate particle acceleration and magnetic field dissipation at the termination shock of a relativistic striped wind. At the shock, the flow compresses and the alternating fields annihilate by driven magnetic reconnection. Irrespective of the stripe wavelength "lambda" or the wind magnetization "sigma" (in the regime sigma>>1 of magnetically-dominated flows), shock-driven reconnection transfers all the magnetic energy of alternating fields to the particles, whose average Lorentz factor increases by a factor of sigma with respect to the pre-shock value. The shape of the post-shock spectrum depends primarily on the ratio lambda/(r_L*sigma), where r_L is the relativistic Larmor radius in the wind. The spectrum becomes broader as the value of lambda/(r_L*sigma) increases, passing from a relativistic Maxwellian to a flat power-law tail with slope around -1.5, populated by particles accelerated by the reconnection electric field. Close to the equatorial plane of the wind, where the stripes are symmetric, the highest energy particles resulting from magnetic reconnection can escape ahead of the shock, and be injected into a Fermi-like acceleration process. In the post-shock spectrum, they populate a power-law tail with slope around -2.5, that extends beyond the flat component produced by reconnection. Our study suggests that the spectral break between the radio and the optical band in Pulsar Wind Nebulae can be a natural consequence of particle acceleration at the termination shock of striped pulsar winds.

We investigate shock structure and particle acceleration in relativistic magnetized collisionless electron-ion shocks by means of 2.5D particle-in-cell simulations with ion-to-electron mass ratios (m_i/m_e) ranging from 16 to 1000. We explore a range of inclination angles between the pre-shock magnetic field and the shock normal. In "subluminal" shocks, where relativistic particles can escape ahead of the shock along the magnetic field lines, ions are efficiently accelerated via a Fermi-like mechanism. The downstream ion spectrum consists of a relativistic Maxwellian and a high-energy power-law tail, which contains ~5% of ions and ~30% of ion energy. Its slope is -2.1. Upstream electrons enter the shock with lower energy than ions, so they are more strongly tied to the field. As a result, only ~1% of the incoming electrons are Fermi-accelerated at the shock before being advected downstream, where they populate a steep power-law tail (with slope -3.5). For "superluminal" shocks, where relativistic particles cannot outrun the shock along the field, the self-generated turbulence is not strong enough to permit efficient Fermi acceleration, and the ion and electron downstream spectra are consistent with thermal distributions. The incoming electrons are heated up to equipartition with ions, due to strong electromagnetic waves emitted by the shock into the upstream. Thus, efficient electron heating (>15% of the upstream ion energy) is the universal property of relativistic electron-ion shocks, but significant nonthermal acceleration of electrons (>2% by number, >10% by energy, with slope flatter than -2.5) is hard to achieve in magnetized flows and requires weakly magnetized shocks (magnetization <1e-3). These findings place important constraints on the models of AGN jets and Gamma Ray Bursts that invoke particle acceleration in relativistic magnetized electron-ion shocks.

We present a feedback mechanism for supermassive black holes and their host bulges that operates during epochs of radio-loud quasar activity. In the radio cores of relativistic quasar jets, internal shocks convert a fraction of ordered bulk kinetic energy into randomized relativistic ions, or in other words cosmic rays. By employing a phenomenologically-motivated jet model, we show that enough 1-100 GeV cosmic rays escape the radio core into the host galaxy to break the Eddington limit in cosmic rays. As a result, hydrostatic balance is lost and a cosmic ray momentum-driven wind develops, expelling gas from the host galaxy and thus self-limiting the black hole and bulge growth. Although the interstellar cosmic ray power is much smaller than the quasar photon luminosity, cosmic rays provide a stronger feedback than UV photons, since they exchange momentum with the galactic gas much more efficiently. The amount of energy released into the host galaxy as cosmic rays, per unit of black hole rest mass energy, is independent of black hole mass. It follows that radio-loud jets should be more prevalent in relatively massive systems since they sit in galaxies with relatively deep gravitational potentials. Therefore, jet-powered cosmic ray feedback not only self-regulates the black hole and bulge growth, but also provides an explanation for the lack of radio-loud activity in relatively small galaxies. By employing basic known facts regarding the physical conditions in radio cores, we approximately reproduce both the slope and the normalization of the M_{BH}-M_{Bulge} relation.

We extract synthetic photon spectra from first-principles particle-in-cell simulations of relativistic shocks propagating in unmagnetized pair plasmas. The two basic ingredients for the radiation, namely accelerated particles and magnetic fields, are produced self-consistently as part of the shock evolution. We use the method of Hededal & Nordlund (2005) and compute the photon spectrum via Fourier transform of the electric far-field from a large number of particles, sampled directly from the simulation. We find that the spectrum from relativistic collisionless shocks is entirely consistent with synchrotron radiation in the magnetic fields generated by Weibel instability. We can recover the so-called "jitter'' regime only if we artificially reduce the strength of the electromagnetic fields, such that the wiggler parameter K = qB lambda/mc^2 becomes much smaller than unity ("B" and "lambda" are the strength and scale of the magnetic turbulence, respectively). These findings may place constraints on the origin of non-thermal emission in astrophysics, especially for the interpretation of the hard (harder than synchrotron) low-frequency spectrum of Gamma-Ray Bursts.

We study the polarized infrared emission by Polycyclic Aromatic Hydrocarbons (PAHs), when anisotropically illuminated by UV photons. PAH molecules are modeled as planar disks with in-plane and out-of-plane vibrational dipoles. As first pointed out by Leger (1988), infrared emission features resulting from in-plane and out-of-plane modes should have orthogonal polarization directions. We show analytically how the degree of polarization depends on the viewing geometry and the molecule's internal alignment between principal axis of inertia and angular momentum, which gets worse after photon absorption. Longer wavelength features, emitted after better internal alignment is recovered, should be more strongly polarized. The degree of polarization for uni-directional illumination (e.g., by a star) is larger than for diffuse illumination (e.g., by a disk galaxy), all else being equal. For PAHs in the Cold Neutral Medium, the predicted polarization is probably too small to distinguish from the contribution of linear dichroism by aligned foreground dust. The level of polarization predicted for PAH emission from the Orion Bar is only ~0.06% at 3.3 microns; Sellgren et al. (1988) report a much larger value, 0.86+-0.28%, which suggests that the smallest PAHs may have moderately suprathermal rotation rates. Future observations of (or upper limits on) the degree of polarization for the Orion Bar or for dust above edge-on galaxies (e.g., NGC 891 or M82) may constrain the internal alignment of emitting PAHs, thus providing clues to their rotational dynamics.

We investigate shock structure and particle acceleration in relativistic magnetized collisionless pair shocks by means of 2.5D and 3D particle-in-cell simulations. We explore a range of inclination angles between the pre-shock magnetic field and the shock normal. We find that only magnetic inclinations corresponding to "subluminal" shocks, where relativistic particles following the magnetic field can escape ahead of the shock, lead to particle acceleration. The downstream spectrum in such shocks consists of a relativistic Maxwellian and a high-energy power-law tail with exponential cutoff. For increasing magnetic inclination in the subluminal range, the high-energy tail accounts for an increasing fraction of particles (from ~1% to ~2%) and energy (from ~4% to ~12%). The spectral index of the power law increases with angle from -2.8+-0.1 to -2.3+-0.1. Particle energization is driven by the Diffusive Shock Acceleration process for nearly parallel shocks, and switches to Shock-Drift Acceleration for larger subluminal inclinations. For "superluminal" shocks, the downstream particle spectrum does not show any significant suprathermal tail. As seen from the upstream frame, efficient acceleration in relativistic (Lorentz factor gamma0 > 5) magnetized (sigma > 0.03) flows exists only for a very small range of magnetic inclination angles (< 34/gamma0 degrees), so relativistic astrophysical pair shocks have to be either nearly parallel or weakly magnetized to generate nonthermal particles. These findings place constraints on the models of AGN jets, Pulsar Wind Nebulae and Gamma Ray Bursts that invoke particle acceleration in relativistic magnetized shocks. (Abridged)

Afterglows of gamma-ray bursts are believed to require magnetic fields much stronger than that of the compressed pre-shock medium. As an alternative to microscopic plasma instabilities, we propose amplification of the field by macroscopic turbulence excited by the interaction of the shock with a clumpy pre-shock medium, for example a stellar wind. Using a recently developed formalism for localized perturbations to an ultra-relativistic shock, we derive constraints on the lengthscale, amplitude, and volume filling factor of density clumps required to produce a given magnetic energy fraction within the expansion time of the shock, assuming that the energy in the field achieves equipartion with the turbulence. Stronger and smaller-scale inhomogeneities are required for larger shock Lorentz factors. Hence it is likely that the magnetic energy fraction evolves as the shock slows. This could be detected by monitoring the synchrotron cooling frequency if the radial density profile ahead of the shock, smoothed over clumps, is known.