We analyze the stochastic acceleration of particles inside a fully developed turbulent plasma. It is well known that large-amplitude magnetic fluctuations and coherent structures in such an environment obey a fractal scaling, and our specific aim is to study for the first time the effects of the fractality of these environments on stochastic acceleration. We have shown that an injected Maxwellian energy distribution is heated and forms a high energy tail in a very short time. Using standard parameters for the low solar corona, the injected Maxwellian distribution of electrons gets heated from the initial 100 eV to 10 KeV, and the power-law index of the high energy tail is about -2.3-4.0. The high energy tail starts around 100 keV, and reaches 10 MeV. The index of the power-law tail depends on the system size, and it is in good agreement with observed values for realistic system sizes. The heating and acceleration process is very fast (\sim 2 s). The reason why the acceleration time is so short is that the particles are trapped within small scale parts of the fractal environment, and their scattering mean free path reduces drastically. The presence of small scale activity also pulls easily particles from the thermal pool, so there is no need for a seed population. The mean square displacement in space and energy is superdiffusive for the high energy particles.

We analyze the transport properties of charged particles (ions and electrons) interacting with randomly formed magnetic scatterers (e.g.\ large scale local ``magnetic fluctuations'' or ``coherent magnetic irregularities'' usually present in strongly turbulent plasmas), using the energization processes proposed initially by Fermi in 1949. The scatterers are formed by large scale local fluctuations ($\delta B/B \approx 1$) and are randomly distributed inside the unstable magnetic topology. We construct a 3D grid on which a small fraction of randomly chosen grid points are acting as scatterers. In particular, we study how a large number of test particles are accelerated and transported inside a collection of scatterers in a finite volume. Our main results are: (1) The spatial mean-square displacement $<(\Delta r)^2>$ inside the stochastic Fermi accelerator is superdiffusive, $<(\Delta r)^2> \sim t^{a_{r}},$ with $a_r \sim 1.2-1.6$, for the high energy electrons with kinetic energy $(W)$ larger than $1 MeV$, and it is normal ($a_r=1$) for the heated low energy $(W< 10 keV)$ electrons. (2) The transport properties of the high energy particles are closely related with the mean-free path that the particles travel in-between the scatterers ($\lambda_{sc}$). The smaller $\lambda_{sc}$ is, the faster the electrons and ions escape from the acceleration volume. (3) The mean displacement in energy $<\Delta W> \sim t^{a_{W}}$ is strongly enhanced inside the acceleration volume $(a_W=1.5- 2.5)$ for the high energy particles compared to the thermal low energy particles ($a_W=0.4$), i.e.\ high energy particles undergo an enhanced systematic gain in energy.(4) The mean-square displacement in energy $<W^2>$ is superdiffusive for the high energy particles and normal for the low energy, heated particles.

The interaction between emerging and pre-existing magnetic fields in the solar atmosphere can trigger several dynamic phenomena, such as eruptions and jets. A key element during this interaction is the formation of large scale current sheets and, eventually, their fragmentation that leads to the creation of a strongly turbulent environment. In this paper, we study the kinetic aspects of the interaction (reconnection) between emerging and ambient magnetic fields. We show that the statistical properties of the spontaneously fragmented and fractal electric fields are responsible for the efficient heating and acceleration of charged particles, which form a power law tail at high energies on sub-second time scales. A fraction of the energized particles escapes from the acceleration volume, with a super-hot component with temperature close to $150\,$MK, and with a power law high energy tail with index between -2 and -3. We estimate the transport coefficients in energy space from the dynamics of the charged particles inside the fragmented and fractal electric fields, and the solution of a fractional transport equation, as appropriate for a strongly turbulent plasma, agrees with the test particle simulations. We also show that the acceleration mechanism is not related to Fermi acceleration, and the Fokker Planck equation is inconsistent and not adequate as a transport model. Finally, we address the problem of correlations between spatial transport and transport in energy space. Our results confirm the observations reported for high energy particles (hard X-rays, type III bursts and solar energetic particles) during the emission of solar jets.

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.

Turbulence, magnetic reconnection, and shocks can be present in explosively unstable plasmas, forming a new electromagnetic environment, which we call here turbulent reconnection, and where spontaneous formation of current sheets takes place. We will show that the heating and the acceleration of particles is the result of the synergy of stochastic (second order Fermi) and systematic (first order Fermi) acceleration inside fully developed turbulence. The solar atmosphere is magnetically coupled to a turbulent driver (the convection zone), therefore the appearance of turbulent reconnection in the solar atmosphere is externally driven. Turbulent reconnection, once it is established in the solar corona, drives the coronal heating and particle acceleration.

