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.

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.

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.