We model the time dependent radio emission from a disk accretion event in a T-Tauri star using 3D, ideal magnetohydrodynamic simulations combined with a gyrosynchrotron emission and radiative transfer model. We predict for the first time, the multi-frequency (1$-$1000 GHz) intensity and circular polarisation from a flaring T-Tauri star. A flux tube, connecting the star with its circumstellar disk, is populated with a distribution of non-thermal electrons which is allowed to decay exponentially after a heating event in the disk and the system is allowed to evolve. The energy distribution of the electrons, as well as the non-thermal power law index and loss rate, are varied to see their effect on the overall flux. Spectra are generated from different lines of sight, giving different views of the flux tube and disk. The peak flux typically occurs around 20$-$30 GHz and the radio luminosity is consistent with that observed from T-Tauri stars. For all simulations, the peak flux is found to decrease and move to lower frequencies with elapsing time. The frequency-dependent circular polarisation can reach 10$-$30$\%$ but has a complex structure which evolves as the flare evolves. Our models show that observations of the evolution of the spectrum and its polarisation can provide important constraints on physical properties of the flaring environment and associated accretion event.

The goal of this study is to explore a novel method for the solar photospheric magnetic field diagnostics using Stokes V widths of different magnetosensitive Fe~I spectral lines. We calculate Stokes I and V profiles of several Fe I lines based on a one-dimensional photospheric model VAL C using the NICOLE radiative transfer code. These profiles are used to produce calibration curves linking the intrinsic magnetic field values with the widths of blue peaks of Stokes V profiles. The obtained calibration curves are then tested using the Stokes profiles calculated for more realistic photospheric models based on MHD models of magneto-convection. It is shown that the developed Stokes V widths (SVW) method can be used with various optical and near-infrared lines. Out of six lines considered in this study, FeI 6301 line appears to be the most effective: it is sensitive to fields over ~200G and does not show any saturation up to ~2kG. Other lines considered can also be used for the photospheric field diagnostics with this method, however, only in narrower field value ranges, typically from about 100G to 700-1000G. The developed method can be a useful alternative to the classical magnetic line ratio method, particularly when the choice of lines is limited.

The observed properties (i.e., source size, source position, time duration, decay time) of solar radio emission produced through plasma processes near the local plasma frequency, and hence the interpretation of solar radio bursts, are strongly influenced by propagation effects in the inhomogeneous turbulent solar corona. In this work, a 3D stochastic description of the propagation process is presented, based on the Fokker-Planck and Langevin equations of radio-wave transport in a medium containing anisotropic electron density fluctuations. Using a numerical treatment based on this model, we investigate the characteristic source sizes and burst decay times for Type III solar radio bursts. Comparison of the simulations with the observations of solar radio bursts shows that predominantly perpendicular density fluctuations in the solar corona are required, with an anisotropy factor $\sim 0.3$ for sources observed at around 30~MHz. The simulations also demonstrate that the photons are isotropized near the region of primary emission, but the waves are then focused by large-scale refraction, leading to plasma radio emission directivity that is characterized by a half-width-half-maximum of about 40~degrees near 30~MHz. The results are applicable to various solar radio bursts produced via plasma emission.

Forced magnetic reconnection, a reconnection event triggered by external perturbation, should be ubiquitous in the solar corona. Energy released during such cases can be much greater than that which was introduced by the perturbation. It is unclear how the properties of the external perturbation and the initial current sheet affect the reconnection region properties, and thereby the reconnection dynamics and energy release profile. We investigate the effect of the form of the external perturbation and initial current sheet on the evolution of the reconnection region and the energy release process. Chiefly we explore the non-linear interactions between multiple, simultaneous perturbations, which represent more realistic scenarios. Future work will use these results in test particle simulations to investigate particle acceleration over multiple reconnection events. Simulations are performed using Lare2d, a 2.5D Lagrangian-remap solver for the visco-resistive MHD equations. The model of forced reconnection is extended to include superpositions of sinusoidal driving disturbances, including localised Gaussian perturbations. A transient perturbation is applied to the boundaries of a region containing a force-free current sheet. The simulation domain is sufficiently wide to allow multiple magnetic islands to form and coalesce. Island coalescence contributes significantly to energy release and involves rapid reconnection. Long wavelength modes in perturbations dominate the evolution, without the presence of which reconnection is either slow, as in the case of short wavelength modes, or the initial current sheet remains stable, as in the case of noise perturbations. Multiple perturbations combine in a highly non-linear manner: reconnection is typically faster than when either disturbance is applied individually, with multiple low-energy events contributing to the same total energy release.

