We study the solar eruptive event on 2017 September 10 that produced long-lasting $>$100 MeV $\gamma$-ray emission and a ground level enhancement (GLE72). The origin of the high-energy ions producing late-phase gamma-ray emission (LPGRE) is still an open question, but a possible explanation is proton acceleration at coronal shocks produced by coronal mass ejections. We examine a common shock acceleration origin for both the LPGRE and GLE72. The $\gamma$-ray emission observed by the Fermi-Large Area Telescope exhibits a weak impulsive phase, consistent with that observed in hard X-and $\gamma$-ray line flare emissions, and what appear to be two distinct stages of LPGRE. From a detailed modeling of the shock wave, we derive the 3D distribution and temporal evolution of the shock parameters, and we examine the shock wave magnetic connection with the visible solar disk. The evolution of shock parameters on field lines returning to the visible disk, mirrors the two stages of LPGRE. We find good agreement between the time history of $>$100 MeV $\gamma$-rays and one produced by a basic shock acceleration model. The time history of shock parameters magnetically mapped to Earth agrees with the rates observed by the Fort Smith neutron monitor during the first hour of the GLE72 if we include a 30% contribution of flare-accelerated protons during the first 10 minutes, having a release time following the time history of nuclear $\gamma$-rays. Our analysis provides compelling evidence for a common shock origin for protons producing the LPGRE and most of the particles observed in GLE72.

This review aims at providing an up-to-date status and a general introduction to the subject of the numerical study of energetic particle acceleration and transport in turbulent astrophysical flows. The subject is also complemented by a short overview of recent progresses obtained in the domain of laser plasma experiments. We review the main physical processes at the heart of the production of a non-thermal distribution in both Newtonian and relativistic astrophysical flows, namely the first and second order Fermi acceleration processes. We also discuss shock drift and surfing acceleration, two processes important in the context of particle injection in shock acceleration. We analyze with some details the particle-in-cell (PIC) approach used to describe particle kinetics. We review the main results obtained with PIC simulations in the recent years concerning particle acceleration at shocks and in reconnection events. The review discusses the solution of Fokker-Planck problems with application to the study of particle acceleration at shocks but also in hot coronal plasmas surrounding compact objects. We continue by considering large scale physics. We describe recent developments in magnetohydrodynamic (MHD) simulations. We give a special emphasize on the way energetic particle dynamics can be coupled to MHD solutions either using a multi-fluid calculation or directly coupling kinetic and fluid calculations. This aspect is mandatory to investigate the acceleration of particles in the deep relativistic regimes to explain the highest Cosmic Ray energies.

Weakly magnetized, relativistic collisionless shock waves are not only the natural offsprings of relativistic jets in high-energy astrophysical sources, they are also associated with some of the most outstanding displays of energy dissipation through particle acceleration and radiation. Perhaps their most peculiar and exciting feature is that the magnetized turbulence that sustains the acceleration process, and (possibly) the secondary radiation itself, is self-excited by the accelerated particles themselves, so that the phenomenology of these shock waves hinges strongly on the microphysics of the shock. In this review, we draw a status report of this microphysics, benchmarking analytical arguments with particle-in-cell simulations, and extract consequences of direct interest to the phenomenology, regarding in particular the so-called microphysical parameters used in phenomenological studies.

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.

The gyro-resonant cosmic-ray (CR) streaming instability is believed to play a crucial role in CR transport, leading to growth of Alfv\'en waves at small scales that scatter CRs, and impacts the interaction of CRs with the ISM on large scales. However, extreme scale separation ($\lambda \ll \rm pc$), low cosmic ray number density ($n_{\rm CR}/n_{\rm ISM} \sim 10^{-9}$), and weak CR anisotropy ($\sim v_A/c$) pose strong challenges for proper numerical studies of this instability on a microphysical level. Employing the recently developed magnetohydrodynamic-particle-in-cell (MHD-PIC) method, which has unique advantages to alleviate these issues, we conduct one-dimensional simulations that quantitatively demonstrate the growth and saturation of the instability in the parameter regime consistent with realistic CR streaming in the large-scale ISM. Our implementation of the $\delta f$ method dramatically reduces Poisson noise and enables us to accurately capture wave growth over a broad spectrum, equally shared between left and right handed Alfv\'en modes. We are also able to accurately follow the quasi-linear diffusion of CRs subsequent to wave growth, which is achieved by employing phase randomization across periodic boundaries. Full isotropization of the CRs in the wave frame requires pitch angles of most CRs to efficiently cross $90^\circ$, and can be captured in simulations with relatively high wave amplitude and/or high spatial resolution. We attribute this crossing to non-linear wave-particle interaction (rather than mirror reflection) by investigating individual CR trajectories. We anticipate our methodology will open up opportunities for future investigations that incorporate additional physics.

