Fully kinetic two-dimensional particle-in-cell simulations are used to study electron acceleration at high-Mach-number nonrelativistic perpendicular shocks. SNR shocks are mediated by the Weibel instability which is excited because of an interaction between shock-reflected and upstream ions. Nonlinear evolution of the Weibel instability leads to the formation of current sheets. At the turbulent shock ramp the current sheets decay through magnetic reconnection. The number of reconnection sites strongly depends on the ion-to-electron mass ratio and the Alfv\'enic Mach number of the simulated shock. Electron acceleration is observed at locations where magnetic reconnection operates. For the highest mass ratios almost all electrons are involved in magnetic reconnection, which makes the magnetic reconnection the dominant acceleration process for electrons at these shocks. We discuss the relevance of our results for 3D systems with realistic ion-to-electron mass ratio.

We explore electron pre-acceleration at high Mach-number nonrelativistic perpendicular shocks at, e.g., young supernova remnants, which are a prerequisite of further acceleration to very high energies via diffusive shock acceleration. Using fully kinetic particle-in-cell simulations of shocks and electron dynamics in them, we investigate the influence of shock-surfing acceleration at the shock foot on the nonthermal population of electrons downstream of the shock. The shock-surfing acceleration is followed by further energization at the shock ramp where the Weibel instability spawns a type of second-order Fermi acceleration. The combination of these two processes leads to the formation of a nonthermal electron population, but the importance of shock-surfing acceleration becomes smaller for larger ion-to-electron mass ratio in the simulation. We discuss the resulting electron spectra and the relevance of our results to the physics of systems with real ion-to-electron mass ratio and fully three-dimensional behavior.

We revisit the effect of non-linear Landau (NL) damping on the electrostatic instability of blazar-induced pair beams, using a realistic pair-beam distribution. We employ a simplified 2D model in ${\bf k}$-space to study the evolution of the electric-field spectrum and to calculate the relaxation time of the beam. We demonstrate that the 2D model is an adequate representation of the 3D physics. We find that non-linear Landau damping, once it operates efficiently, transports essentially the entire wave energy to small wavenumbers where wave driving is weak or absent. The relaxation time also strongly depends on the IGM temperature, $T_\mathrm{IGM}$, and for $T_\mathrm{IGM}\ll10$ eV, and in the absence of any other damping mechanism, the relaxation time of the pair beam is longer than the inverse Compton (IC) scattering time. The weak late-time beam energy losses arise from the accumulation of wave energy at small $k$, that non-linearly drains the wave energy at the resonant $\mathbf{k}$ of the pair-beam instability. Any other dissipation process operating at small $k$ would reduce that wave-energy drain and hence lead to stronger pair-beam energy losses. As an example, collisions reduce the relaxation time by an order of magnitude, although their rate is very small. Other non-linear processes, such as the modulation instability, could provide additional damping of the non-resonant waves and dramatically reduce the relaxation time of the pair beam. An accurate description of the spectral evolution of the electrostatic waves is crucial for calculating the relaxation time of the pair beam.

Electron injection at high Mach-number nonrelativistic perpendicular shocks is studied here for parameters that are applicable to young SNR shocks. Using high-resolution large-scale two-dimensional fully kinetic particle-in-cell (PIC) simulations and tracing individual particles we in detail analyze the shock surfing acceleration (SSA) of electrons at the leading edge of the shock foot. The central question is to what degree the process can be captured in 2D3V simulations. We find that the energy gain in SSA always arises from the electrostatic field of a Buneman wave. Electron energization is more efficient in the out-of-plane orientation of the large-scale magnetic field because both the phase speed and the amplitude of the waves are higher than for the in-plane scenario. Also, a larger number of electrons is trapped by the waves compared to the in-plane configuration. We conclude that significant modifications of the simulation parameters are needed to reach the same level of SSA efficiency as in simulations with out-of-plane magnetic field or 3D simulations.

