Understanding how the release of stored magnetic energy contributes to the generation of non-thermal high energy particles during solar flares is an important open problem in solar physics. Magnetic reconnection plays a fundamental role in the energy release and conversion processes taking place during flares. A common approach for investigating particle acceleration is to use test particles in fields derived from magnetohydrodynamic (MHD) simulations of reconnection. These MHD simulations use anomalous resistivities that are much larger than the Spitzer resistivity based on Coulomb collisions. The processes leading to enhanced resistivity should also affect the test particles. We explore the link between resistivity and particle orbits building on a previous study using a 2D MHD simulation of magnetic reconnection. This paper extends the previous investigation to a 3D magnetic reconnection configuration and to study the effect on test particle orbits. We carried out orbit calculations using a 3D MHD simulation of separator reconnection. We use the relativistic guiding centre approximation including stochastic pitch angle scattering. The effects of varying the resistivity and the models for pitch angle scattering on particle orbit trajectories, final positions, energy spectra, final pitch angle distribution, and orbit duration are all studied in detail. Pitch angle scattering widens collimated beams of orbit trajectories, allowing orbits to access previously unaccessible field lines; this causes final positions to spread to topological structures that were previously inaccessible. Scattered orbit energy spectra are found to be predominantly affected by the level of anomalous resistivity, with the pitch angle scattering model only playing a role in isolated cases. Scattering is found to play a crucial role in determining the pitch angle and orbit duration distributions.

Zipper reconnection has been proposed as a mechanism for creating most of the twist in the flux tubes that are present prior to eruptive flares and coronal mass ejections. We have conducted a first numerical experiment on this new regime of reconnection, where two initially untwisted parallel flux tubes are sheared and reconnected to form a large flux rope. We describe the properties of this experiment, including the linkage of magnetic flux between concentrated flux sources at the base of the simulation, the twist of the newly formed flux rope and the conversion of mutual magnetic helicity in the sheared pre-reconnection state into the self-helicity of the newly formed flux rope.

The aim of this work is to investigate and characterise non-thermal particle behaviour in a three-dimensional (3D) magnetohydrodynamical (MHD) model of unstable multi-threaded flaring coronal loops. We have used a numerical scheme which solves the relativistic guiding centre approximation to study the motion of electrons and protons. The scheme uses snapshots from high resolution numerical MHD simulations of coronal loops containing two threads, where a single thread becomes unstable and (in one case) destabilises and merges with an additional thread. The particle responses to the reconnection and fragmentation in MHD simulations of two loop threads are examined in detail. We illustrate the role played by uniform background resistivity and distinguish this from the role of anomalous resistivity using orbits in an MHD simulation where only one thread becomes unstable without destabilising further loop threads. We examine the (scalable) orbit energy gains and final positions recovered at different stages of a second MHD simulation wherein a secondary loop thread is destabilised by (and merges with) the first thread. We compare these results with other theoretical particle acceleration models in the context of observed energetic particle populations during solar flares.

The conversion of magnetic energy into other forms during solar flares is one of the outstanding open problems in solar physics. It is generally accepted that magnetic reconnection plays a crucial role in these conversion processes. To achieve the rapid energy release required in solar flares, an anomalous resistivity, orders of magnitude higher than the Spitzer resistivity, is often used in MHD simulations of reconnection. Spitzer resistivity is based on Coulomb scattering, which becomes negligible at the high energies achieved by accelerated particles. As a result, simulations of particle acceleration in reconnection events are often performed in the absence of any interaction between accelerated particles and any background plasma. This need not be the case for scattering associated with anomalous resistivity caused by turbulence within solar flares, as the higher resistivity implies an elevated scattering rate. We present results of test particle calculations, with and without pitch angle scattering, subject to fields derived from MHD simulations of two-dimensional (2D) X-point reconnection. Scattering rates proportional to the ratio of the anomalous resistivity to the local Spitzer resistivity, as well as at fixed values, are considered. Pitch angle scattering, which is independent of the anomalous resistivity, causes higher maximum energies in comparison to those obtained without scattering. Scattering rates which are dependent on the local anomalous resistivity tend to produce fewer highly energised particles due to weaker scattering in the separatrices, even though scattering in the current sheet may be stronger when compared to resistivity-independent scattering. Strong scattering also causes an increase in the number of particles exiting the computational box in the reconnection outflow region, as opposed to along the separatrices as is the case in the absence of scattering.