Diffusive Shock Acceleration (DSA) cannot efficiently accelerate particles without the presence of self-consistently generated or pre-existing strong turbulence ($ \delta B/B \sim 1 $) in the vicinity of the shock. The problem we address in this article is: if large amplitude magnetic disturbances are present upstream and downstream of a shock then Turbulent Reconnection (TR) will set in and will participate not only in the elastic scattering of particles but also in their heating and acceleration. We demonstrate that large amplitude magnetic disturbances and Unstable Current Sheets (UCS), spontaneously formed in the strong turbulence in the vicinity of a shock, can accelerate particles as efficiently as DSA in {\bf large scale systems and on long time scales}. We start our analysis with "elastic" scatterers upstream and downstream and estimate the energy distribution of particles escaping from the shock, recovering the well known results from the DSA theory. Next we analyze the additional interaction of the particles with active scatterers (magnetic disturbances and UCS) upstream and downstream of the shock. We show that the asymptotic energy distribution of the particles accelerated by DSA/TR has very similar characteristics with the one due to DSA alone, but the synergy of DSA with TR is much more efficient: The acceleration time is an order of magnitude shorter and the maximum energy reached two orders of magnitude higher. We claim that DSA is the dominant acceleration mechanism in a short period before TR is established, and then strong turbulence will dominate the heating and acceleration of the particles. In other words, the shock serves as the mechanism to set up a strongly turbulent environment, in which the acceleration mechanism will ultimately be the synergy of DSA and TR.

The important characteristic of turbulent reconnection is that it combines large scale magnetic disturbances $(\delta B/B \sim 1)$ with randomly distributed Unstable Current Sheets (UCSs). Many well known non linear MHD structures (strong turbulence, current sheet(s), shock(s)) lead asymptotically to the state of turbulent reconnection. We analyze in this article, for the first time, the energization of electrons and ions in a {\bf large scale} environment that {\bf combines} large amplitude disturbances propagating with sub-Alfv\'enic speed with UCSs. The magnetic disturbances interact stochastically (second order Fermi) with the charged particles and they play a crucial role in the heating of the particles, while the UCS interact systematically (first order Fermi) and play a crucial role in the formation of the high energy tail. The synergy of stochastic and systematic acceleration provided by the mixture of magnetic disturbances and UCSs influences the energetics of the thermal and non-thermal particles, the power law index, and the time the particles remain inside the energy release volume. We show that this synergy can explain the observed very fast and impulsive particle acceleration and the slightly delayed formation of a super-hot particle population.

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.

The limitation of the Quasilinear Theory (QLT) to describe the diffusion of electrons and ions in velocity space when interacting with a spectrum of large amplitude electrostatic Langmuir, Upper and Lower hybrid waves, is analyzed. We analytically and numerically estimate the threshold for the amplitude of the waves above which the QLT breaks down, using a test particle code. The evolution of the velocity distribution, the velocity-space diffusion coefficients, the driven current, and the heating of the particles are investigated, for the interaction with small and large amplitude electrostatic waves, i.e. in both regimes, there where QLT is valid and there where it clearly breaks down.

Coronal heating through the explosive release of magnetic energy remains an open problem in solar physics. Recent hydrodynamical models attempt an investigation by placing swarms of 'nanoflares' at random sites and times in modeled one-dimensional coronal loops. We investigate the problem in three dimensions, using extrapolated coronal magnetic fields of observed solar active regions. We apply a nonlinear force-free field extrapolation above an observed photospheric magnetogram of NOAA active region (AR) 11158. We then determine the locations, energy contents, and volumes of 'unstable' areas, namely areas prone to releasing magnetic energy due to locally accumulated electric current density. Statistical distributions of these volumes and their fractal dimension are inferred, investigating also their dependence on spatial resolution. Further adopting a simple resistivity model, we infer the properties of the fractally distributed electric fields in these volumes. Next, we monitor the evolution of 10^5 particles (electrons and ions) obeying an initial Maxwellian distribution with a temperature of 10 eV, by following their trajectories and energization when subjected to the resulting electric fields. For computational convenience, the length element of the magnetic-field extrapolation is 1 arcsec, much coarser than the particles collisional mean free path in the low corona. The presence of collisions traps the bulk of the plasma around the unstable volumes, or current sheets (UCS), with only a tail of the distribution gaining substantial energy. Assuming that the distance between UCS is similar to the collisional mean free path we find that the low active-region corona is heated to 100-200 eV, corresponding to temperatures exceeding 2 MK, within tens of seconds for electrons and thousands of seconds for ions. Fractally distributed, nanoflare-triggening fragmented UCS ...

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.