Perpendicular relativistic ($\gamma_0=10$) shocks in magnetized pair plasmas are investigated using two dimensional Particle-in-Cell simulations. A systematic survey, from unmagnetized to strongly magnetized shocks, is presented accurately capturing the transition from Weibel-mediated to magnetic-reflection-shaped shocks. This transition is found to occur for upstream flow magnetizations $10^{-3}<\sigma<10^{-2}$ at which a strong perpendicular net current is observed in the precursor, driving the so-called current-filamentation instability. The global structure of the shock and shock formation time are discussed. The MHD shock jump conditions are found in good agreement with the numerical results, except for $10^{-4} < \sigma < 10^{-2}$ where a deviation up to 10\% is observed. The particle precursor length converges toward the Larmor radius of particles injected in the upstream magnetic field at intermediate magnetizations. For $\sigma>10^{-2}$, it leaves place to a purely electromagnetic precursor following from the strong emission of electromagnetic waves at the shock front. Particle acceleration is found to be efficient in weakly magnetized perpendicular shocks in agreement with previous works, and is fully suppressed for $\sigma > 10^{-2}$. Diffusive Shock Acceleration is observed only in weakly magnetized shocks, while a dominant contribution of Shock Drift Acceleration is evidenced at intermediate magnetizations. The spatial diffusion coefficients are extracted from the simulations allowing for a deeper insight into the self-consistent particle kinematics and scale with the square of the particle energy in weakly magnetized shocks. These results have implications for particle acceleration in the internal shocks of AGN jets and in the termination shocks of Pulsar Wind Nebulae.

Context. The observation of >100 MeV {\gamma}-rays in the minutes to hours following solar flares suggests that high-energy particles interacting in the solar atmosphere can be stored and/or accelerated for long time periods. The occasions when {\gamma}-rays are detected even when the solar eruptions occurred beyond the solar limb as viewed from Earth provide favorable viewing conditions for studying the role of coronal shocks driven by coronal mass ejections (CMEs) in the acceleration of these particles. Aims: In this paper, we investigate the spatial and temporal evolution of the coronal shocks inferred from stereoscopic observations of behind-the-limb flares to determine if they could be the source of the particles producing the {\gamma}-rays. Methods: We analyzed the CMEs and early formation of coronal shocks associated with {\gamma}-ray events measured by the Fermi-Large Area Telescope (LAT) from three eruptions behind the solar limb as viewed from Earth on 2013 Oct. 11, 2014 Jan. 06 and Sep. 01. We used a 3D triangulation technique, based on remote-sensing observations to model the expansion of the CME shocks from above the solar surface to the upper corona. Coupling the expansion model to various models of the coronal magnetic field allowed us to derive the time-dependent distribution of shock Mach numbers and the magnetic connection of particles produced by the shock to the solar surface visible from Earth. Results: The reconstructed shock fronts for the three events became magnetically connected to the visible solar surface after the start of the flare and just before the onset of the >100 MeV {\gamma}-ray emission. The shock surface at these connections also exhibited supercritical Mach numbers required for significant particle energization. [...] (Abridged)