We perform particle-in-cell simulations of perpendicular nonrelativistic collisionless shocks to study electron heating and pre-acceleration for parameters that permit extrapolation to the conditions at young supernova remnants. Our high-resolution large-scale numerical experiments sample a representative portion of the shock surface and demonstrate that the efficiency of electron injection is strongly modulated with the phase of the shock reformation. For plasmas with low and moderate temperature (plasma beta $\beta_{\rm p}=5\cdot 10^{-4}$ and $\beta_{\rm p}=0.5$), we explore the nonlinear shock structure and electron pre-acceleration for various orientations of the large-scale magnetic field with respect to the simulation plane while keeping it at $90^\circ$ to the shock normal. Ion reflection off the shock leads to the formation of magnetic filaments in the shock ramp, resulting from Weibel-type instabilities, and electrostatic Buneman modes in the shock foot. In all cases under study, the latter provides first-stage electron energization through the shock-surfing acceleration (SSA) mechanism. The subsequent energization strongly depends on the field orientation and proceeds through adiabatic or second-order Fermi acceleration processes for configurations with the out-of-plane and in-plane field components, respectively. For strictly out-of-plane field the fraction of supra-thermal electrons is much higher than for other configurations, because only in this case the Buneman modes are fully captured by the 2D simulation grid. Shocks in plasma with moderate $\beta_{\rm p}$ provide more efficient pre-acceleration. The relevance of our results to the physics of fully three-dimensional systems is discussed.

We present a detailed analysis of the best-quality multi-wavelength data gathered for the large-scale jet in the core-dominated quasar 3C 273. We analyze all the archival observations of the target with the Chandra X-ray Observatory, the far-ultraviolet observations with the Hubble Space Telescope, and the 8.4 GHz map obtained with the Very Large Array. In our study we focus on investigating the morphology of the outflow at different frequencies, and therefore we apply various techniques for the image deconvolution, paying particular attention to a precise modeling of the Chandra and Hubble point spread functions. We find that the prominent brightness enhancements in the X-ray and far-ultraviolet jet of 3C 273 - the "knots" - are not point-like, and can be resolved transversely as extended features with sizes of about $\simeq 0.5$ kpc. Also, the radio outflow is wider than the deconvolved X-ray/ultraviolet jet. We have also found circumstantial evidence that the intensity peaks of the X-ray knots are located systematically upstream of the corresponding radio intensity peaks, with the projected spatial offsets along the jet ranging from $\lesssim 0.2$ kpc up to $\simeq 1$ kpc. We discuss our findings in the wider context of multi-component models for the emission and structure of large-scale quasar jets, and speculate on the physical processes enabling an efficient acceleration of the emitting ultra-relativistic electrons along the entire jet length that exceeds 100 kpc.

A nonresonant cosmic-ray-current-driven instability may operate in the shock precursors of young supernova remnants and be responsible for magnetic-field amplification, plasma heating and turbulence. Earlier simulations demonstrated magnetic-field amplification, and in kinetic studies a reduction of the relative drift between cosmic rays and thermal plasma was observed as backreaction. However, all published simulations used periodic boundary conditions, which do not account for mass conservation in decelerating flows and only allow the temporal development to be studied. Here we report results of fully kinetic Particle-In-Cell simulations with open boundaries that permit inflow of plasma on one side of the simulation box and outflow at the other end, hence allowing an investigation of both the temporal and the spatial development of the instability. Magnetic-field amplification proceeds as in studies with periodic boundaries and, observed here for the first time, the reduction of relative drifts causes the formation of a shock-like compression structure at which a fraction of the plasma ions are reflected. Turbulent electric field generated by the nonresonant instability inelastically scatters cosmic rays, modifying and anisotropizing their energy distribution. Spatial CR scattering is compatible with Bohm diffusion. Electromagnetic turbulence leads to significant nonadiabatic heating of the background plasma maintaining bulk equipartition between ions and electrons. The highest temperatures are reached at sites of large-amplitude electrostatic fields. Ion spectra show supra-thermal tails resulting from stochastic scattering in the turbulent electric field. Together, these modifications in the plasma flow will affect the properties of the shock and particle acceleration there.