In this review many of the well known tools for the analysis of Complex systems are used in order to study the global coupling of the turbulent convection zone with the solar atmosphere where the magnetic energy is dissipated explosively. Several well documented observations are not easy to interpret with the use of Magnetohydrodynamic (MHD) and/or Kinetic numerical codes. Such observations are: (1) The size distribution of the Active Regions (AR) on the solar surface, (2) The fractal and multi fractal characteristics of the observed magnetograms, (3) The Self-Organised characteristics of the explosive magnetic energy release and (4) the very efficient acceleration of particles during the flaring periods in the solar corona. We review briefly the work published the last twenty five years on the above issues and propose solutions by using methods borrowed from the analysis of complex systems. The scenario which emerged is as follows: (a) The fully developed turbulence in the convection zone generates and transports magnetic flux tubes to the solar surface. Using probabilistic percolation models we were able to reproduce the size distribution and the fractal properties of the emerged and randomly moving magnetic flux tubes. (b) Using a Non Linear Force Free (NLFF) magnetic extrapolation numerical code we can explore how the emerged magnetic flux tubes interact nonlinearly and form thin and Unstable Current Sheets (UCS) inside the coronal part of the AR. (c) The fragmentation of the UCS and the redistribution of the magnetic field locally, when the local current exceeds a Critical threshold, is a key process which drives avalanches and forms coherent structures. This local reorganization of the magnetic field, enhances the energy dissipation and influences the global evolution of the complex magnetic topology.

We interpret solar flares as events originating from active regions that have reached the Self Organized Critical state, by using a refined Cellular Automaton model with initial conditions derived from observations. Aims: We investigate whether the system, with its imposed physical elements,reaches a Self Organized Critical state and whether well-known statistical properties of flares, such as scaling laws observed in the distribution functions of characteristic parameters, are reproduced after this state has been reached. Results: Our results show that Self Organized Criticality is indeed reached when applying specific loading and relaxation rules. Power law indices obtained from the distribution functions of the modeled flaring events are in good agreement with observations. Single power laws (peak and total flare energy) as well as power laws with exponential cutoff and double power laws (flare duration) are obtained. The results are also compared with observational X-ray data from GOES satellite for our active-region sample. Conclusions: We conclude that well-known statistical properties of flares are reproduced after the system has reached Self Organized Criticality. A significant enhancement of our refined Cellular Automaton model is that it commences the simulation from observed vector magnetograms, thus facilitating energy calculation in physical units. The model described in this study remains consistent with fundamental physical requirements, and imposes physically meaningful driving and redistribution rules.

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.

We study the interaction of a gravitational wave (GW) with a plasma that is strongly magnetized. The GW is considered a small disturbance, and the plasma is modeled by the general relativistic analogue of the induction equation of ideal MHD and the single fluid equations. The equations are derived without neglecting any of the non-linear interaction terms, and the non-linear equations are integrated numerically. We find that for strong magnetic fields of the order of $10^{15} $G the GW excites electromagnetic plasma waves very close to the magnetosonic mode. The magnetic and electric field oscillations have very high amplitude, and a large amount of energy is absorbed from the GW by the electromagnetic oscillations, of the order of $10^{23} $erg/cm$^3$ in the case presented here. The absorbed energy is proportional to $B_0^2$, with $B_0$ the background magnetic field. The energization of the plasma takes place on fast time scales of the order of milliseconds. The amount of absorbed energy is comparable to the energies emitted in the most energetic astrophysical events, such as giant flares on magnetars and possibly even short Gamma ray bursts (GRB), for which the mechanism analyzed here also has the fast time-scales required.