Solar flares are explosive events in the solar corona, representing fast conversion of magnetic energy into thermal and kinetic energy, and hence radiation, due to magnetic reconnection. Modelling is essential for understanding and predicting these events. However, self-consistent modelling is extremely difficult due to the vast spatial and temporal scale separation between processes involving thermal plasma (normally considered using magnetohydrodynamic (MHD) approach) and non-thermal plasma (requiring a kinetic approach). In this mini-review we consider different approaches aimed at bridging the gap between fluid and kinetic modelling of solar flares. Two types of approaches are discussed: combined MHD/test-particle (MHDTP) models, which can be used for modelling the flaring corona with relatively small numbers of energetic particles, and hybrid fluid-kinetic methods, which can be used for modelling stronger events with higher numbers of energetic particles. Two specific examples are discussed in more detail: MHDTP models of magnetic reconnection and particle acceleration in kink-unstable twisted coronal loops, and a novel reduced-kinetic model of particle transport.

Early magnetographic observations indicated that magnetic field in the solar photosphere has unresolved small-scale structure. Near-infrared and optical data with extremely high spatial resolution show that these structures have scales of few tens of kilometres, which are not resolved in the majority of solar observations. The goal of this study is to establish the effect of unresolved photospheric magnetic field structure on Stokes profiles observed with relatively low spatial resolution. Ultimately, we aim to develop methods for fast estimation of the photospheric magnetic filling factor and line-of-sight gradient of the photospheric magnetic field, which can be applied to large observational data sets. We exploit 3D MHD models of magneto-convection developed using MURAM code. Corresponding profiles of Fe I 6301.5 and 6302.5 $\mathrm{\AA}$ spectral lines are calculated using NICOLE radiative transfer code. The resulting I and V Stokes [x,y,$\lambda$] cubes with reduced spatial resolution of 150 km are used to calculate magnetic field values as they would be obtained in observations with Hinode/SOT or SDO/HMI. Three different methods of the magnetic filling factor estimation are considered: the magnetic line ratio method, Stokes V width method and a simple statistical method. We find that the statistical method and the Stokes V width method are sufficiently reliable for fast filling factor estimations. Furthermore, we find that Stokes $I\pm V$ bisector splitting gradient can be used for fast estimation of line-of-sight gradient of the photospheric magnetic field.

Magnetic reconnection is believed to be the primary mechanism by which non-potential energy stored in coronal magnetic fields is rapidly released during solar eruptive events. Unfortunately, owing to the small spatial scales on which reconnection is thought to occur, it is not directly observable in the solar corona. However, larger scale processes, such as associated inflow and outflow, and signatures of accelerated particles have been put forward as evidence of reconnection. We explore the origin of a persistent Type I radio source that accompanies a coronal X--shaped structure during its passage across the disk. Of particular interest is the time range around a partial collapse of the structure that is associated with inflow, outflow, and signatures of particle acceleration. Imaging radio observations from the NRH were used to localise the radio source. SDO AIA EUV observations from the same time period were analysed, looking for evidence of inflows and outflows. Further \texttt{mpole} magnetic reconstructions using SDO HMI observations allowed the magnetic connectivity associated with the radio source to be determined. The Type I radio source was well aligned with a magnetic separator identified in the extrapolations. During the partial collapse, gradual (1 km/s) and fast (5 km/s) inflow phases and fast (30 km/s) and rapid (80-100 km/s) outflow phases were observed, resulting in an estimated reconnection rate of ~ 0.06. The radio source brightening and dimming was found to be co-temporal with increased soft x-ray emission in both RHESSI and GOES. We interpret the brightening and dimming of the radio emission as evidence for accelerated electrons in the reconnection region responding to a gradual fall and rapid rise in electric drift velocity, in response to the inflowing and outflowing field lines. These results present a comprehensive example of 3D null-point reconnection.

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.

