Solar Energetic Particles (SEP) are an integral part of the physical processes related with Space Weather. We present a review for the acceleration mechanisms related to the explosive phenomena (flares and/or CMEs) inside the solar corona. For more than 40 years, the main 2D cartoon representing our understanding of the explosive phenomena inside the solar corona remained almost unchanged. The acceleration mechanisms related to solar flares and CMEs also remained unchanged and were part of the same cartoon. In this review, we revise the standard cartoon and present evidence from recent global MHD simulations that supports the argument that explosive phenomena will lead to the spontaneous formation of current sheets in different parts of the erupting magnetic structure. The evolution of the large scale current sheets and their fragmentation will lead to strong turbulence and turbulent reconnection during solar flares and turbulent shocks. In other words, the acceleration mechanism in flares and CME-driven shocks may be the same, and their difference will be the overall magnetic topology, the ambient plasma parameters, and the duration of the unstable driver.

We consider a large scale environment of turbulent reconnection that is fragmented into a number of randomly distributed Unstable Current Sheets (UCS), and we statistically analyze the acceleration of particles within this environment. We address two important cases of acceleration mechanisms when the particles interact with the UCS: (a) electric field acceleration, and (b) acceleration through reflection at contracting islands. Electrons and ions are accelerated very efficiently, attaining an energy distribution of power-law shape with an index $1-2$, depending on the acceleration mechanism. The transport coefficients in energy space are estimated from the test-particle simulation data, and we show that the classical Fokker-Planck (FP) equation fails to reproduce the simulation results when the transport coefficients are inserted into it and it is solved numerically. The cause for this failure is that the particles perform Levy flights in energy space, the distributions of energy increments exhibit power-law tails. We then use the fractional transport equation (FTE) derived by Isliker et al., 2017, whose parameters and the order of the fractional derivatives are inferred from the simulation data, and, solving the FTE numerically, we show that the FTE successfully reproduces the kinetic energy distribution of the test-particles. We discuss in detail the analysis of the simulation data and the criteria that allow judging the appropriateness of either an FTE or a classical FP equation as a transport model.

The aim of this study is to analyze the interaction of charged particles (ions and electrons) with randomly formed particle scatterers (e.g.\ large scale local "magnetic fluctuations" or "coherent magnetic irregularities"), using the set up proposed initially by \cite{Fermi49}. These scatterers are formed by the explosive magnetic energy release and propagate with the Alfv\'en speed along the irregular magnetic fields. They are large scale local fluctuations ($\delta B/B \approx 1$), randomly distributed inside the unstable magnetic topology and will here be called {\bf Alfv\'enic Scatterers (AS)}. We constructed a 3D grid on which a small fraction of randomly chosen grid points are acting as AS. In particular, we study how a large number of test particles evolve inside a collection of AS, analyzing the evolution of their energy distribution and their escape time distribution. We use a well established method to estimate the transport coefficients directly from the trajectories of the particles. Using the estimated transport coefficients and solving the Fokker-Planck (FP) equation numerically, we can recover the energy distribution of the particles. We have shown that the Stochastic Fermi Energization (SFE) of mildly relativistic and relativistic plasma can heat and accelerate the tail of the ambient particle distribution as predicted by \citet{Parker58} and \citet{Ramaty79}. The temperature of the hot plasma and the tail of the energetic particles depend on the mean free path ($\lambda_{\rm sc}$) of the particles between the scatterers inside the energization volume.

Turbulent flows in the solar wind, large scale current sheets, multiple current sheets, and shock waves lead to the formation of environments in which a dense network of current sheets is established and sustains "turbulent reconnection". We constructed a 2D grid on which a number of randomly chosen grid points are acting as scatterers (i.e. magnetic clouds or current sheets). Our goal is to examine how test particles respond inside this large scale collection of scatterers. We study the energy gain of individual particles, the evolution of their energy distribution and their escape time distribution. We have developed a new method to estimate the transport coefficients from the dynamics of the interaction of the particles with the scatterers. Replacing the "magnetic clouds" with current sheets, we have proven that the energization processes can be more efficient depending on the strength of the effective electric fields inside the current sheets and their statistical properties. Using the estimated transport coefficients and solving the Fokker-Planck (FP) equation we can recover the energy distribution of the particles only for the sstochastic Fermi process. We have shown that the evolution of the particles inside a turbulent reconnecting volume is not a solution of the FP equation, since the interaction of the particles with the current sheets is "anomalous", in contrast to the case of the second order Fermi process.

