No Arabic abstract
While the exact acceleration mechanism of energetic particles during solar flares is (as yet) unknown, magnetic reconnection plays a key role both in the release of stored magnetic energy of the solar corona and the magnetic restructuring during a flare. Recent work has shown that special field lines, called separators, are common sites of reconnection in 3D numerical experiments. To date, 3D separator reconnection sites have received little attention as particle accelerators. We investigate the effectiveness of separator reconnection as a particle acceleration mechanism for electrons and protons. We study the particle acceleration using a relativistic guiding-centre particle code in a time-dependent kinematic model of magnetic reconnection at a separator. The effect upon particle behaviour of initial position, pitch angle and initial kinetic energy are examined in detail, both for specific (single) particle examples and for large distributions of initial conditions. The separator reconnection model contains several free parameters and we study the effect of changing these parameters upon particle acceleration, in particular in view of the final particle energy ranges which agree with observed energy spectra.
The aim of this work is to investigate and characterise particle behaviour in a 3D magnetohydrodynamic (MHD) model of a reconnecting magnetic separator. We use a relativistic guiding-centre test-particle code to investigate electron and proton acceleration in snapshots from 3D MHD separator reconnection experiments, and compare the results with findings from an analytical separator reconnection model studied in a previous investigation. The behaviour and acceleration of large distributions of particles are examined in detail for both analytical and numerical separator reconnection models. Accelerated particle orbit trajectories are shown to follow the separator before leaving the system along the separatrix surface of one of the nulls (determined by particle species) in the system of both models. A sufficiently localised electric field about the separator causes the orbits to appear to follow the spine bounding the separatrix surface field lines instead. We analyse and discuss the locations and spread of accelerated particle orbit final positions, which are seen to change over time in the numerical separator reconnection model. We deduce a simple relationship between the final energy range of particle orbits and the model dimensions, and discuss its implications for observed magnetic separators in the solar corona.
Twisted coronal loops should be ubiquitous in the solar corona. Twisted magnetic fields contain excess magnetic energy, which can be released during magnetic reconnection, causing solar flares. The aim of this work is to investigate magnetic reconnection, and particle acceleration and transport in kink-unstable twisted coronal loops, with a focus on the effects of resistivity, loop geometry and atmospheric stratification. Another aim is to perform forward-modelling of bremsstrahlung emission and determine the structure of hard X-ray sources. We use a combination of magnetohydrodynamic (MHD) and test-particle methods. First, the evolution of the kinking coronal loop is considered using resistive MHD model, incorporating atmospheric stratification and loop curvature. Then, the obtained electric and magnetic fields and density distributions are used to calculate electron and proton trajectories using a guiding-centre approximation, taking into account Coulomb collisions. It is shown that electric fields in twisted coronal loops can effectively accelerate protons and electrons to energies up to 10 MeV. High-energy particles have hard, nearly power-law energy spectra. The volume occupied by high-energy particles demonstrates radial expansion, which results in the expansion of the visible hard X-ray loop and a gradual increase in hard X-ray footpoint area. Synthesised hard X-ray emission reveals strong footpoint sources and the extended coronal source, whose intensity strongly depends on the coronal loop density.
Using analytical and numerical methods (fluid and particle-in-cell simulations) we study a number of model problems involving merger of magnetic flux tubes in relativistic magnetically-dominated plasma. Mergers of current-carrying flux tubes (exemplified by the two dimensional `ABC structures) and zero total current magnetic flux tubes are considered. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration during flux mergers: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. During the X-point collapse particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $sigma$. For plasma magnetization $sigma leq 10^2$ the spectrum power law index is $p> 2$; in this case the maximal energy depends linearly on the size of the reconnecting islands. For higher magnetization, $sigma geq 10^2$, the spectra are hard, $p< 2$, yet the maximal energy $gamma_{max}$ can still exceed the average magnetic energy per particle, $ sim sigma$, by orders of magnitude (if $p$ is not too close to unity). The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced magnetic reconnection, further accelerating the particles, but proceeds at a slower reconnection rate.
Using particle-in-cell (PIC) numerical simulations with electron-positron pair plasma, we study how the efficiencies of magnetic dissipation and particle acceleration scale with the initial coherence length $lambda_0$ in relation to the system size $L$ of the two-dimensional (2D) `Arnold-Beltrami-Childress (ABC) magnetic field configurations. Topological constraints on the distribution of magnetic helicity in 2D systems, identified earlier in relativistic force-free (FF) simulations, that prevent the high-$(L/lambda_0)$ configurations from reaching the Taylor state, limit the magnetic dissipation efficiency to about $epsilon_{rm diss} simeq 60%$. We find that the peak growth time scale of the electric energy $tau_{rm E,peak}$ scales with the characteristic value of initial Alfven velocity $beta_{rm A,ini}$ like $tau_{rm E,peak} propto (lambda_0/L)beta_{rm A,ini}^{-3}$. The particle energy change is decomposed into non-thermal and thermal parts, with non-thermal energy gain dominant only for high initial magnetisation. The most robust description of the non-thermal high-energy part of the particle distribution is that the power-law index is a linear function of the initial magnetic energy fraction.
Magnetic reconnection, especially in the relativistic regime, provides an efficient mechanism for accelerating relativistic particles and thus offers an attractive physical explanation for nonthermal high-energy emission from various astrophysical sources. I present a simple analytical model that elucidates key physical processes responsible for reconnection-driven relativistic nonthermal particle acceleration (NTPA) in the large-system, plasmoid-dominated regime in two dimensions. The model aims to explain the numerically-observed dependencies of the power-law index $p$ and high-energy cutoff $gamma_c$ of the resulting nonthermal particle energy spectrum $f(gamma)$ on the ambient plasma magnetization $sigma$, and (for $gamma_c$) on the system size $L$. In this self-similar model, energetic particles are continuously accelerated by the out-of-plane reconnection electric field $E_{rm rec}$ until they become magnetized by the reconnected magnetic field and eventually trapped in plasmoids large enough to confine them. The model also includes diffusive Fermi acceleration by particle bouncing off rapidly moving plasmoids. I argue that the balance between electric acceleration and magnetization controls the power-law index, while trapping in plasmoids governs the cutoff, thus tying the particle energy spectrum to the plasmoid distribution.