We present a statistical analysis of near-relativistic (NR) solar energetic electron event spectra near 1au. We use measurements of the STEREO Solar Electron and Proton Telescope (SEPT) in the energyrange of 45-425 keV and utilize the SEPT electron event list containing all electron events observed bySTEREO A and STEREO B from 2007 through 2018. We select 781 events with significant signal tonoise ratios for our analysis and fit the spectra with single or broken power law functions of energy.We find 437 (344) events showing broken (single) power laws in the energy range of SEPT. The eventswith broken power laws show a mean break energy of about 120 keV. We analyze the dependence ofthe spectral index on the rise times and peak intensities of the events as well as on the presence ofrelativistic electrons. The results show a relation between the power law spectral index and the risetimes of the events with softer spectra belonging to rather impulsive events. Long rise-time events areassociated with hard spectra as well as with the presence of higher energy (>0.7 MeV) electrons. Thisgroup of events cannot be explained by a pure flare scenario but suggests an additional accelerationmechanism, involving a prolonged acceleration and/or injection of the particles. A dependence of thespectral index on the longitudinal separation from the parent solar source region was not found. Astatistical analysis of the spectral indices during impulsively rising events (rise times<20 minutes) isalso shown.

We calculate the interplanetary magnetic field path lengths traveled by electrons in solar electron events detected by the WIND 3DP instrument from $1994$ to $2016$. The velocity dispersion analysis method is applied for electrons at energies of $\sim$ $27$ keV to $310$ keV. Previous velocity dispersion analyses employ the onset times, which are often affected by instrumental effects and the pre-existing background flux, leading to large uncertainties. We propose a new method here. Instead of using the peak or onset time, we apply the velocity dispersion analysis to the times that correspond to the rising phase of the fluxes that are a fraction, $\eta$, of the peak flux. We perform statistical analysis on selected events whose calculated path lengths have uncertainties smaller than $0.1$ AU. The mean and standard deviation, ($\mu$, $\sigma$), of the calculated path lengths corresponding to $\eta=$ $3/4$, $1/2$, and $1/3$ of the peak flux is ($1.17$ AU, $0.17$ AU), ($1.11$ AU, $0.14$ AU), and ($1.06$ AU, $0.15$ AU). The distribution of the calculated path lengths is also well fitted by a Gaussian distribution for the $\eta=3/4$ and $1/3$ cases. These results suggest that in these electron events the interplanetary magnetic field topology is close to the nominal Parker spiral with little field line meandering. Our results have important implications for particles' perpendicular diffusion.

A knowledge of the particle escape time from the acceleration regions of many space and astrophysical sources is of critical importance in the analysis of emission signatures produced by these particles and in the determination of the acceleration and transport mechanisms at work. This paper addresses this general problem, in particular in solar flares, where in addition to scattering by turbulence, the magnetic field convergence from the acceleration region towards its boundaries also influences the particle escape. We test an (approximate) analytic relation between escape and scattering times, and the field convergence rate, based on the work of Malyshkin and Kulsrud (2001), valid for both strong and weak diffusion limits and isotropic pitch angle distribution of the injected particles, with a numerical model of particle transport. To this end, a kinetic Fokker-Planck transport model of particles is solved with a stochastic differential equation scheme assuming different initial pitch angle distributions. This approach enables further insights into the phase-space dynamics of the transport process, which would otherwise not be accessible. We find that in general the numerical results agree well with the analytic equation for the isotropic case, however, there are significant differences weak diffusion regime for non-isotopic cases, especially for distributions beamed along the magnetic field lines. The results are important in the interpretation of observations of energetic particles in solar flares, and other similar space and astrophysical acceleration sites and for the determination of acceleration-transport coefficients, commonly used in Fokker-Planck type kinetic equations.

