Highly energetic, relativistic electrons are commonly present in many astrophysical systems, from solar flares to the intra-cluster medium, as indicated by observed electromagnetic radiation. However, open questions remain about the mechanisms responsible for their acceleration, and possible re-acceleration. Ubiquitous plasma turbulence is one of the possible universal mechanisms. We study the energization of transrelativistic electrons in turbulence using hybrid particle-in-cell, which provide a realistic model of Alfv\'{e}nic turbulence from MHD to sub-ion scales, and test particle simulations for electrons. We find that, depending on the electron initial energy and turbulence strength, electrons may undergo a fast and efficient phase of energization due to the magnetic curvature drift during the time they are trapped in dynamic magnetic structures. In addition, electrons are accelerated stochastically which is a slower process that yields lower maximum energies. The combined effect of these two processes determines the overall electron acceleration. With appropriate turbulence parameters, we find that superthermal electrons can be accelerated up to relativistic energies. For example, with heliospheric parameters and a relatively high turbulence level, rapid energization to MeV energies is possible.

The first two orbits of the Parker Solar Probe (PSP) spacecraft have enabled the first in situ measurements of the solar wind down to a heliocentric distance of 0.17 au (or 36 Rs). Here, we present an analysis of this data to study solar wind turbulence at 0.17 au and its evolution out to 1 au. While many features remain similar, key differences at 0.17 au include: increased turbulence energy levels by more than an order of magnitude, a magnetic field spectral index of -3/2 matching that of the velocity and both Elsasser fields, a lower magnetic compressibility consistent with a smaller slow-mode kinetic energy fraction, and a much smaller outer scale that has had time for substantial nonlinear processing. There is also an overall increase in the dominance of outward-propagating Alfv\'enic fluctuations compared to inward-propagating ones, and the radial variation of the inward component is consistent with its generation by reflection from the large-scale gradient in Alfv\'en speed. The energy flux in this turbulence at 0.17 au was found to be ~10% of that in the bulk solar wind kinetic energy, becoming ~40% when extrapolated to the Alfv\'en point, and both the fraction and rate of increase of this flux towards the Sun is consistent with turbulence-driven models in which the solar wind is powered by this flux.

Magnetospheric Multiscale (MMS) observations of plasma turbulence generated by a Kelvin-Helmholtz (KH) event at the Earth's magnetopause are compared with a high-resolution two-dimensional (2D) hybrid direct numerical simulation (DNS) of decaying plasma turbulence driven by large-scale balanced Alfv\'enic fluctuations. The simulation, set up with four observation-driven physical parameters (ion and electron betas, turbulence strength, and injection scale) exhibits a quantitative agreement on the spectral, intermittency, and cascade-rate properties with in situ observations, despite the different driving mechanisms. Such agreement demonstrates a certain universality of the turbulent cascade from magnetohydrodynamic (MHD) to sub-ion scales, whose properties are mainly determined by the selected parameters, also indicating that the KH instability-driven turbulence has a quasi-2D nature. The validity of the Taylor hypothesis in the sub-ion spatial range suggests that the fluctuations at sub-ion scales have predominantly low frequencies, consistent with a kinetic Alfv\'en wave-like nature or with quasi-static structures. Finally, the third-order structure function analysis indicates that the cascade rate of the turbulence generated by a KH event in the magnetopause is an order of magnitude larger than in the ambient magnetosheath.

The statistical properties of ions in two-dimensional fully developed turbulence have been compared between two different numerical algorithms. In particular, we compare Hybrid Particle In Cell (hybrid PIC with fluid electrons) and full PIC simulations, focusing on particle diffusion and acceleration phenomena. To investigate several heliospheric plasma conditions, a series of numerical simulations has been performed by varying the plasma $\beta$ - the ratio between kinetic and magnetic pressure. These numerical studies allow the exploration of different scenarios, going from the solar corona (low $\beta$) to the solar wind ($\beta \sim 1$), as well as the Earth's magnetosheath (high $\beta$). It has been found that the two approaches compare pretty well, especially for the spectral properties of the magnetic field and the ion diffusion statistics. Small differences among the models have been found regarding the electric field behaviour at sub-ion scales and the acceleration statistics, due evidently to the more consistent treatment of the plasma in the full PIC approach.