Four different approximate approaches used to model the stressing of coronal magnetic fields due to an imposed photospheric motion are compared with each other and the results from a full time-dependent magnetohydrodynamic (MHD) code. The assumptions used for each of the approximate methods are tested by considering large photospheric footpoint displacements. We consider a simple model problem, comparing the full nonlinear magnetohydrodynamic evolution, determined with the Lare2D numerical code, with four approximate approaches. Two of these, magneto-frictional relaxation and a quasi-1D Grad-Shafranov approach, assume sequences of equilibria, whilst the other two methods, a second-order linearisation of the MHD equations and Reduced MHD, are time-dependent. The relaxation method is very accurate compared to full MHD for force-free equilibria for all footpoint displacements but has significant errors when the plasma $\beta_0$ is of order unity. The 1D approach gives an extremely accurate description of the equilibria away from the photospheric boundary layers, and agrees well with Lare2D for all parameter values tested. The linearised MHD equations correctly predict the existence of photospheric boundary layers that are present in the full MHD results. As soon as the footpoint displacement becomes a significant fraction of the loop length, the RMHD method fails to model the sequences of equilibria correctly. The full numerical solution is interesting in its own right, and care must be taken for low $\beta_0$ plasmas if the viscosity is too large.

Magnetic reconnection and particle acceleration due to the kink instability in twisted coronal loops can be a viable scenario for confined solar flares. Detailed investigation of this phenomenon requires reliable methods for observational detection of magnetic twist in solar flares, which may not be possible solely through extreme UV and soft X-ray thermal emission. Polarisation of microwave emission in flaring loops can be used as one of the detection criteria. The aim of this study is to investigate the effect of magnetic twist in flaring coronal loops on the polarisation of gyro-synchrotron microwave (GSMW) emission, and determine whether it could provide a means for magnetic twist detection. We use time-dependent magnetohydrodynamic and test-particle models developed using LARE3D and GCA codes to investigate twisted coronal loops relaxing following the kink-instability. Synthetic GSMW emission maps (I and V Stokes components) are calculated using GX simulator. It is found that flaring twisted coronal loops produce GSMW radiation with a gradient of circular polarisation across the loop. However, these patterns may be visible only for a relatively short period of time due to fast magnetic reconfiguration after the instability. Their visibility also depends on the orientation and position of the loop on solar disk. Typically, it would be difficult to see these characteristic polarisation pattern in a twisted loop seen from the top (close to the centre of the solar disk), but easier in a twisted loop seen from the side (i.e. observed very close to the limb).