Electron and proton acceleration in three-dimensional electric and magnetic fields is studied through test particle simulations. The fields are obtained by a three-dimensional magnetohydrodynamic simulation of magnetic reconnection in slab geometry. The nonlinear evolution of the system is characterized by the growth of many unstable modes and the initial current sheet is fragmented with formation of small scale structures. We inject at random points inside the evolving current sheet a Maxwellian distribution of particles. In relatively short time (less than a millisecond) the particles develop a power law tail. The acceleration is extremely efficient and the electrons absorb a large percentage of the available energy in a small fraction of the characteristic time of the MHD simulation, suggesting that resistive MHD codes, used extensively in the current literature, are unable to represent the full extent of particle acceleration in 3D reconnection.

We propose a model for the acceleration of charged particles in interplanetary space that appear during quiet time periods, that is, not associated with solar activity events like intense flares or coronal mass ejections. The interaction of charged particles with modeled turbulent electromagnetic fields, which mimic the fields observed in the interplanetary medium, is studied. The turbulence is modeled by means of a dynamical system,the Gledzer-Ohkitani-Yamada (GOY) shell model, which describes the gross features of the Navier-Stokes equations. The GOY model is used to build a 3-D velocity field, which in turn is used to numerically solve the ideal magneto-hydrodynamic (MHD) induction equation, while the electric field is calculated from the ideal Ohm's law. Particle acceleration in such an environment is investigated by test particle simulations, and the resulting energy distributions are discussed and compared to observations of suprathermal electrons and ions during quiet periods in interplanetary space.

We study the acceleration of electrons and protons interacting with localized, multiple, small-scale dissipation regions inside an evolving, turbulent active region. The dissipation regions are Unstable Current Sheets (UCS), and in their ensemble they form a complex, fractal, evolving network of acceleration centers. Acceleration and energy dissipation are thus assumed to be fragmented. A large-scale magnetic topology provides the connectivity between the UCS and determines in this way the degree of possible multiple acceleration. The particles travel along the magnetic field freely without loosing or gaining energy, till they reach a UCS. In a UCS, a variety of acceleration mechanisms are active, with the end-result that the particles depart with a new momentum. The stochastic acceleration process is represented in the form of Continuous Time Random Walk (CTRW), which allows to estimate the evolution of the energy distribution of the particles. It is found that under certain conditions electrons are heated and accelerated to energies above 1 MeV in much less than a second. Hard X-ray (HXR) and microwave spectra are calculated from the electrons' energy distributions, and they are found to be compatible with the observations. Ions (protons) are also heated and accelerated, reaching energies up to 10 MeV almost simultaneously with the electrons. The diffusion of the particles inside the active region is extremely fast (anomalous super-diffusion). Although our approach does not provide insight into the details of the specific acceleration mechanisms involved, its benefits are that it relates acceleration to the energy release, and it well describes the stochastic nature of the acceleration process.

The formation and evolution of active regions is an inherently complex phenomenon. Magnetic fields generated at the base of the convection zone follow a chaotic evolution before reaching the solar surface. In this article, we use a 2-D probabilistic Cellular Automaton (CA) to model the statistical properties of the magnetic patterns formed on the solar surface and to estimate the magnetic energy released in the interaction of opposite polarities. We assume that newly emerged magnetic flux tubes stimulate the emergence of new magnetic flux in their neighborhood. The flux-tubes move randomly on the surface of the sun, and they cancel and release their magnetic energy when they collide with magnetic flux of opposite polarity, or diffuse into the ``empty'' photosphere. We assume that cancellation of magnetic flux in collisions causes "flares" and determine the released energy as the difference in the square of the magnetic field flux. The statistics of the simulated "flares" follow a power-law distribution in energy, with power-law index a=2.2 +/- 0.1. The size distribution function of the simulated active regions exhibits a power law behavior with index k = 1.93 +/- 0.08, and the fractal dimension of the magnetized areas on the simulated solar surface is close to D = 1.42 +/- 0.12. Both quantities, D and k, are inside the range of the observed values.

Narrow-band radio spikes have been recorded during a solar flare with unprecedented resolution. This unique example allows to study the effect of low resolution in previously published peak-flux distributions of radio spikes. We give a general, analytical expression for how an actual peak-flux distribution is changed in shape if the peaks are determined with low temporal and/or frequency resolution. It turns out that, generally, low resolution tends to cause an exponential behavior at large flux values if the actual distribution is of power-law shape. The distribution may be severely altered if the burst-duration depends on the peak-flux. The derived expression is applicable also to peak-flux distributions derived at other wavelengths (e.g. soft and hard X-rays, EUV). We show that for the analyzed spike-event the resolution was sufficient for a reliable peak flux distribution. It can be fitted by generalized power-laws or by an exponential.

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.