A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics.

We present the first observational evidence of the irregular surface of interplanetary (IP) shocks by using multi-spacecraft observations of the Cluster mission. In total we discuss observations of four IP shocks that exhibit moderate Alfv\'enic Mach numbers (M$_A\leq$6.5). Three of them are high-$\beta$ shocks with upstream $\beta$ = 2.2--3.7. During the times when these shocks were observed, the Cluster spacecraft formed constellations with inter-spacecraft separations ranging from less than one upstream ion inertial length (d$_i$) up to 100~d$_i$. Expressed in kilometers, the distances ranged between 38~km and $\sim$10$^4$~km. We show that magnetic field profiles and the local shock normals of observed shocks are very similar when the spacecraft are of the order of one d$_i$ apart, but are strikingly different when the distances increase to ten or more d$_i$. We interpret these differences to be due to the irregular surface of IP shocks and discuss possible causes for such irregularity. We strengthen our interpretation by comparing observed shock profiles with profiles of simulated shocks. The latter had similar characteristics (M$_A$, $\theta_{BN}$, upstream ion $\beta$) as observed shocks and the profiles were obtained at separations across the simulation domain equivalent to the Cluster inter-spacecraft distances.

Shock accelerated electrons are found in many astrophysical environments, and the mechanisms by which they are accelerated to high energies are still not completely clear. For relatively high Mach numbers, the shock is supercritical, and its front exhibit broadband fluctuations, or ripples. Shock surface fluctuations have been object of many observational and theoretical studies, and are known to be important for electron acceleration. We employ a combination of hybrid Particle-In-Cell and test-particle methods to study how shock surface fluctuations influence the acceleration of suprathermal electrons in fully three dimensional simulations, and we give a complete comparison for the 2D and 3D cases. A range of different quasi-perpendicular shocks in 2D and 3D is examined, over a range of parameters compatible with the ones observed in the solar wind. Initial electron velocity distributions are taken as kappa functions, consistent with solar wind \emph{in-situ} measurements. Electron acceleration is found to be enhanced in the supercritical regime compared to subcritical. When the fully three-dimensional structure of the shock front is resolved, slightly larger energisation for the electrons is observed, and we suggest that this is due to the possibility for the electrons to interact with more than one surface fluctuation per interaction. In the supecritical regime, efficient electron energisation is found also at shock geometries departing from $\theta_{Bn}$ very close to 90$^\circ$. Two dimensional simulations show indications of unrealistic electron trapping, leading to slightly higher energisation in the subcritical cases.

A leading explanation for the origin of Galactic cosmic rays is acceleration at high-Mach number shock waves in the collisionless plasma surrounding young supernova remnants. Evidence for this is provided by multi-wavelength non-thermal emission thought to be associated with ultrarelativistic electrons at these shocks. However, the dependence of the electron acceleration process on the orientation of the upstream magnetic field with respect to the local normal to the shock front (quasi-parallel/quasi-perpendicular) is debated. Cassini spacecraft observations at Saturn's bow shock has revealed examples of electron acceleration under quasi-perpendicular conditions, and the first in situ evidence of electron acceleration at a quasi-parallel shock. Here we use Cassini data to make the first comparison between energy spectra of locally accelerated electrons under these differing upstream magnetic field regimes. We present data taken during a quasi-perpendicular shock crossing on 2008 March 8 and during a quasi-parallel shock crossing on 2007 February 3, highlighting that both were associated with electron acceleration to at least MeV energies. The magnetic signature of the quasi-perpendicular crossing has a relatively sharp upstream-downstream transition, and energetic electrons were detected close to the transition and immediately downstream. The magnetic transition at the quasi-parallel crossing is less clear, energetic electrons were encountered upstream and downstream, and the electron energy spectrum is harder above ~100 keV. We discuss whether the acceleration is consistent with diffusive shock acceleration theory in each case, and suggest that the quasi-parallel spectral break is due to an energy-dependent interaction between the electrons and short, large-amplitude magnetic structures.