Collapsing magnetic traps (CMTs) are one proposed mechanism for generating non-thermal particle populations in solar flares. CMTs occur if an initially stretched magnetic field structure relaxes rapidly into a lower-energy configuration, which is believed to happen as a by-product of magnetic reconnection. A similar mechanism for energising particles has also been found to operate in the Earth's magnetotail. One particular feature proposed to be of importance for particle acceleration in the magnetotail is that of a braking plasma jet, i.e. a localised region of strong flow encountering stronger magnetic field which causes the jet to slow down and stop. Such a feature has not been included in previously proposed analytical models of CMTs for solar flares. In this work we incorporate a braking plasma jet into a well studied CMT model for the first time. We present results of test particle calculations in this new CMT model. We observe and characterise new types of particle behaviour caused by the magnetic structure of the jet braking region, which allows electrons to be trapped both in the braking jet region and the loop legs. We compare and contrast the behaviour of particle orbits for various parameter regimes of the underlying trap by examining particle trajectories, energy gains and the frequency with which different types of particle orbit are found for each parameter regime.

The aim of this work is to investigate and characterise particle behaviour in an (observationally-driven) 3D MHD model of the solar atmosphere above a slowly evolving, non-flaring active region. We use a relativistic guiding-centre particle code to investigate the behaviour of selected particle orbits, distributed throughout a single snapshot of the 3D MHD simulation. Two distinct particle acceleration behaviours are recovered, which affect both electrons and protons: (i) direct acceleration along field lines and (ii) tangential drifting of guiding centres with respect to local magnetic field. However, up to 40\% of all particles actually experience a form of (high energy) particle trap, because of changes in the direction of the electric field and unrelated to the strength of the magnetic field; such particles are included in the first category. Additionally, category (i) electron and proton orbits undergo surprisingly strong acceleration to non-thermal energies ($\lesssim42$MeV), because of the strength and extent of super-Dreicer electric fields created by the MHD simulation. Reducing the electric field strength of the MHD model does not significantly affect the efficiency of the (electric field-based) trapping mechanism, but does reduce the peak energies gained by orbits. We discuss the implications for future experiments, which aim to simulate non-flaring active region heating and reconnection.

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.

Magnetic fields permeate the entire solar atmosphere weaving an extremely complex pattern on both local and global scales. In order to understand the nature of this tangled web of magnetic fields, its magnetic skeleton, which forms the boundaries between topologically distinct flux domains, may be determined. The magnetic skeleton consists of null points, separatrix surfaces, spines and separators. The skeleton is often used to clearly visualize key elements of the magnetic configuration, but parts of the skeleton are also locations where currents and waves may collect and dissipate. In this review, the nature of the magnetic skeleton on both global and local scales, over solar cycle time scales, is explained. The behaviour of wave pulses in the vicinity of both nulls and separators is discussed and so too is the formation of current layers and reconnection at the same features. Each of these processes leads to heating of the solar atmosphere, but collectively do they provide enough heat, spread over a wide enough area, to explain the energy losses throughout the solar atmosphere? Here, we consider this question for the three different solar regions: active regions, open-field regions and the quiet Sun. We find that the heating of active regions and open-field regions is highly unlikely to be due to reconnection or wave dissipation at topological features, but it is possible that these may play a role in the heating of the quiet Sun. In active regions, the absence of a complex topology may play an important role in allowing large energies to build up and then, subsequently, be explosively released in the form of a solar flare. Additionally, knowledge of the intricate boundaries of open-field regions (which the magnetic skeleton provides) could be very important in determining the main acceleration mechanism(s) of the solar wind.

Recent observations have demonstrated that waves which are capable of carrying large amounts of energy are ubiquitous throughout the solar corona. However, the question of how this wave energy is dissipated (on which time and length scales) and released into the plasma remains largely unanswered. Both analytic and numerical models have previously shown that Alfv\'enic turbulence may play a key role not only in the generation of the fast solar wind, but in the heating of coronal loops. In an effort to bridge the gap between theory and observations, we expand on a recent study [De Moortel et al., ApJL, 782:L34, 2014] by analyzing thirty-seven clearly isolated coronal loops using data from the Coronal Multi-channel Polarimeter (CoMP) instrument. We observe Alfv\'enic perturbations with phase speeds which range from 250-750 km/s and periods from 140-270 s for the chosen loops. While excesses of high frequency wave-power are observed near the apex of some loops (tentatively supporting the onset of Alfv\'enic turbulence), we show that this excess depends on loop length and the wavelength of the observed oscillations. In deriving a proportional relationship between the loop length/wavelength ratio and the enhanced wave power at the loop apex, and from the analysis of the line-widths associated with these loops, our findings are supportive of the existence of Alfv\'enic turbulence in coronal loops.

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 presence of turbulent phenomena in the outer solar atmosphere is a given. However, because we are reduced to remotely sensing the atmosphere of a star with instruments of limited spatial and/or spectral resolution, we can only infer the physical progression from macroscopic to microscopic phenomena. Even so, we know that many, if not all, of the turbulent phenomena that pervade interplanetary space have physical origins at the Sun and so in this brief article we consider some recent measurements which point to sustained potential source(s) of heliospheric turbulence in the magnetic and thermal domains. In particular, we look at the scales of magnetism that are imprinted on the outer solar atmosphere by the relentless magneto-convection of the solar interior and combine state-of-the-art observations from the Solar Dynamics Observatory (SDO) and the Coronal Multi-channel Polarimeter (CoMP) which are beginning to hint at the origins of the wave/plasma interplay prevalent closer to the Earth. While linking these disparate scales of observation and understanding of their connection is near to impossible, it is clear that the constant evolution of subsurface magnetism on a host of scales guides and governs the flow of mass and energy at the smallest scales. In the near future significant progress in this area will be made by linking observations from high resolution platforms like the Interface Region Imaging Spectrograph (IRIS) and Advanced Technology Solar Telescope (ATST) with full-disk synoptic observations such as those presented herein.