Observation of coronal EUV spectral lines sensitive to different temperatures offers an opportunity to evaluate the thermal structure and flows in flaring atmospheres. This can be used to estimate the partitioning between the thermal and kinetic energies. Our aim is to model large-scale (50-10000km) velocity distributions in order to interpret non-thermal broadening of different spectral EUV lines. The developed models allow us to understand the origin of the observed spectral line shifts and broadening, and link these features to particular physical effects in flaring atmospheres. We use ideal MHD to derive twisted magnetic fluxtube configurations in a stratified atmosphere. The evolution of these fluxtubes is followed using resistive MHD, with anomalous resistivity depending on the local density and temperature. The model also takes into account the thermal conduction and radiative losses. The model allows us to evaluate average velocities and velocity dispersions, which would be interpreted as `non-thermal' velocities in observations, at different temperatures for different parts of the models. Our models show qualitative and quantitative agreement with observations. The line-of-sight (LOS) velocity dispersions demonstrate substantial correlation with the temperature, increasing from about 20-30 km/s around 1 MK to about 200-400 km/s near 10-20 MK. The average LOS velocities also correlate with velocity dispersions, although they demonstrate a very strong scattering, compared to observations. We also note that near foot-points the velocity dispersions across the magnetic field are systematically lower than those along the field. We conclude, that the correlation between the flow velocities, velocity dispersions and temperatures are likely to indicate that the same heating mechanism is responsible for heating the plasma, its turbulisation and expansion/evaporation.

The goal of this study is to develop an approximate self-consistent description of particle motion in strongly magnetised solar corona. We derive a set of reduced drift-kinetic equations based on the assumption that the gyro-velocity distribution is Maxwellian. The equations are tested using simple 1D models.

Twisted magnetic fields should be ubiquitous in flare-producing active regions where the magnetic fields are strongly non-potential. It has been shown that reconnection in helical magnetic coronal loops results in plasma heating and particle acceleration distributed within a large volume, including the lower coronal and chromospheric sections of the loops. This scenario can be an alternative to the standard flare model, where particles are accelerated only in a small volume located in the upper corona. We use a combination of MHD simulations and test-particle methods, which describe the development of kink instability and magnetic reconnection in twisted coronal loops using resistive compressible MHD, and incorporate atmospheric stratification and large-scale loop curvature. The resulting distributions of hot plasma let us estimate thermal X-ray emission intensities. The electric and magnetic fields obtained are used to calculate electron trajectories using the guiding-centre approximation. These trajectories combined with the MHD plasma density distributions let us deduce synthetic HXR bremsstrahlung intensities. Our simulations emphasise that the geometry of the emission patterns produced by hot plasma in flaring twisted coronal loops can differ from the actual geometry of the underlying magnetic fields. The twist angles revealed by the emission threads (SXR) are consistently lower than the field-line twist present at the onset of the kink-instability. HXR emission due to the interaction of energetic electrons with the stratified background are concentrated at the loop foot-points in these simulations, even though the electrons are accelerated everywhere within the coronal volume of the loop. The maximum of HXR emission consistently precedes that of SXR emission, with the HXR light-curve being approximately proportional to the temporal derivative of the SXR light-curve.

For the first time, we demonstrate how an MHD avalanche might occur in a multi-threaded coronal loop. Considering 23 non-potential magnetic threads within a loop, we use 3D MHD simulations to show that only one thread needs to be unstable in order to start an avalanche even when the others are below marginal stability. This has significant implications for coronal heating in that it provides for energy dissipation with a trigger mechanism. The instability of the unstable thread follows the evolution determined in many earlier investigations. However, once one stable thread is disrupted, it coalesces with a neighbouring thread and this process disrupts other nearby threads. Coalescence with these disrupted threads then occurs leading to the disruption of yet more threads as the avalanche develops. Magnetic energy is released in discrete bursts as the surrounding stable threads are disrupted. The volume integrated heating, as a function of time, shows short spikes suggesting that the temporal form of the heating is more like that of \textit{nanoflares} than of constant heating.

The solar corona, the tenuous outer atmosphere of the Sun, is orders of magnitude hotter than the solar surface. This 'coronal heating problem' requires the identification of a heat source to balance losses due to thermal conduction, radiation and (in some locations) convection. The review papers in this Theo Murphy meeting issue present an overview of recent observational findings, large- and small-scale numerical modelling of physical processes occurring in the solar atmosphere and other aspects which may affect our understanding of the proposed heating mechanisms. At the same time, they also set out the directions and challenges which must be tackled by future research. In this brief introduction, we summarize some of the issues and themes which reoccur throughout this issue.

The fields of solar radiophysics and solar system radio physics, or radio heliophysics, will benefit immensely from an instrument with the capabilities projected for SKA. Potential applications include interplanetary scintillation (IPS), radio-burst tracking, and solar spectral radio imaging with a superior sensitivity. These will provide breakthrough new insights and results in topics of fundamental importance, such as the physics of impulsive energy releases, magnetohydrodynamic oscillations and turbulence, the dynamics of post-eruptive processes, energetic particle acceleration, the structure of the solar wind and the development and evolution of solar wind transients at distances up to and beyond the orbit of the Earth. The combination of the high spectral, time and spatial resolution and the unprecedented sensitivity of the SKA will radically advance our understanding of basic physical processes operating in solar and heliospheric plasmas and provide a solid foundation for the forecasting of space weather events.

Context. Heating the solar corona to several million degrees requires the conversion of magnetic energy into thermal energy. In this paper, we investigate whether an unstable magnetic thread within a coronal loop can destabilise a neighbouring magnetic thread. Aims. By running a series of simulations, we aim to understand under what conditions the destabilisation of a single magnetic thread can also trigger a release of energy in a nearby thread. Methods. The 3D magnetohydrodynamics code, Lare3d, is used to simulate the temporal evolution of coronal magnetic fields during a kink instability and the subsequent relaxation process. We assume that a coronal magnetic loop consists of non-potential magnetic threads that are initially in an equilibrium state. Results. The non-linear kink instability in one magnetic thread forms a helical current sheet and initiates magnetic reconnection. The current sheet fragments, and magnetic energy is released throughout that thread. We find that, under certain conditions, this event can destabilise a nearby thread, which is a necessary requirement for starting an avalanche of energy release in magnetic threads. Conclusions. It is possible to initiate an energy release in a nearby, non-potential magnetic thread, because the energy released from one unstable magnetic thread can trigger energy release in nearby threads, provided that the nearby structures are close to marginal stability.

In the present study we investigate magnetic reconnection in twisted magnetic fluxtubes with different initial configurations. In all considered cases, energy release is triggered by the ideal kink instability, which is itself the result of applying footpoint rotation to an initially potential field. The main goal of this work is to establish the influence of the field topology and various thermodynamic effects on the energy release process. Specifically, we investigate convergence of the magnetic field at the loop footpoints, atmospheric stratification, as well as thermal conduction. In all cases, the application of vortical driving at the footpoints of an initally potential field leads to an internal kink instability. With the exception of the curved loop with high footpoint convergence, the global geometry of the loop change little during the simulation. Footpoint convergence, curvature and atmospheric structure clearly influences the rapidity with which a loop achieves instability as well as the size of the subsequent energy release. Footpoint convergence has a stabilising influence and thus the loop requires more energy for instability, which means that the subsequent relaxation has a larger heating effect. Large-scale curvature has the opposite result: less energy is needed for instability and so the amount of energy released from the field is reduced. Introducing a stratified atmosphere gives rise to decaying wave phenomena during the driving phase, and also results in a loop that is less stable.

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.

We study a model of particle acceleration coupled with an MHD model of magnetic reconnection in unstable twisted coronal loops. The kink instability leads to the formation of helical currents with strong parallel electric fields resulting in electron acceleration. The motion of electrons in the electric and magnetic fields of the reconnecting loop is investigated using a test-particle approach taking into account collisional scattering. We discuss the effects of Coulomb collisions and magnetic convergence near loop footpoints on the spatial distribution and energy spectra of high-energy electron populations and possible implications on the hard X-ray emission in solar flares.

In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment. In this paper, 2D fluid simulations are presented of this merging process in order to understand the underlying physics, and better interpret the experimental data. These simulations examine the individual and combined effects of tight-aspect ratio geometry and two-fluid physics on the merging. The ideal self-driven flux-rope dynamics are coupled to the diffusion layer physics, resulting in a large range of phenomena. For resistive MHD simulations, the flux-ropes enter the sloshing regime for normalised resistivity eta < 1E-5. In Hall-MHD three regimes are found for the qualitative behaviour of the current sheet, depending on the ratio of the current sheet width to the ion-sound radius. These are a stable collisional regime, an open X-point regime, and an intermediate regime that is highly unstable to tearing-type instabilities. In toroidal axisymmetric geometry, the final state after merging is a MAST-like spherical tokamak with nested flux-surfaces. It is also shown that the evolution of simulated 1D radial density profiles closely resembles the Thomson scattering electron density measurements in MAST. An intuitive explanation for the origin of the measured density structures is proposed, based upon the results of the toroidal Hall-MHD simulations.

Context: Relaxation theory offers a straightforward method for estimating the energy that is released when a magnetic field becomes unstable, as a result of continual convective driving. Aims: We present new results obtained from nonlinear magnetohydrodynamic (MHD) simulations of idealised coronal loops. The purpose of this work is to determine whether or not the simulation results agree with Taylor relaxation, which will require a modified version of relaxation theory applicable to unbounded field configurations. Methods: A three-dimensional (3D) MHD Lagrangian-remap code is used to simulate the evolution of a line-tied cylindrical coronal loop model. This model comprises three concentric layers surrounded by a potential envelope; hence, being twisted locally, each loop configuration is distinguished by a piecewise-constant current profile. Initially, all configurations carry zero-net-current fields and are in ideally unstable equilibrium. The simulation results are compared with the predictions of helicity conserving relaxation theory. Results: For all simulations, the change in helicity is no more than 2% of the initial value; also, the numerical helicities match the analytically-determined values. Magnetic energy dissipation predominantly occurs via shock heating associated with magnetic reconnection in distributed current sheets. The energy release and final field profiles produced by the numerical simulations are in agreement with the predictions given by a new model of partial relaxation theory: the relaxed field is close to a linear force free state; however, the extent of the relaxation region is limited, while the loop undergoes some radial expansion. Conclusions: The results presented here support the use of partial relaxation theory, specifically, when calculating the heating-event distributions produced by ensembles of kink-unstable loops.

Context: The strong electric fields associated with magnetic reconnection in solar flares are a plausible mechanism to accelerate populations of high energy, non-thermal particles. One such reconnection scenario occurs at a 3D magnetic null point, where global plasma flows give rise to strong currents in the spine axis or fan plane. Aims: To understand the mechanism of charged particle energy gain in both the external drift region and the diffusion region associated with 3D magnetic reconnection. In doing so we evaluate the efficiency of resistive spine and fan models for particle acceleration, and find possible observables for each. Method: We use a full orbit test particle approach to study proton trajectories within electromagnetic fields that are exact solutions to the steady and incompressible magnetohydrodynamic equations. We study single particle trajectories and find energy spectra from many particle simulations. The scaling properties of the accelerated particles with respect to field and plasma parameters is investigated. Results: For fan reconnection, strong non-uniform electric drift streamlines can accelerate the bulk of the test particles. The highest energy gain is for particles that enter the current sheet, where an increasing "guide field" stabilises particles against ejection. The energy is only limited by the total electric potential energy difference across the fan current sheet. The spine model has both slow external electric drift speed and weak energy gain for particles reaching the current sheet. Conclusions: The electromagnetic fields of fan reconnection can accelerate protons to the high energies observed in solar flares, gaining up to 0.1 GeV for anomalous values of resistivity. However, the spine model, which gave a harder energy spectrum in the ideal case, is not an efficient accelerator after pressure constraints in the resistive model are included.

We review basic theoretical concepts in particle acceleration, with particular emphasis on processes likely to occur in regions of magnetic reconnection. Several new developments are discussed, including detailed studies of reconnection in three-dimensional magnetic field configurations (e.g., current sheets, collapsing traps, separatrix regions) and stochastic acceleration in a turbulent environment. Fluid, test-particle, and particle-in-cell approaches are used and results compared. While these studies show considerable promise in accounting for the various observational manifestations of solar flares, they are limited by a number of factors, mostly relating to available computational power. Not the least of these issues is the need to explicitly incorporate the electrodynamic feedback of the accelerated particles themselves on the environment in which they are accelerated. A brief prognosis for future advancement is offered.

It has been proposed that the million degree temperature of the corona is due to the combined effect of barely-detectable energy releases, so called nanoflares, that occur throughout the solar atmosphere. Alas, the nanoflare density and brightness implied by this hypothesis means that conclusive verification is beyond present observational abilities. Nevertheless, we investigate the plausibility of the nanoflare hypothesis by constructing a magnetohydrodynamic (MHD) model that can derive the energy of a nanoflare from the nature of an ideal kink instability. The set of energy-releasing instabilities is captured by an instability threshold for linear kink modes. Each point on the threshold is associated with a unique energy release and so we can predict a distribution of nanoflare energies. When the linear instability threshold is crossed, the instability enters a nonlinear phase as it is driven by current sheet reconnection. As the ensuing flare erupts and declines, the field transitions to a lower energy state, which is modelled by relaxation theory, i.e., helicity is conserved and the ratio of current to field becomes invariant within the loop. We apply the model so that all the loops within an ensemble achieve instability followed by energy-releasing relaxation. The result is a nanoflare energy distribution. Furthermore, we produce different distributions by varying the loop aspect ratio, the nature of the path to instability taken by each loop and also the level of radial expansion that may accompany loop relaxation. The heating rate obtained is just sufficient for coronal heating. In addition, we also show that kink instability cannot be associated with a critical magnetic twist value for every point along the instability threshold.

Context: It is thought likely that vast numbers of nanoflares are responsible for the corona having a temperature of millions of degrees. Current observational technologies lack the resolving power to confirm the nanoflare hypothesis. An alternative approach is to construct a magnetohydrodynamic coronal loop model that has the ability to predict nanoflare energy distributions. Aims: This paper presents the initial results generated by such a model. It predicts heating events with a range of sizes, depending on where the instability threshold for linear kink modes is encountered. The aims are to calculate the distribution of event energies and to investigate whether kink instability can be predicted from a single parameter. Methods: The loop is represented as a straight line-tied cylinder. The twisting caused by random photospheric motions is captured by two parameters, representing the ratio of current density to field strength for specific regions of the loop. Dissipation of the loop's magnetic energy begins during the nonlinear stage of the instability, which develops as a consequence of current sheet reconnection. After flaring, the loop evolves to the state of lowest energy where, in accordance with relaxation theory, the ratio of current to field is constant throughout the loop and helicity is conserved. Results: The results suggest that instability cannot be predicted by any simple twist-derived property reaching a critical value. The model is applied such that the loop undergoes repeated episodes of instability followed by energy-releasing relaxation. Hence, an energy distribution of the nanoflares produced is collated. Conclusions: The final energy distribution features two nanoflare populations that follow different power laws. The power law index for the higher energy population is more than sufficient for coronal heating.

High cadence, multiwavelength, optical observations of a solar active region are presented, obtained with the Swedish Solar Telescope. Two magnetic bright points are seen to separate in opposite directions at a constant velocity of 2.8km/s. After a separation distance of approximately 4400km is reached, multiple Ellerman bombs are observed in both H-alpha and Ca-K images. As a result of the Ellerman bombs, periodic velocity perturbations in the vicinity of the magnetic neutral line, derived from simultaneous MDI data, are generated with amplitude +/- 6km/s and wavelength 1000km. The velocity oscillations are followed by an impulsive brightening visible in H-alpha and Ca-K, with a peak intensity enhancement of 63%. We interpret these velocity perturbations as the magnetic field deformation necessary to trigger forced reconnection. A time delay of approximately 3min between the H-alpha wing and Ca-K observations indicate that the observed magnetic reconnection occurs at a height of 200km above the solar surface. These observations are consistent with theoretical predictions and provide the first observational evidence of microflare activity driven by forced magnetic reconnection.

Context. The primary energy release in solar flares is almost certainly due to magnetic reconnection, making this a strong candidate as a mechanism for particle acceleration. While particle acceleration in 2D geometries has been widely studied, investigations in 3D are a recent development. Two main classes of reconnection regimes at a 3D magnetic null point have been identified: fan and spine reconnection Aims. Here we investigate particle trajectories and acceleration during reconnection at a 3D null point, using a test particle numerical code, and compare the efficiency of the fan and spine regimes in generating an energetic particle population. Methods. We calculated the time evolution of the energy spectra. We discuss the geometry of particle escape from the two configurations and characterise the trapped and escaped populations. Results. We find that fan reconnection is less efficent than spine reconnection in providing seed particles to the region of strong electric field where acceleration is possible. The establishment of a steady-state spectrum requires approximately double the time in fan reconnection. The steady-state energy spectrum at intermediate energies (protons 1 keV to 0.1 MeV) is comparable in the fan and spine regimes. While in spine reconnection particle escape takes place in two symmetric jets along the spine, in fan reconnection no jets are produced and particles escape in the fan plane, in a ribbon-like structure.