In the outer heliosphere, beyond the solar wind termination shock, it is expected that the warped heliospheric current sheet forms a region of closely-packed, multiple, thin current sheets. Such a system may be subject to the ion-kinetic tearing instability, and hence generate magnetic islands and hot populations of ions associated with magnetic reconnection. Reconnection processes in this environment have important implications for local particle transport, and for particle acceleration at reconnection sites and in turbulence. We study this complex environment by means of three-dimensional hybrid simulations over long time scales, in order to capture the evolution from linear growth of the tearing instability to a fully developed turbulent state at late times. The final state develops from the highly ordered initial state via both forward and inverse cascades. Component and spectral anisotropy in the magnetic fluctuations is present when a guide field is included. The inclusion of a population of new-born interstellar pickup protons does not strongly affect these results. Finally, we conclude that reconnection between multiple current sheets can act as an important source of turbulence in the outer heliosphere, with implications for energetic particle acceleration and propagation.

Collisionless shocks are efficient particle accelerators. At Earth, ions with energies exceeding 100 keV are seen upstream of the bow shock when the magnetic geometry is quasi-parallel, and large-scale supernova remnant shocks can accelerate ions into cosmic rays energies. This energization is attributed to diffusive shock acceleration, however, for this process to become active the ions must first be sufficiently energized. How and where this initial acceleration takes place has been one of the key unresolved issues in shock acceleration theory. Using Cluster spacecraft observations, we study the signatures of ion reflection events in the turbulent transition layer upstream of the terrestrial bow shock, and with the support of a hybrid simulation of the shock, we show that these reflection signatures are characteristic of the first step in the ion injection process. These reflection events develop in particular in the region where the trailing edge of large-amplitude upstream waves intercept the local shock ramp and the upstream magnetic field changes from quasi-perpendicular to quasi-parallel. The dispersed ion velocity signature observed can be attributed to a rapid succession of ion reflections at this wave boundary. After the ions' initial interaction with the shock, they flow upstream along the quasi-parallel magnetic field. Each subsequent wave front in the upstream region will sweep the ions back toward the shock, where they gain energy with each transition between the upstream and the shock wave frames. Within three to five gyroperiods, some ions have gained enough parallel velocity to escape upstream, thus completing the injection process.

Supercritical collisionless perpendicular shocks have an average macrostructure determined primarily by the dynamics of ions specularly reflected at the magnetic ramp. Within the overall macrostructure, instabilities, both linear and nonlinear, generate fluctuations and microstructure. To identify the sources of such microstructure, high-resolution two- and three-dimensional simulations have been carried out using the hybrid method, wherein the ions are treated as particles and the electron response is modelled as a massless fluid. We confirm the results of earlier 2-D simulations showing both field-parallel aligned propagating fluctuations and fluctuations carried by the reflected-gyrating ions. In addition, it is shown that, for 2-D simulations of the shock coplanarity plane, the presence of short-wavelength fluctuations in all magnetic components is associated with the ion Weibel instability driven at the upstream edge of the foot by the reflected-gyrating ions. In 3-D simulations we show for the first time that the dominant microstructure is due to a coupling between field-parallel propagating fluctuations in the ramp and the motion of the reflected ions. This results in a pattern of fluctuations counter-propagating across the surface of the shock at an angle inclined to the magnetic field direction, due to a combination of field-parallel motion at the Alfv\'en speed of the ramp, and motion in the sense of gyration of the reflected ions.