Waves have long been thought to contribute to the heating of the solar corona and the generation of the solar wind. Recent observations have demonstrated evidence of quasi-periodic longitudinal disturbances and ubiquitous transverse wave propagation in many different coronal environments. This paper investigates signatures of different types of oscillatory behaviour, both above the solar limb and on-disk, by comparing findings from the Coronal Multi-channel Polarimeter (CoMP) and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO) for the same active region. We study both transverse and longitudinal motion by comparing and contrasting time-distance images of parallel and perpendicular cuts along/across active region fan loops. Comparisons between parallel space-time features in CoMP Doppler velocity and transverse oscillations in AIA images are made, together with space-time analysis of propagating quasi-periodic intensity features seen near the base of loops in AIA. Signatures of transverse motions are observed along the same magnetic structure using CoMP Doppler velocity (Vphase=600-750km/s, P=3-6mins) and in AIA/SDO above the limb (P=3-8mins). Quasi-periodic intensity features (Vphase=100-200km/s, P=6-11mins) also travel along the base of the same structure. On the disk, signatures of both transverse and longitudinal intensity features were observed by AIA; both show similar properties to signatures found along structures anchored in the same active region three days earlier above the limb. Correlated features are recovered by space-time analysis of neighbouring tracks over perpendicular distances of <2.6Mm.

The highly dynamical, complex nature of the solar atmosphere naturally implies the presence of waves in a topologically varied magnetic environment. Here, the interaction of waves with topological features such as null points is inevitable and potentially important for energetics. The low resistivity of the solar coronal plasma implies that non-MHD effects should be considered in studies of magnetic energy release in this environment. This paper investigates the role of the Hall term in the propagation and dissipation of waves, their interaction with 2D magnetic X-points and the nature of the resulting reconnection. A Lagrangian remap shock-capturing code (Lare2d) is used to study the evolution of an initial fast magnetoacoustic wave annulus for a range of values of the ion skin depth in resistive Hall MHD. A magnetic null-point finding algorithm is also used to locate and track the evolution of the multiple null-points that are formed in the system. Depending on the ratio of ion skin depth to system size, our model demonstrates that Hall effects can play a key role in the wave-null interaction. In particular, the initial fast-wave pulse now consists of whistler and ion-cyclotron components; the dispersive nature of the whistler wave leads to (i) earlier interaction with the null, (ii) the creation of multiple additional, transient nulls and, hence, an increased number of energy release sites. In the Hall regime, the relevant timescales (such as the onset of reconnection and the period of the oscillatory relaxation) of the system are reduced significantly, and the reconnection rate is enhanced.

Aims. To determine the effect of the Hall term in the generalised Ohm's law on the damping and phase mixing of Alfven waves in the ion cyclotron range of frequencies in uniform and non-uniform equilibrium plasmas. Methods. Wave damping in a uniform plasma is treated analytically, whilst a Lagrangian remap code (Lare2d) is used to study Hall effects on damping and phase mixing in the presence of an equilibrium density gradient. Results. The magnetic energy associated with an initially Gaussian field perturbation in a uniform resistive plasma is shown to decay algebraically at a rate that is unaffected by the Hall term to leading order in k^2di^2 where k is wavenumber and di is ion skin depth. A similar algebraic decay law applies to whistler perturbations in the limit k^2di^2>>1. In a non-uniform plasma it is found that the spatially-integrated damping rate due to phase mixing is lower in Hall MHD than it is in MHD, but the reduction in the damping rate, which can be attributed to the effects of wave dispersion, tends to zero in both the weak and strong phase mixing limits.