The remoteness of the Sun and the harsh conditions prevailing in the solar corona have so far limited the observational data used in the study of solar physics to remote-sensing observations taken either from the ground or from space. In contrast, the `solar wind laboratory' is directly measured in situ by a fleet of spacecraft measuring the properties of the plasma and magnetic fields at specific points in space. Since 2007, the solar-terrestrial relations observatory (STEREO) has been providing images of the solar wind that flows between the solar corona and spacecraft making in-situ measurements. This has allowed scientists to directly connect processes imaged near the Sun with the subsequent effects measured in the solar wind. This new capability prompted the development of a series of tools and techniques to track heliospheric structures through space. This article presents one of these tools, a web-based interface called the 'Propagation Tool' that offers an integrated research environment to study the evolution of coronal and solar wind structures, such as Coronal Mass Ejections (CMEs), Corotating Interaction Regions (CIRs) and Solar Energetic Particles (SEPs). These structures can be propagated from the Sun outwards to or alternatively inwards from planets and spacecraft situated in the inner and outer heliosphere. In this paper, we present the global architecture of the tool, discuss some of the assumptions made to simulate the evolution of the structures and show how the tool connects to different databases.

We study the link between an expanding coronal shock and the energetic particles measured near Earth during the Ground Level Enhancement (GLE) of 17 May 2012. We developed a new technique based on multipoint imaging to triangulate the 3-D expansion of the shock forming in the corona. It uses images from three vantage points by mapping the outermost extent of the coronal region perturbed by the pressure front. We derive for the first time the 3-D velocity vector and the distribution of Mach numbers, $M_{FM}$, of the entire front as a function of time. Our approach uses magnetic field reconstructions of the coronal field, full magneto-hydrodynamic simulations and imaging inversion techniques. We find that the highest $M_{FM}$ values appear near the coronal neutral line within a few minutes of the Coronal Mass Ejection (CME) onset; this neutral line is usually associated with the source of the heliospheric current and plasma sheet. We illustrate the variability of the shock speed, shock geometry and Mach number along different modeled magnetic field lines. Despite the level of uncertainty in deriving the shock Mach numbers, all employed reconstruction techniques show that the release time of GeV particles occurs when the coronal shock becomes super-critical ($M_{FM}>3$). Combining in-situ measurements with heliospheric imagery, we also demonstrate that magnetic connectivity between the accelerator (the coronal shock of 17 May 2012) and the near-Earth environment is established via a magnetic cloud that erupted from the same active region roughly five days earlier.

The origin of the slow solar wind is still a topic of much debate. The continual emergence of small transient structures from helmet streamers is thought to constitute one of the main sources of the slow wind. Determining the height at which these transients are released is an important factor in determining the conditions under which the slow solar wind forms. To this end, we have carried out a multipoint analysis of small transient structures released from a north-south tilted helmet streamer into the slow solar wind over a broad range of position angles during Carrington Rotation 2137. Combining the remote-sensing observations taken by the Solar-TErrestrial RElations Observatory (STEREO) mission with coronagraphic observations from the Solar and Heliospheric Observatory (SoHO) spacecraft, we show that the release of such small transient structures (often called blobs), which subsequently move away from the Sun, is associated with the concomitant formation of transient structures collapsing back towards the Sun; the latter have been referred to by previous authors as raining inflows. This is the first direct association between outflowing blobs and raining inflows, which locates the formation of blobs above the helmet streamers and gives strong support that the blobs are released by magnetic reconnection.

The systematic monitoring of the solar wind in high-cadence and high-resolution heliospheric images taken by the Solar-Terrestrial Relation Observatory (STEREO) spacecraft permits the study of the spatial and temporal evolution of variable solar wind flows from the Sun out to 1~AU, and beyond. As part of the EU Framework 7 (FP7) Heliospheric Cataloguing, Analysis and Techniques Service (HELCATS) project, we have generated a catalogue listing the properties of 190 corotating structures well-observed in images taken by the Heliospheric Imager instruments on-board STEREO-A. We present here one of very few long-term analyses of solar wind structures advected by the background solar wind. This analysis confirms that most of the corotating density structures detected by the heliospheric imagers comprises a series of density inhomogeneities advected by the slow solar wind that eventually become entrained by stream interaction regions. We have derived the spatial-temporal evolution of each of these corotating density structures by using a well-established fitting technique. The mean radial propagation speed of the corotating structures is found to be $311 \pm 31$ km~s$^{-1}$. We predicted the arrival time of each corotating density structure at different probes. We show that the speeds of the corotating density structures derived using our fitting technique track well the long-term variation of the radial speed of the slow solar wind during solar minimum years (2007--2008). Furthermore, we demonstrate that these features originate near the coronal neutral line that eventually becomes the heliospheric current sheet.

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 report here on a fast current-driven instability at relativistic collisionless shocks, triggered by the perpendicular current carried by the supra-thermal particles as they gyrate around the background magnetic field in the shock precursor. We show that this instability grows faster than any other instability studied so far in this context, and we argue that it is likely to shape the physics of the shock and of particle acceleration in a broad parameter range.

The physics of collisionless relativistic shocks with a moderate magnetization is presented. Micro-physics is relevant to explain the most energetic radiative phenomena of Nature, namely that of the termination shock of Gamma Ray Bursts. A transition towards Fermi process occurs for decreasing magnetization around a critical value which turns out to be the condition for the scattering to break the mean field inhibition. Scattering is produced by magnetic micro-turbulence driven by the current carried by returning particles, which had not been considered till now, but turns out to be more intense than Weibel's one around the transition. The current is also responsible for a buffer effect on the motion of the incoming flow, on which the threshold for the onset of turbulence depends.

The physics of instabilities in the precursor of relativistic collisionless shocks is of broad importance in high energy astrophysics, because these instabilities build up the shock, control the particle acceleration process and generate the magnetic fields in which the accelerated particles radiate. Two crucial parameters control the micro-physics of these shocks: the magnetization of the ambient medium and the Lorentz factor of the shock front; as of today, much of this parameter space remains to be explored. In the present paper, we report on a new instability upstream of electron-positron relativistic shocks and we argue that this instability shapes the micro-physics at moderate magnetization levels and/or large Lorentz factors. This instability is seeded by the electric current carried by the accelerated particles in the shock precursor as they gyrate around the background magnetic field. The compensation current induced in the background plasma leads to an unstable configuration, with the appearance of charge neutral filaments carrying a current of the same polarity, oriented along the perpendicular current. This ``current-driven filamentation'' instability grows faster than any other instability studied so far upstream of relativistic shocks, with a growth rate comparable to the plasma frequency. Furthermore, the compensation of the current is associated with a slow-down of the ambient plasma as it penetrates the shock precursor (as viewed in the shock rest frame). This slow-down of the plasma implies that the ``current driven filamentation'' instability can grow for any value of the shock Lorentz factor, provided the magnetization \sigma <~ 10^{-2}. We argue that this instability explains the results of recent particle-in-cell simulations in the mildly magnetized regime.

Collisionless relativistic shocks have been the focus of intense theoretical and numerical investigations in recent years. The acceleration of particles, the generation of electromagnetic microturbulence and the building up of a shock front are three interrelated essential ingredients of a relativistic collisionless shock wave. In this paper we investigate two issues of importance in this context: (1) the transport of suprathermal particles in the excited microturbulence upstream of the shock and its consequences regarding particle acceleration; (2) the preheating of incoming background electrons as they cross the shock precursor and experience relativistic oscillations in the microturbulent electric fields. We place emphasis on the importance of the motion of the electromagnetic disturbances relatively to the background plasma and to the shock front. This investigation is carried out for the two major instabilities involved in the precursor of relativistic shocks, the filamentation instability and the oblique two stream instability. Finally, we use our results to discuss the maximal acceleration at the external shock of a gamma-ray burst; we find in particular a maximal synchrotron photon energy of the order of a few GeV.

Various astrophysical studies have motivated the investigation of the transport of high energy particles in magnetic turbulence, either in the source or en route to the observation sites. For strong turbulence and large rigidity, the pitch-angle scattering rate is governed by a simple law involving a mean free path that increases proportionally to the square of the particle energy. In this paper, we show that perpendicular diffusion deviates from this behavior in the presence of a mean field. We propose an exact theoretical derivation of the diffusion coefficients and show that a mean field significantly changes the transverse diffusion even in the presence of a stronger turbulent field. In particular, the transverse diffusion coefficient is shown to reach a finite value at large rigidity instead of increasing proportionally to the square of the particle energy. Our theoretical derivation is corroborated by a dedicated Monte Carlo simulation. We briefly discuss several possible applications in astrophysics.