Recent detections of high-energy $\gamma$-rays from behind-the-limb (BTL) solar flares by the \emph{Fermi $\gamma$-ray Space Telescope} pose a puzzle and challenge on the particle acceleration and transport mechanisms. In such events, the $\gamma$-ray emission region is located away from the BTL flare site by up to tens of degrees in heliogrpahic longitude. It is thus hypothesized that particles are accelerated at the shock driven by the coronal mass ejection (CME) and then travel from the shock downstream back to the front side of the Sun to produce the observed $\gamma$-rays. To test this scenario, we performed data-driven, global magnetohydrodynamics simulations of the CME associated with a well-observed BTL flare on 2014 September 1. We found that part of the CME-driven shock develops magnetic connectivity with the $\gamma$-ray emission region, facilitating transport of particles back to the Sun. Moreover, the observed increase in $\gamma$-ray flux is temporally correlated with (1) the increase of the shock compression ratio and (2) the presence of a quasi-perpendicular shock over the area that is magnetically connected to the $\gamma$-ray emitting region, both conditions favoring the diffusive shock acceleration (DSA) of particles. These results support the above hypothesis and can help resolve another puzzle, i.e., long-duration (up to 20 hours) $\gamma$-rays flares. We suggest that, in addition to DSA, stochastic acceleration by plasma turbulence may also play a role, especially in the shock downstream region and during the early stage when the shock Alfv\'{e}n Mach number is small.