Shock waves exist throughout the universe and are fundamental to understanding the nature of collisionless plasmas. Reformation is a process, driven by microphysics, which typically occurs at high Mach number supercritical shocks. While ongoing studies have investigated this process extensively both theoretically and via simulations, their observations remain few and far between. In this letter we present a study of very high Mach number shocks in a parameter space that has been poorly explored and we identify reformation using in situ magnetic field observations from the Cassini spacecraft at 10 AU. This has given us an insight into quasi-perpendicular shocks across two orders of magnitude in Alfven Mach number (MA) which could potentially bridge the gap between modest terrestrial shocks and more exotic astrophysical shocks. For the first time, we show evidence for cyclic reformation controlled by specular ion reflection occurring at the predicted timescale of ~0.3 {\tau}c, where {\tau}c is the ion gyroperiod. In addition, we experimentally reveal the relationship between reformation and MA and focus on the magnetic structure of such shocks to further show that for the same MA, a reforming shock exhibits stronger magnetic field amplification than a shock that is not reforming.

The solar wind magnetic field contains rotations at a broad range of scales, which have been extensively studied in the MHD range. Here we present an extension of this analysis to the range between ion and electron kinetic scales. The distribution of rotation angles was found to be approximately log-normal, shifting to smaller angles at smaller scales almost self-similarly, but with small, statistically significant changes of shape. The fraction of energy in fluctuations with angles larger than $\alpha$ was found to drop approximately exponentially with $\alpha$, with e-folding angle $9.8^\circ$ at ion scales and $0.66^\circ$ at electron scales, showing that large angles ($\alpha > 30^\circ$) do not contain a significant amount of energy at kinetic scales. Implications for kinetic turbulence theory and the dissipation of solar wind turbulence are discussed.

We perform fully-kinetic particle-in-cell simulations of an hot plasma that expands radially in a cylindrical geometry. The aim of the paper is to study the consequent development of the electron temperature anisotropy in an expanding plasma flow as found in a collisionless stellar wind. Kinetic plasma theory and simulations have shown that the electron temperature anisotropy is controlled by fluctuations driven by electromagnetic kinetic instabilities. In this study the temperature anisotropy is driven self-consistently by the expansion. While the expansion favors an increase of parallel anisotropy ($T_\parallel>T_\perp$), the onset of the firehose instability will tend to decrease it. We show the results for a supersonic, subsonic, and static expansion flows, and suggest possible applications of the results for the solar wind and other stellar winds.

Energetic electrons are a common feature of interplanetary shocks and planetary bow shocks, and they are invoked as a key component of models of nonthermal radio emission, such as solar radio bursts. A simulation study is carried out of electron acceleration for high Mach number, quasi-perpendicular shocks, typical of the shocks in the solar wind. Two dimensional self-consistent hybrid shock simulations provide the electric and magnetic fields in which test particle electrons are followed. A range of different shock types, shock normal angles, and injection energies are studied. When the Mach number is low, or the simulation configuration suppresses fluctuations along the magnetic field direction, the results agree with theory assuming magnetic moment conserving reflection (or Fast Fermi acceleration), with electron energy gains of a factor only 2 - 3. For high Mach number, with a realistic simulation configuration, the shock front has a dynamic rippled character. The corresponding electron energization is radically different: Energy spectra display: (1) considerably higher maximum energies than Fast Fermi acceleration; (2) a plateau, or shallow sloped region, at intermediate energies 2 - 5 times the injection energy; (3) power law fall off with increasing energy, for both upstream and downstream particles, with a slope decreasing as the shock normal angle approaches perpendicular; (4) sustained flux levels over a broader region of shock normal angle than for adiabatic reflection. All these features are in good qualitative agreement with observations, and show that dynamic structure in the shock surface at ion scales produces effective scattering and can be responsible for making high Mach number shocks effective sites for electron acceleration.