We model nanoflare heating of extrapolated active-region coronal loops via the acceleration of electrons and protons in Harris-type current sheets. The kinetic energy of the accelerated particles is estimated using semi-analytical and test-particle-tracing approaches. Vector magnetograms and photospheric Doppler velocity maps of NOAA active region 09114, recorded by the Imaging Vector Magnetograph (IVM), were used for this analysis. A current-free field extrapolation of the active-region corona was first constructed. The corresponding Poynting fluxes at the footpoints of 5000 extrapolated coronal loops were then calculated. Assuming that reconnecting current sheets develop along these loops, we utilized previous results to estimate the kinetic-energy gain of the accelerated particles and we related this energy to nanoflare heating and macroscopic loop characteristics. Kinetic energies of 0.1 to 8 keV (for electrons) and 0.3 to 470 keV (for protons) were found to cause heating rates ranging from $10^{-6}$ to 1 $\mathrm{erg\, s^{-1} cm^{-3}}$. Hydrodynamic simulations show that such heating rates can sustain plasma in coronal conditions inside the loops and generate plasma thermal distributions which are consistent with active region observations. We concluded the analysis by computing the form of X-ray spectra generated by the accelerated electrons using the thick target approach that were found to be in agreement with observed X-ray spectra, thus supporting the plausibility of our nanoflare-heating scenario.

This work examines the relation between the fractal properties of the photospheric magnetic patterns and those of the coronal magnetic fields in solar active regions.Our results show no correlation between the fractal dimensions of the observed 2d photospheric structures and the extrapolated unstable volumes in the corona, when nonlinear force-free extrapolation is used. This result is independent of efforts to (1) bring the photospheric magnetic fields closer to a nonlinear force-free equilibrium and (2) omit the lower part of the modeled magnetic field volume that is almost completely filled by unstable volumes. A significant correlation between the fractal dimensions of the photospheric and coronal magnetic features is only observed at the zero level (lower limit) of approximation of a current-free (potential) magnetic field extrapolation. We conclude that the complicated transition from photospheric non-force-free fields to coronal force-free ones hampers any direct correlation between the fractal dimensions of the 2d photospheric patterns and their 3d counterparts in the corona at the nonlinear force-free limit.

The acceleration of charged particles (electrons and protons) in flaring solar active regions is analyzed by numerical experiments. The acceleration is modelled as a stochastic process taking place by the interaction of the particles with local magnetic reconnection sites via multiple steps. Two types of local reconnecting topologies are studied: the Harris-type and the X-point. A formula for the maximum kinetic energy gain in a Harris-type current sheet, found in a previous work of ours, fits well the numerical data for a single step of the process. A generalization is then given approximating the kinetic energy gain through an X-point. In the case of the multiple step process, in both topologies the particles' kinetic energy distribution is found to acquire a practically invariant form after a small number of steps. This tendency is interpreted theoretically. Other characteristics of the acceleration process are given, such as the mean acceleration time and the pitch angle distributions of the particles.

In Isliker et al. (2000b), an extended cellular automaton (X-CA) model for solar flares was introduced. In this model, the interpretation of the model's grid-variable is specified, and the magnetic field, the current, and an approximation to the electric field are yielded, all in a way that is consistent with Maxwell's and the MHD equations. Here, we reveal which relevant plasma physical processes are implemented by the X-CA model and in what form, and what global physical set-up is assumed by this model when it is in its natural state (SOC). The basic results are: (1) On large-scales, all variables show characteristic quasi-symmetries. (2) The global magnetic topology forms either (i) closed magnetic field lines, or (ii) an arcade of field lines above the bottom plane line, if the model is slightly modified. (3) In case of the magnetic topology (ii), loading can be interpreted as if there were a plasma which flows predominantly upwards, whereas in case of the magnetic topology (i), as if there were a plasma flow expanding from the neutral line. (4) The small-scale physics in the bursting phase represent localized diffusive processes. (5) The local diffusivity usually has a value which is effectively zero, and it turns locally to an anomalous value if a threshold is exceeded, whereby diffusion dominates the quiet evolution (loading). (6) Flares (avalanches) are accompanied by the appearance of localized, intense electric fields. (7) In a variant on the X-CA model, the magnitude of the current is used directly in the instability criterion. First results indicate that the SOC state persists. (8) The current-dissipation during flares is spatially fragmented into a large number of dissipative current-surfaces of varying sizes, which show a highly dynamic temporal evolution.

A set-up is introduced which can be superimposed onto the existing solar flare cellular automata (CA) models, and which specifies the interpretation of the model's variables. It extends the CA models, yielding the magnetic field, the current, and an approximation to the electric field, in a way that is consistent with Maxwell's and the MHD equations. Applications to several solar flare CA models during their natural state (self-organized criticality (SOC)) show, among others, that (1) the magnetic field exhibits characteristic large-scale organization over the entire modeled volume; (2) the magnitude of the current seems spatially dis-organized, with no obvious tendency towards large-scale structures or even local organization; (3) bursts occur at sites with increased current, and after a burst the current is relaxed; (4) by estimating the energy released in individual bursts with the use of the current as Ohmic dissipation, it turns out that the power-law distributions of the released energy persist. The CA models, extended with the set-up, can thus be considered as models for energy-release through current-dissipation. The concepts of power-law loading and anisotropic events (bursts) in CA models are generalized to 3--D vector-field models, and their effect on the magnetic field topology is demonstrated.