Particles are accelerated to very high, non-thermal energies in solar and space plasma environments. While energy spectra of accelerated electrons often exhibit a power law, it remains unclear how electrons are accelerated to high energies and what processes determine the power-law index $\delta$. Here, we review previous observations of the power-law index $\delta$ in a variety of different plasma environments with a particular focus on sub-relativistic electrons. It appears that in regions more closely related to magnetic reconnection (such as the `above-the-looptop' solar hard X-ray source and the plasma sheet in Earth's magnetotail), the spectra are typically soft ($\delta \gtrsim$ 4). This is in contrast to the typically hard spectra ($\delta \lesssim$ 4) that are observed in coincidence with shocks. The difference implies that shocks are more efficient in producing a larger non-thermal fraction of electron energies when compared to magnetic reconnection. A caveat is that during active times in Earth's magnetotail, $\delta$ values seem spatially uniform in the plasma sheet, while power-law distributions still exist even in quiet times. The role of magnetotail reconnection in the electron power-law formation could therefore be confounded with these background conditions. Because different regions have been studied with different instrumentations and methodologies, we point out a need for more systematic and coordinated studies of power-law distributions for a better understanding of possible scaling laws in particle acceleration as well as their universality.

Insights into the processes of Solar Energetic Particle (SEP) propagation are essential for understanding how solar eruptions affect the radiation environment of near-Earth space. SEP propagation is influenced by turbulent magnetic fields in the solar wind, resulting in stochastic transport of SEPs to Earth. Multi-spacecraft observations suggest that the cross-field propagation shapes the SEP fluxes at Earth strongly. However, modelling SEP cross-field transport as spatial diffusion has been shown to be insufficient without use of unrealistically large cross-field diffusion coefficients. Recent work has shown that the early-time propagation of energetic particles across the mean field direction in turbulent fields is not diffusive, as the particles propagating along meandering field lines. This early-time transport mode results in fast access of the particles across the mean field direction, in agreement with the SEP observations. In this work, we demonstrate the significance of turbulence strength on evolution of the SEP radiation environment near Earth. We calculate the transport parameters with a turbulence transport model, parametrised by the SEP parallel scattering mean free path at 1~AU, $\lambda_{\parallel}^{*}$, and show that the parallel and cross-field transport are connected, with conditions resulting in slow parallel transport corresponding to wider events. We find a scaling $\sigma_{\phi,\, \mathrm{max}}\propto (1/\lambda_{\parallel}^{*})^{1/4}$ for the Gaussian fitting of the longitudinal distribution of maximum intensities. The longitudes with highest intensities are shifted towards the west for strong scattering conditions. Our results emphasise the importance of understanding both the SEP transport and the interplanetary turbulence conditions for modelling and predicting the SEP radiation environment at Earth.

In this review, an overview of the recent history of stochastic differential equations (SDEs) in application to particle transport problems in space physics and astrophysics is given. The aim is to present a helpful working guide to the literature and at the same time introduce key principles of the SDE approach via "toy models". Using these examples, we hope to provide an easy way for newcomers to the field to use such methods in their own research. Aspects covered are the solar modulation of cosmic rays, diffusive shock acceleration, galactic cosmic ray propagation and solar energetic particle transport. We believe that the SDE method, due to its simplicity and computational efficiency on modern computer architectures, will be of significant relevance in energetic particle studies in the years to come.

Flares close to the solar limb, where the footpoints are occulted, can reveal the spectrum and structure of the coronal loop-top source in X-rays. We aim at studying the properties of the corresponding energetic electrons near their acceleration site, without footpoint contamination. To this end, a statistical study of partially occulted flares observed with RHESSI is presented here, covering a large part of solar cycles 23 and 24. We perform a detailed spectra, imaging and light curve analysis for 116 flares and include contextual observations from SDO and STEREO when available, providing further insights into flare emission that was previously not accessible. We find that most spectra are fitted well with a thermal component plus a broken power-law, non-thermal component. A thin-target kappa distribution model gives satisfactory fits after the addition of a thermal component. X-rays imaging reveals small spatial separation between the thermal and non-thermal components, except for a few flares with a richer coronal source structure. A comprehensive light curve analysis shows a very good correlation between the derivative of the soft X-ray flux (from GOES) and the hard X-rays for a substantial number of flares, indicative of the Neupert effect. The results confirm that non-thermal particles are accelerated in the corona and estimated timescales support the validity of a thin-target scenario with similar magnitudes of thermal and non-thermal energy fluxes.

At hard X-ray energies, the bright footpoint emission from solar flare loops often prevents a detailed analysis of the weaker loop-top source morphology due to the limited dynamic range available for X-ray imaging. Here, we study the X1.3 April 25, 2014 flare with the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI). This partially occulted limb flare allows the analysis of the loop-top emission in isolation. We present results on the flare light curve at different energies, the source morphology from X-ray imaging and a detailed spectral analysis of the different source components by imaging spectroscopy. The loop-top source, a likely site of particle acceleration, shows a clear composition of different emission components. The results indicate the opportunities that detailed imaging of hard X-rays can provide to learn about particle acceleration, transport and heating processes in solar flares.

Context. Current solar energetic particle (SEP) propagation models describe the effects of interplanetary plasma turbulence on SEPs as diffusion, using a Fokker-Planck (FP) equation. However, FP models cannot explain the observed fast access of SEPs across the average magnetic field to regions that are widely separated in longitude within the heliosphere without using unrealistically strong cross-field diffusion. Aims. We study whether the recently suggested early non-diffusive phase of SEP propagation can explain the wide SEP events with realistic particle transport parameters. Methods. We used a novel model that accounts for the SEP propagation along field lines that meander as a result of plasma turbulence. Such a non-diffusive propagation mode has been shown to dominate the SEP cross-field propagation early in the SEP event history. We compare the new model to the traditional approach, and to SEP observations. Results. Using the new model, we reproduce the observed longitudinal extent of SEP peak fluxes that are characterised by a Gaussian profile with $\sigma=30-50^\circ$, while current diffusion theory can only explain extents of 11$^\circ$ with realistic diffusion coefficients. Our model also reproduces the timing of SEP arrival at distant longitudes, which cannot be explained using the diffusion model. Conclusions. The early onset of SEPs over a wide range of longitudes can be understood as a result of the effects of magnetic field-line random walk in the interplanetary medium and requires an SEP transport model that properly describes the non-diffusive early phase of SEP cross-field propagation.

Starting from an exact, steady-state, force-free solution of the magnetohydrodynamic (MHD) equations, we investigate how resistive current layers are induced by perturbing line-tied three-dimensional magnetic equilibria. This is achieved by the superposition of a weak perturbation field in the domain, in contrast to studies where the boundary is driven by slow motions, like those present in photospheric active regions. Our aim is to quantify how the current structures are altered by the contribution of so called quasi-separatrix layers (QSLs) as the null point is shifted outside the computational domain. Previous studies based on magneto-frictional relaxation have indicated that, despite the severe field line gradients of the QSL, the presence of a null is vital in maintaining fast reconnection. Here, we explore this notion using highly resolved simulations of the full MHD evolution. We show that for the null-point configuration, the resistive scaling of the peak current density is close to $J\sim\eta^{-1}$, while the scaling is much weaker, i.e. $J\sim\eta^{-0.4}$, when only the QSL connectivity gradients provide a site for the current accumulation.

In-situ spacecraft observations recently suggested that the transport of energetic particles accelerated at heliospheric shocks can be anomalous, i.e. the mean square displacement can grow non-linearly in time. In particular, a new analysis technique has permitted the study of particle transport properties from energetic particle time profiles upstream of interplanetary shocks. Indeed, the time/spatial power laws of the differential intensity upstream of several shocks are indicative of superdiffusion. A complete determination of the key parameters of superdiffusive transport comprises the power-law index, the superdiffusion coefficient, the related transition scale at which the energetic particle profiles turn to decay as power laws, and the energy spectral index of the shock accelerated particles. Assuming large-scale spatial homogeneity of the background plasma, the power-law behaviour can been derived from both a (microscopic) propagator formalism and a (macroscopic) fractional transport equation. We compare the two approaches and find a relation between the diffusion coefficients used in the two formalisms. Based on the assumption of superdiffusive transport, we quantitatively derive these parameters by studying energetic particle profiles observed by the Ulysses and Voyager 2 spacecraft upstream of shocks in the heliosphere, for which a superdiffusive particle transport has previously been observed. Further, we have jointly studied the electron energy spectra, comparing the values of the spectral indices observed with those predicted by the standard diffusive shock acceleration theory and by a model based on superdiffusive transport. For a number of interplanetary shocks and for the solar wind termination shock, for the first time we obtain the anomalous diffusion constants and the scale at which the probability of particle free paths changes to a power-law...

Diffusive cosmic-ray transport in nonuniform large-scale magnetic fields in the presence of boundaries is considered. Reflecting and absorbing boundary conditions are derived for a modified telegraph equation with a convective term. Analytical and numerical solutions of illustrative boundary problems are presented. The applicability and accuracy of the telegraph approximation for focused cosmic-ray transport in the presence of boundaries are discussed, and potential applications to modeling cosmic-ray transport are noted.

The open problem of how singular current structures form in line-tied, three-dimensional magnetic fields is addressed. A Lagrangian magneto-frictional relaxation method is employed to model the field evolution towards the final near-singular state. Our starting point is an exact force-free solution of the governing magnetohydrodynamic (MHD) equations which is sufficiently general to allow for topological features like magnetic nulls to be inside or outside the computational domain, depending on a simple set of parameters. Quasi-separatrix layers (QSLs) are present in these structures and together with the magnetic nulls, they significantly influence the accumulation of current. It is shown that perturbations affecting the lateral boundaries of the configuration lead not only to collapse around the magnetic null, but also to significant QSL currents. Our results show that once a magnetic null is present, the developing currents are always attracted to that specific location and show a much stronger scaling with resolution than the currents which form along the QSL. In particular, the null point scalings can be consistent with models of "fast" reconnection. The QSL currents also appear to be unbounded but give rise to weaker singularities, independent of the perturbation amplitude.

The diffusion approximation to the Fokker-Planck equation is commonly used to model the transport of solar energetic particles in interplanetary space. In this study, we present exact analytical predictions of a higher order telegraph approximation for particle transport and compare them with the corresponding predictions of the diffusion approximation and numerical solutions of the full Fokker-Planck equation. We specifically investigate the role of the adiabatic focusing effect of a spatially varying magnetic field on an evolving particle distribution. Comparison of the analytical and numerical results shows that the telegraph approximation reproduces the particle intensity profiles much more accurately than does the diffusion approximation, especially when the focusing is strong. However, the telegraph approximation appears to offer no significant advantage over the diffusion approximation for calculating the particle anisotropy. The telegraph approximation can be a useful tool for describing both diffusive and wavelike aspects of the cosmic-ray transport.

Motivated by recent applications of superdiffusive transport models to shock-accelerated particle distributions in the heliosphere, we solve analytically a one-dimensional fractional diffusion-advection equation for the particle density. We derive an exact Fourier transform solution, simplify it in a weak diffusion approximation, and compare the new solution with previously available analytical results and with a semi-numerical solution based on a Fourier series expansion. We apply the results to the problem of describing the transport of energetic particles, accelerated at a traveling heliospheric shock. Our analysis shows that significant errors may result from assuming an infinite initial distance between the shock and the observer. We argue that the shock travel time should be a parameter of a realistic superdiffusive transport model.

Context. Spacecraft observations have motivated the need for a refined description of the phase-space distribution function. Of particular importance is the pitch-angle diffusion coefficient that occurs in the Fokker-Planck transport equation. Aims. Simulations and analytical test-particle theories are compared to verify the diffusion description of particle transport, which does not allow for non-Markovian behavior. Methods. A Monte-Carlo simulation code was used to trace the trajectories of test particles moving in turbulent magnetic fields. From the ensemble average, the pitch-angle Fokker-Planck coefficient is obtained via the mean square displacement. Results. It is shown that, while excellent agreement with analytical theories can be obtained for slab turbulence, considerable deviations are found for isotropic turbulence. In addition, all Fokker-Planck coefficients tend to zero for high time values.

Galactic transport models for cosmic rays involve the diffusive motion of these particles in the interstellar medium. Due to the large-scale structured galactic magnetic field this diffusion is anisotropic with respect to the local field direction. We included this transport effect along with continuous loss processes in a quantitative model of galactic propagation for cosmic ray protons which is based on stochastic differential equations. We calculated energy spectra at different positions along the Sun's galactic orbit and compared them to the isotropic diffusion case. The results show that a larger amplitude of variation as well as different spectral shapes are obtained in the introduced anisotropic diffusion scenario and emphasize the need for accurate galactic magnetic field models.

The spatial diffusion of cosmic rays in turbulent magnetic fields can, in the most general case, be fully anisotropic, i.e. one has to distinguish three diffusion axes in a local, field-aligned frame. We reexamine the transformation for the diffusion tensor from this local to a global frame, in which the Parker transport equation for energetic particles is usually formulated and solved. Particularly, we generalize the transformation formulas to allow for an explicit choice of two principal local perpendicular diffusion axes. This generalization includes the 'traditional' diffusion tensor in the special case of isotropic perpendicular diffusion. For the local frame, we motivate the choice of the Frenet-Serret trihedron which is related to the intrinsic magnetic field geometry. We directly compare the old and the new tensor elements for two heliospheric magnetic field configurations, namely the hybrid Fisk and the Parker field. Subsequently, we examine the significance of the different formulations for the diffusion tensor in a standard 3D model for the modulation of galactic protons. For this we utilize a numerical code to evaluate a system of stochastic differential equations equivalent to the Parker transport equation and present the resulting modulated spectra. The computed differential fluxes based on the new tensor formulation deviate from those obtained with the 'traditional' one (only valid for isotropic perpendicular diffusion) by up to 60% for energies below a few hundred MeV depending on heliocentric distance.