Modelling the solar Transition Region with the use of an Adaptive Conduction (TRAC) method permits fast and accurate numerical solutions of the field-aligned hydrodynamic equations, capturing the enthalpy exchange between the corona and transition region, when the corona undergoes impulsive heating. The TRAC method eliminates the need for highly resolved numerical grids in the transition region and the commensurate very short time steps that are required for numerical stability. When employed with coarse spatial resolutions, typically achieved in multi-dimensional magnetohydrodynamic codes, the errors at peak density are less than 5% and the computation time is three orders of magnitude faster than fully resolved field-aligned models. This paper presents further examples that demonstrate the versatility and robustness of the method over a range of heating events, including impulsive and quasi-steady footpoint heating. A detailed analytical assessment of the TRAC method is also presented, showing that the approach works through all phases of an impulsive heating event because (i) the total radiative losses and (ii) the total heating when integrated over the transition region are both preserved at all temperatures under the broadening modifications of the method. The results from the numerical simulations complement this conclusion.

Context: This paper investigates the effectiveness of phase mixing as a coronal heating mechanism. A key quantity is the wave damping rate, $\gamma$, defined as the ratio of the heating rate to the wave energy. Aims: We investigate whether or not laminar phase-mixed Alfv\'en waves can have a large enough value of $\gamma$ to heat the corona. We also investigate the degree to which the $\gamma$ of standing Alfv\'en waves which have reached steady-state can be approximated with a relatively simple equation. Further foci of this study are the cause of the reduction of $\gamma$ in response to leakage of waves out of a loop, the quantity of this reduction, and how increasing the number of excited harmonics affects $\gamma$. Results: We find that at observed frequencies $\gamma$ is too small to heat the corona by approximately three orders of magnitude. Therefore, we believe that laminar phase mixing is not a viable stand-alone heating mechanism for coronal loops. We show that $\gamma$ is largest at resonance. We find our simple equation provides a good estimate for the damping rate (within approximately 10% accuracy) for resonant field lines. However, away from resonance, the equation provides a poor estimate, predicting $\gamma$ to be orders of magnitude too large. We find that leakage acts to reduce $\gamma$ but plays a negligible role if $\gamma$ is of the order required to heat the corona. If the wave energy follows a power spectrum with slope -5/3 then $\gamma$ grows logarithmically with the number of excited harmonics. If the number of excited harmonics is increased by much more than 100, then the heating is mainly caused by gradients that are parallel to the field rather than perpendicular to it. Therefore, in this case, the system is not heated mainly by phase mixing.

Context. Some of the most dynamic active regions are associated with complex photospheric magnetic configurations such as quadrupolar regions, and especially ones with a $\delta$-spot configuration and a strong Polarity Inversion Line (PIL). Aims. We study the formation and eruption of magnetic flux ropes in quadrupolar regions. Methods. 3D MHD simulations of the partial emergence of a highly twisted flux tube from the solar interior into a non-magnetized stratified atmosphere. Results. We emerge two $\Omega$-shaped loops forming a quadrupolar region. The emerging flux forms two initially separated bipoles, that later come in contact creating a $\delta$-spot central region. Above the two bipoles, two magnetic lobes expand and interact through a series of current sheets at the interface between them. Two recurrent confined eruptions are produced. In both cases, reconnection between sheared low-lying field lines forms a flux rope. Reconnection between the two lobes high in the atmosphere forms retracting field lines that push against the flux rope, creating a current sheet between them. It also forms a third magnetic lobe between the two emerged ones, that later acts as a strapping field. The flux rope eruptions are triggered when the reconnection between the flux ropes and the field above them becomes efficient enough to remove the tension of the overlying field. These reconnection events occur internally in the system. The first erupting flux rope almost fully reconnects with the overlying field. The second eruption is confined by the overlying field. During the confined eruption, the flux rope is enhanced in size, flux and twist, similar to confined-flare-to-flux-rope observations. Proxies of the emission reveal the two erupting filaments channels. A flare arcade is formed only in the second eruption due to the longer-lasting and more efficient reconnection below the flux rope.

Dependable predictions of the X-ray burst ashes and light curves require a stringent constraint on the composition of the accreted material as an input parameter. Lower metallicity models are generally based on a metal deficient donor and all metals are summed up in CNO abundances or solar metal distribution is assumed. In this work, we study the alteration of accreted composition due to spallation in the atmosphere of accreting neutron stars considering a cascading destruction process. We find that the inclusion of the cascading process brings the replenishment of CNO elements and overall survival probability is higher compared to isolated destruction of CNO elements. Spallation model provides the distribution of metals as a function of mass accretion rate. Multi-zone X-ray burst models calculated with reduced metallicities have enhanced abundances for high-mass nuclei in X-ray burst ashes. The increased metallicity due to the replenishment of CNO elements changes the composition of burst ashes compared to lower metallicity conditions. This will modify the thermal and compositional structure of accreted neutron star crusts.

We consider a series of MHD simulations in which a small amplitude, transverse velocity perturbation is introduced into a complex magnetic field. We analysed the deformation of the wave fronts as the perturbation propagates through the braided magnetic structures and explore the nature of Alfv\'enic wave phase mixing in this regime. Spatial gradients in the local Alfv\'en speed and variations in the length of magnetic field lines ensure that small scales form throughout the propagating wave front due to phase mixing. Additionally, the presence of complex, intricate current sheets associated with the background field locally modifies the polarisation of the wave front. The combination of these two effects enhances the rate of viscous dissipation, particularly in more complex field configurations. Unlike in classical phase mixing configurations, the greater spatial extent of Alfv\'en speed gradients ensures that wave energy is deposited over a larger cross-section of the magnetic structure. Further, the complexity of the background magnetic field ensures that small gradients in a wave driver can map to large gradients within the coronal plasma. The phase mixing of MHD waves in a complex magnetic field will progress throughout the braided volume. As a result, in a non-ideal regime wave energy will be dissipated over a greater cross-section than in classical phase mixing models. The formation rate of small spatial scales in a propagating wave front is a function of the complexity of the background magnetic field. As such, if the coronal field is sufficiently complex it remains plausible that phase mixing induced wave heating can contribute to maintaining observed temperatures. Furthermore, the weak compressibility of the transverse wave and the observed phase mixing pattern may provide seismological information about the nature of the background plasma.

Thermal non-equilibrium (TNE) is believed to be a potentially important process in understanding some properties of the magnetically closed solar corona. Through one-dimensional hydrodynamic models, this paper addresses the importance of the numerical spatial resolution, footpoint heating timescales and background heating on TNE. Inadequate transition region (TR) resolution can lead to significant discrepancies in TNE cycle behaviour, with TNE being suppressed in under-resolved loops. A convergence on the periodicity and plasma properties associated with TNE required spatial resolutions of less than 2 km for a loop of length 180 Mm. These numerical problems can be resolved using an approximate method that models the TR as a discontinuity using a jump condition, as proposed by Johnston et al. (2017a,b). The resolution requirements (and so computational cost) are greatly reduced while retaining good agreement with fully resolved results. Using this approximate method we (i) identify different regimes for the response of coronal loops to time-dependent footpoint heating including one where TNE does not arise and (ii) demonstrate that TNE in a loop with footpoint heating is suppressed unless the background heating is sufficiently small. The implications for the generality of TNE are discussed.

Aims: This paper presents 2.5D numerical experiments of Alfv\'en wave phase mixing and aims to assess the effects of nonlinearities on wave behaviour and dissipation. In addition, this paper aims to quantify how effective the model presented in this work is at providing energy to the coronal volume. Methods: The model is presented and explored through the use of several numerical experiments which were carried out using the Lare2D code. The experiments study footpoint driven Alfv\'en waves in the neighbourhood of a two-dimensional x-type null point with initially uniform density and plasma pressure. A continuous sinusoidal driver with a constant frequency is used. Each experiment uses different driver amplitudes to compare weakly nonlinear experiments with linear experiments. Results: We find that the wave trains phase-mix owing to variations in the length of each field line and variations in the field strength. The nonlinearities reduce the amount of energy entering the domain, as they reduce the effectiveness of the driver, but they have relatively little effect on the damping rate (for the range of amplitudes studied). The nonlinearities produce density structures which change the natural frequencies of the field lines and hence cause the resonant locations to move. The shifting of the resonant location causes the Poynting flux associated with the driver to decrease. Reducing the magnetic diffusivity increases the energy build-up on the resonant field lines, however, it has little effect on the total amount of energy entering the system. From an order of magnitude estimate, we show that the Poynting flux in our experiments is comparable to the energy requirements of the quiet Sun corona. However a (possibly unphysically) large amount of magnetic diffusion was used however and it remains unclear if the model is able to provide enough energy under actual coronal conditions.

We report on our three-dimensional (3D) magnetohydrodynamic (MHD) simulations of cylindrical weakly twisted flux tubes emerging from 18 Mm below the photosphere. We perform a parametric study, by varying the initial magnetic field strength ($B_0$), radius ($R$), twist ($\alpha)$ and length of the emerging part of the flux tube ($\lambda$) to investigate how these parameters affect the transfer of the magnetic field from the convection zone to the photosphere. We show that the efficiency of emergence at the photosphere (i.e. how strong the photospheric field will be in comparison to $B_0$) depends not only on the $B_0$ but also the morphology of the emerging field and the twist. We show that parameters such as $B_0$ and magnetic flux cannot alone determine whether a flux tube will emerge to the solar surface. For instance, high-$B_0$ (weak-$B_0$) fields may fail (succeed) to emerge at the photosphere, depending on their geometrical properties. We also show that the photospheric magnetic field strength can vary greatly for flux tubes with the same $B_0$ but different geometric properties. Moreover, in some cases we have found scaling laws, whereby the magnetic field strength scales with the local density as $B\propto \rho^\kappa$, where $\kappa \approx 1$ deeper in the convection zone and $\kappa <1$, close to the photosphere. The transition between the two values occurs approximately when the local pressure scale ($H_p$) becomes comparable to the diameter of the flux tube ($H_p\approx2R$). We derive forms to explain how and when these scaling laws appear and compare them with the numerical simulations.

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.

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.

We proposed that the use of an approximate "jump condition" at the solar transition region permits fast and accurate numerical solutions of the one dimensional hydrodynamic equations when the corona undergoes impulsive heating. In particular, it eliminates the need for the very short timesteps imposed by a highly resolved numerical grid. This paper presents further examples of the applicability of the method for cases of non-uniform heating, in particular, nanoflare trains (uniform in space but non-uniform in time) and spatially localised impulsive heating, including at the loop apex and base of the transition region. In all cases the overall behaviour of the coronal density and temperature shows good agreement with a fully resolved one dimensional model and is significantly better than the equivalent results from a 1D code run without using the jump condition but with the same coarse grid. A detailed assessment of the errors introduced by the jump condition is presented showing that the causes of discrepancy with the fully resolved code are (i) the neglect of the terms corresponding to the rate of change of total energy in the unresolved atmsophere, (ii) mass motions at the base of the transition region and (iii) for some cases with footpoint heating, an over-estimation of the radiative losses in the transition region.

We present a new computational approach that addresses the difficulty of obtaining the correct interaction between the solar corona and the transition region in response to rapid heating events. In the coupled corona, transition region and chromosphere system, an enhanced downward conductive flux results in an upflow (chromospheric evaporation). However, obtaining the correct upflow generally requires high spatial resolution in order to resolve the transition region. With an unresolved transition region, artificially low coronal densities are obtained because the downward heat flux jumps across the unresolved region to the chromosphere, underestimating the upflows. Here, we treat the lower transition region as a discontinuity that responds to changing coronal conditions through the imposition of a jump condition that is derived from an integrated form of energy conservation. To illustrate and benchmark this approach against a fully resolved one-dimensional model, we present field-aligned simulations of coronal loops in response to a range of impulsive (spatially uniform) heating events. We show that our approach leads to a significant improvement in the coronal density evolution than just when using coarse spatial resolutions insufficient to resolve the lower transition region. Our approach compensates for the jumping of the heat flux by imposing a velocity correction that ensures that the energy from the heat flux goes into driving the transition region dynamics, rather than being lost through radiation. Hence, it is possible to obtain improved coronal densities. The advantages of using this approach in both one-dimensional hydrodynamic and three-dimensional magnetohydrodynamic simulations are discussed.

We perform MHD modeling of a single bright coronal loop to include the interaction with a non-uniform magnetic field. The field is stressed by random footpoint rotation in the central region and its energy is dissipated into heating by growing currents through anomalous magnetic diffusivity that switches on in the corona above a current density threshold. We model an entire single magnetic flux tube, in the solar atmosphere extending from the high-beta chromosphere to the low-beta corona through the steep transition region. The magnetic field expands from the chromosphere to the corona. The maximum resolution is ~30 km. We obtain an overall evolution typical of loop models and realistic loop emission in the EUV and X-ray bands. The plasma confined in the flux tube is heated to active region temperatures (~3 MK) after ~2/3 hr. Upflows from the chromosphere up to ~100 km/s fill the core of the flux tube to densities above 10^9 cm^-3. More heating is released in the low corona than the high corona and is finely structured both in space and time.

Context. Observations of flux emergence indicate that rotational velocities may develop within sunspots. However, the dependence of this rotation on sub-photospheric field strength and twist remains largely unknown. Aims. We investigate the effects of varying the initial field strength and twist of an emerging sub-photospheric magnetic flux tube on the rotation of the sunspots at the photosphere. Methods. We consider a simple model of a stratified domain with a sub-photospheric interior layer and three overlying atmospheric layers. A twisted arched flux tube is inserted in the interior and is allowed to rise into the atmosphere. To achieve this, the MHD equations are solved using the Lagrangian-remap code, Lare3d. We perform a parameter study by independently varying the sub-photospheric magnetic field strength and twist. Results. Altering the initial field strength and twist significantly affects the tube's evolution and the rotational motions that develop at the photosphere. The rotation angle, vorticity, and current show a direct dependence on the initial field strength. We find that an increase in field strength increases the angle through which the fieldlines rotate, the length of fieldlines extending into the atmosphere, and the magnetic energy transported to the atmosphere. This also affects the amount of residual twist in the interior. The length of the fieldlines is crucial as we predict the twist per unit length equilibrates to a lower value on longer fieldlines. No such direct dependence is found when we modify the twist owing to the complex effect this has on the tension force acting on the tube. However, there is still a clear ordering in quantities such as the rotation angle, helicity, and free energy with higher initial twist cases being related to sunspots that rotate more rapidly, transporting more helicity and magnetic energy to the atmosphere.

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.

This short review paper gives an overview of recently observed transverse, propagating velocity perturbations in coronal loops. These ubiquitous perturbations are observed to undergo strong damping as they propagate. Using 3D numerical simulations of footpoint-driven transverse waves propagating in a coronal plasma with a cylindrical density structure, in combination with analytical modelling, it is demonstrated that the observed velocity perturbations can be understood in terms of coupling of different wave modes in the inhomogeneous boundaries of the loops. Mode coupling in the inhomogeneous boundary layers of the loops leads to the coupling of the transversal (kink) mode to the azimuthal (Alfven) mode, observed as the decay of the transverse kink oscillations. Both the numerical and analytical results show the spatial profile of the damped wave has a Gaussian shape to begin with, before switching to exponential decay at large heights. In addition, recent analysis of CoMP (Coronal Multi-channel Polarimeter) Doppler shift observations of large, off- limb, trans-equatorial loops shows that Fourier power at the apex appears to be higher in the high-frequency part of the spectrum than expected from theoretical models. This excess high-frequency FFT power could be tentative evidence for the onset of a cascade of the low-to-mid frequency waves into (Alfvenic) turbulence.

Context. Solar eruptions and high flare activity often accompany the rapid rotation of sunspots. The study of sunspot rotation and the mechanisms driving this motion are therefore key to our understanding of how the solar atmosphere attains the conditions necessary for large energy release. Aims. We aim to demonstrate and investigate the rotation of sunspots in a 3D numerical experiment of the emergence of a magnetic flux tube as it rises through the solar interior and emerges into the atmosphere. Furthermore, we seek to show that the sub-photospheric twist stored in the interior is injected into the solar atmosphere by means of a definitive rotation of the sunspots. Methods. A numerical experiment is performed to solve the 3D resistive magnetohydrodynamic (MHD) equations using a Lagrangian-Remap code. We track the emergence of a toroidal flux tube as it rises through the solar interior and emerges into the atmosphere investigating various quantities related to both the magnetic field and plasma. Results. Through detailed analysis of the numerical experiment, we find clear evidence that the photospheric footprints or sunspots of the flux tube undergo a rotation. Significant vertical vortical motions are found to develop within the two polarity sources after the field emerges. These rotational motions are found to leave the interior portion of the field untwisted and twist up the atmospheric portion of the field. This is shown by our analysis of the relative magnetic helicity as a significant portion of the interior helicity is transported to the atmosphere. In addition, there is a substantial transport of magnetic energy to the atmosphere. Rotation angles are also calculated by tracing selected fieldlines; the fieldlines threading through the sunspot are found to rotate through angles of up to 353 degrees over the course of the experiment.

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.

We report on numerical experiment of the recurrent onset of helical "blowout" jets in an emerging flux region (EFR). We find that these jets are running with velocities of $\sim 100-250\,\mathrm{kms^{-1}}$ and they transfer a vast amount of heavy plasma into the outer solar atmosphere. During their emission, they undergo an untwisting motion as a result of reconnection between the twisted emerging and the non-twisted pre-existing magnetic field in the solar atmosphere. For the first time in the context of blowout jets, we provide a direct evidence that their untwisting motion is associated with the propagation of torsional Alfv\'en waves in the corona.

We report on three-dimensional (3D) MHD simulations of recurrent mini CME-like eruptions in a small active region (AR), which is formed by the dynamical emergence of a twisted (not kink unstable) flux tube from the solar interior. The eruptions are developed as a result of repeated formation and expulsion of new flux ropes due to continous emergence and reconnection of sheared fieldlines along the polarity inversion line (PIL) of the AR. The acceleration of the eruptions is triggered by tether-cutting reconnection at the current sheet underneath the erupting field. We find that each explosive eruption is followed by reformation of a sigmoidal structure and a subsequent sigmoid-to-flare arcade transformation in the AR. These results might have implications for recurrent CMEs and eruptive sigmoids/flares observations and theoretical studies.

We provide a valid magnetohydrostatic equilibrium from the collapse of a 2D X-point in the presence of a finite plasma pressure, in which the current density is not simply concentrated in an infinitesimally thin, one-dimensional current sheet, as found in force-free solutions. In particular, we wish to determine if a finite pressure current sheet will still involve a singular current, and if so, what is the nature of the singularity. We use a full MHD code, with the resistivity set to zero, so that reconnection is not allowed, to run a series of experiments in which an X-point is perturbed and then is allowed to relax towards an equilibrium, via real, viscous damping forces. Changes to the magnitude of the perturbation and the initial plasma pressure are investigated systematically. The final state found in our experiments is a "quasi-static" equilibrium where the viscous relaxation has completely ended, but the peak current density at the null increases very slowly following an asymptotic regime towards an infinite time singularity. Using a high grid resolution allows us to resolve the current structures in this state both in width and length. In comparison with the well known pressureless studies, the system does not evolve towards a thin current sheet, but concentrates the current at the null and the separatrices. The growth rate of the singularity is found to be tD, with 0 < D < 1. This rate depends directly on the initial plasma pressure, and decreases as the pressure is increased. At the end of our study, we present an analytical description of the system in a quasi-static non-singular equilibrium at a given time, in which a finite thick current layer has formed at the null.

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.

We present a comprehensive review of MHD wave behaviour in the neighbourhood of coronal null points: locations where the magnetic field, and hence the local Alfven speed, is zero. The behaviour of all three MHD wave modes, i.e. the Alfven wave and the fast and slow magnetoacoustic waves, has been investigated in the neighbourhood of 2D, 2.5D and (to a certain extent) 3D magnetic null points, for a variety of assumptions, configurations and geometries. In general, it is found that the fast magnetoacoustic wave behaviour is dictated by the Alfven-speed profile. In a $\beta=0$ plasma, the fast wave is focused towards the null point by a refraction effect and all the wave energy, and thus current density, accumulates close to the null point. Thus, null points will be locations for preferential heating by fast waves. Independently, the Alfven wave is found to propagate along magnetic fieldlines and is confined to the fieldlines it is generated on. As the wave approaches the null point, it spreads out due to the diverging fieldlines. Eventually, the Alfven wave accumulates along the separatrices (in 2D) or along the spine or fan-plane (in 3D). Hence, Alfven wave energy will be preferentially dissipated at these locations. It is clear that the magnetic field plays a fundamental role in the propagation and properties of MHD waves in the neighbourhood of coronal null points. This topic is a fundamental plasma process and results so far have also lead to critical insights into reconnection, mode-coupling, quasi-periodic pulsations and phase-mixing.

Magnetic neutral points, where the magnitude of the magnetic field vanishes locally, are potential locations for energy conversion in the solar corona. The fact that the magnetic field is identically zero at these points suggests that for the study of current sheet formation and of any subsequent resistive dissipation phase, a finite beta plasma should be considered, rather than neglecting the plasma pressure as has often been the case in the past. The rapid dissipation of a finite current layer in non-force-free equilibrium is investigated numerically, after the sudden onset of an anomalous resistivity. The aim of this study is to determine how the energy is redistributed during the initial diffusion phase, and what is the nature of the outward transmission of information and energy. The resistivity rapidly diffuses the current at the null point. The presence of a plasma pressure allows the vast majority of the free energy to be transferred into internal energy. Most of the converted energy is used in direct heating of the surrounding plasma, and only about 3% is converted into kinetic energy, causing a perturbation in the magnetic field and the plasma which propagates away from the null at the local fast magnetoacoustic speed. The propagating pulses show a complex structure due to the highly non-uniform initial state. It is shown that this perturbation carries no net current as it propagates away from the null. The fact that, under the assumptions taken in this paper, most of the magnetic energy released in the reconnection converts internal energy of the plasma, may be highly important for the chromospheric and coronal heating problem.

Context. For the last thirty years, most of the studies on the relaxation of stressed magnetic fields in the solar environment have onlyconsidered the Lorentz force, neglecting plasma contributions, and therefore, limiting every equilibrium to that of a force-free field. Aims. Here we begin a study of the non-resistive evolution of finite beta plasmas and their relaxation to magnetohydrostatic states, where magnetic forces are balanced by plasma-pressure gradients, by using a simple 2D scenario involving a hydromagnetic disturbance to a uniform magnetic field. The final equilibrium state is predicted as a function of the initial disturbances, with aims to demonstrate what happens to the plasma during the relaxation process and to see what effects it has on the final equilibrium state. Methods. A set of numerical experiments are run using a full MHD code, with the relaxation driven by magnetoacoustic waves damped by viscous effects. The numerical results are compared with analytical calculations made within the linear regime, in which the whole process must remain adiabatic. Particular attention is paid to the thermodynamic behaviour of the plasma during the relaxation. Results. The analytical predictions for the final non force-free equilibrium depend only on the initial perturbations and the total pressure of the system. It is found that these predictions hold surprisingly well even for amplitudes of the perturbation far outside the linear regime. Conclusions. Including the effects of a finite plasma beta in relaxation experiments leads to significant differences from the force-free case.

This review provides an introduction to the generation and evolution of the Sun's magnetic field, summarising both observational evidence and theoretical models. The eleven year solar cycle, which is well known from a variety of observed quantities, strongly supports the idea of a large-scale solar dynamo. Current theoretical ideas on the location and mechanism of this dynamo are presented. The solar cycle influences the behaviour of the global coronal magnetic field and it is the eruptions of this field that can impact on the Earth's environment. These global coronal variations can be modelled to a surprising degree of accuracy. Recent high resolution observations of the Sun's magnetic field in quiet regions, away from sunspots, show that there is a continual evolution of a small-scale magnetic field, presumably produced by small-scale dynamo action in the solar interior. Sunspots, a natural consequence of the large-scale dynamo, emerge, evolve and disperse over a period of several days. Numerical simulations can help to determine the physical processes governing the emergence of sunspots. We discuss the interaction of these emerging fields with the pre-existing coronal field, resulting in a variety of dynamic phenomena.

Recent studies of NOAA active region 10953, by Okamoto {\it et al.} ({\it Astrophys. J. Lett.} {\bf 673}, 215, 2008; {\it Astrophys. J.} {\bf 697}, 913, 2009), have interpreted photospheric observations of changing widths of the polarities and reversal of the horizontal magnetic field component as signatures of the emergence of a twisted flux tube within the active region and along its internal polarity inversion line (PIL). A filament is observed along the PIL and the active region is assumed to have an arcade structure. To investigate this scenario, MacTaggart and Hood ({\it Astrophys. J. Lett.} {\bf 716}, 219, 2010) constructed a dynamic flux emergence model of a twisted cylinder emerging into an overlying arcade. The photospheric signatures observed by Okamoto {\it et al.} (2008, 2009) are present in the model although their underlying physical mechanisms differ. The model also produces two additional signatures that can be verified by the observations. The first is an increase in the unsigned magnetic flux in the photosphere at either side of the PIL. The second is the behaviour of characteristic photospheric flow profiles associated with twisted flux tube emergence. We look for these two signatures in AR 10953 and find negative results for the emergence of a twisted flux tube along the PIL. Instead, we interpret the photospheric behaviour along the PIL to be indicative of photospheric magnetic cancellation driven by flows from the dominant sunspot. Although we argue against flux emergence within this particular region, the work demonstrates the important relationship between theory and observations for the successful discovery and interpretation of signatures of flux emergence.

This paper reviews some of the many 3D numerical experiments of the emergence of magnetic fields from the solar interior and the subsequent interaction with the pre-existing coronal magnetic field. The models described here are idealized, in the sense that the internal energy equation only involves the adiabatic, Ohmic and viscous shock heating terms. However, provided the main aim is to investigate the dynamical evolution, this is adequate. Many interesting observational phenomena are explained by these models in a self-consistent manner.

We investigate the behaviour of nonlinear, nonideal Alfv\'en wave propagation within an inhomogeneous magnetic environment. The governing MHD equations are solved in 1D and 2D using both analytical techniques and numerical simulations. We find clear evidence for the ponderomotive effect and visco-resistive heating. The ponderomotive effect generates a longitudinal component to the transverse Alfv\'en wave, with a frequency twice that of the driving frequency. Analytical work shows the addition of resistive heating. This leads to a substantial increase in the local temperature and thus gas pressure of the plasma, resulting in material being pushed along the magnetic field. In 2D, our system exhibits phase mixing and we observe an evolution in the location of the maximum heating, i.e. we find a drifting of the heating layer. Considering Alfv\'en wave propagation in 2D with an inhomogeneous density gradient, we find that the equilibrium density profile is significantly modified by both the flow of density due to visco-resistive heating and the nonlinear response to the localised heating through phase mixing.

Our aim is to study the photospheric flux distribution of a twisted flux tube that emerges from the solar interior. We also report on the eruption of a new flux rope when the emerging tube rises into a pre-existing magnetic field in the corona. To study the evolution, we use 3D numerical simulations by solving the time-dependent and resistive MHD equations. We qualitatively compare our numerical results with MDI magnetograms of emerging flux at the solar surface. We find that the photospheric magnetic flux distribution consists of two regions of opposite polarities and elongated magnetic tails on the two sides of the polarity inversion line (PIL), depending on the azimuthal nature of the emerging field lines and the initial field strength of the rising tube. Their shape is progressively deformed due to plasma motions towards the PIL. Our results are in qualitative agreement with observational studies of magnetic flux emergence in active regions (ARs). Moreover, if the initial twist of the emerging tube is small, the photospheric magnetic field develops an undulating shape and does not possess tails. In all cases, we find that a new flux rope is formed above the original axis of the emerging tube that may erupt into the corona, depending on the strength of the ambient field.

In this paper we study the effects of a toroidal magnetic flux tube emerging into a magnetized corona, with an emphasis on large-scale eruptions. The orientation of the fields is such that the two flux systems are almost antiparallel when they meet. We follow the dynamic evolution of the system by solving the 3D MHD equations using a Lagrangian remap scheme. Multiple eruptions are found to occur. The physics of the trigger mechanisms are discussed and related to well-known eruption models.

In this paper we study the dynamics of toroidal flux tubes emerging from the solar interior, through the photosphere and into the corona. Many previous theoretical studies of flux emergence use a twisted cylindrical tube in the solar interior as the initial condition. Important insights can be gained from this model, however, it does have shortcomings. The axis of the tube never fully emerges as dense plasma becomes trapped in magnetic dips and restrains its ascent. Also, since the entire tube is buoyant, the main photospheric footpoints (sunspots) continually drift apart. These problems make it difficult to produce a convincing sunspot pair. We aim to address these problems by considering a different initial condition, namely a toroidal flux tube. We perform numerical experiments and solve the 3D MHD equations. The dynamics are investigated through a range of initial field strengths and twists. The experiments demonstrate that the emergence of toroidal flux tubes is highly dynamic and exhibits a rich variety of behaviour. In answer to the aims, however, if the initial field strength is strong enough, the axis of the tube can fully emerge. Also, the sunspot pair does not continually drift apart. Instead, its maximum separation is the diameter of the original toroidal tube.

Mode conversion occurs when a wave passes through a region where the sound and Alfven speeds are equal. At this point there is a resonance, which allows some of the incident wave to be converted into a different mode. We study this phenomenon in the vicinity of a two-dimensional, coronal null point. As a wave approaches the null it passes from low- to high-beta plasma, allowing conversion to take place. We simulate this numerically by sending in a slow magnetoacoustic wave from the upper boundary; as this passes through the conversion layer a fast wave can clearly be seen propagating ahead. Numerical simulations combined with an analytical WKB investigation allow us to determine and track both the incident and converted waves throughout the domain.

This paper extends the models of Craig & McClymont (1991) and McLaughlin & Hood (2004) to include finite $\beta$ and nonlinear effects. We investigate the nature of nonlinear fast magnetoacoustic waves about a 2D magnetic X-point. We solve the compressible and resistive MHD equations using a Lagrangian remap, shock capturing code (Arber et al. 2001) and consider an initial condition in $ {\bf{v}}\times{\bf{B}} \cdot {\hat{\bf{z}}}$ (a natural variable of the system). We observe the formation of both fast and slow oblique magnetic shocks. The nonlinear wave deforms the X-point into a 'cusp-like' point which in turn collapses to a current sheet. The system then evolves through a series of horizontal and vertical current sheets, with associated changes in connectivity, i.e. the system exhibits oscillatory reconnection. Our final state is non-potential (but in force balance) due to asymmetric heating from the shocks. Larger amplitudes in our initial condition correspond to larger values of the final current density left in the system. The inclusion of nonlinear terms introduces several new features to the system that were absent from the linear regime.

Estimating Alfven speeds is of interest in modelling the solar corona, studying the coronal heating problem and understanding the initiation and propagation of coronal mass ejections (CMEs). We assume here that the corona is in a magnetohydrostatic equilibrium and that, because of the low plasma beta, one may decouple the magnetic forces from pressure and gravity. The magnetic field is then described by a force-free field for which we perform a statistical study of the magnetic field strength with height for four different active regions. The plasma along each field line is assumed to be in a hydrostatic equilibrium. As a first approximation, the coronal plasma is assumed to be isothermal with a constant or varying gravity with height. We study a bipolar magnetic field with a ring distribution of currents, and apply this method to four active regions associated with different eruptive events. By studying the global properties of the magnetic field strength above active regions, we conclude that (i) most of the magnetic flux is localized within 50 Mm of the photosphere, (ii) most of the energy is stored in the corona below 150 Mm, (iii) most of the magnetic field strength decays with height for a nonlinear force-free field slower than for a potential field. The Alfven speed values in an isothermal atmosphere can vary by two orders of magnitude (up to 100000 km/s). The global properties of the Alfven speed are sensitive to the nature of the magnetic configuration. For an active region with highly twisted flux tubes, the Alfven speed is significantly increased at the typical height of the twisted flux bundles; in flaring regions, the average Alfven speeds are above 5000 km/s and depart strongly from potential field values.

Mode conversion in the region where the sound and Alfven speeds are equal is a complex process, which has been studied both analytically and numerically, and has been seen in observations. In order to further the understanding of this process we set up a simple, one-dimensional model, and examine wave propagation through this system using a combination of analytical and numerical techniques. Simulations are carried out in a gravitationally stratified atmosphere with a uniform, vertical magnetic field for both isothermal and non-isothermal cases. For the non-isothermal case a temperature profile is chosen to mimic the steep temperature gradient encountered at the transition region. In all simulations, a slow wave is driven on the upper boundary, thus propagating down from low-beta to high-beta plasma across the mode-conversion region. In addition, a detailed analytical study is carried out where we predict the amplitude and phase of the transmitted and converted components of the incident wave as it passes through the mode-conversion region. A comparison of these analytical predictions with the numerical results shows good agreement, giving us confidence in both techniques. This knowledge may be used to help determine wave types observed and give insight into which modes may be involved in coronal heating.

We have simulated the 3D emergence and interaction of two twisted flux tubes, which rise from the interior into the outer atmosphere of the Sun. We present evidence for the multiple formation and eruption of flux ropes inside the emerging flux systems and hot arcade-like structures in between them. Their formation is due to internal reconnection, occurring between oppositely directed, highly stretched and sheared fieldlines at photospheric heights. Most of the eruptions escape into the corona, but some are confined and fade away without leaving the low atmosphere. As these flux ropes erupt, new reconnected fieldlines accumulate around the main axis of the initial magnetic flux systems. We also show the complex 3D fieldline geometry and the structure of the multiple current sheets, which form as a result of the reconnection between the emerging flux systems.

At this time there does not exist a robust set of rules connecting low and high $\beta$ waves across the $\beta \approx 1$ layer. The work here contributes specifically to what happens when a low $\beta$ fast wave crosses the $\beta \approx 1$ layer and transforms into high $\beta$ fast and slow waves. The nature of fast and slow magnetoacoustic waves is investigated in a finite $\beta$ plasma in the neighbourhood of a two-dimensional null point. The linearised equations are solved in both polar and cartesian forms with a two-step Lax-Wendroff numerical scheme. Analytical work (e.g. small $\beta$ expansion and WKB approximation) also complement the work. It is found that when a finite gas pressure is included in magnetic equilibrium containing an X-type null point, a fast wave is attracted towards the null by a refraction effect and that a slow wave is generated as the wave crosses the $\beta \approx 1$ layer. Current accumulation occurs close to the null and along nearby separatrices. The fast wave can now \emph{pass through the origin} due to the non-zero sound speed, an effect not previously seen in related papers but clear seen for larger values of $\beta$. Some of the energy can now leave the region of the null point and there is again generation of a slow wave component (we find that the fraction of the incident wave converted to a slow wave is proportional to $\beta$). We conclude that there are two competing phenomena; the refraction effect (due to the variable Alfv\'en speed) and the contribution from the non-zero sound speed. These experiments illustrate the importance of the magnetic topology and of the location of the $\beta \approx 1$ layer in the system.

The nature of fast magnetoacoustic and Alfv\'en waves is investigated in a zero $\beta$ plasma. This gives an indication of wave propagation in the low $\beta$ solar corona. It is found that for a two-dimensional null point, the fast wave is attracted to that point and the front of the wave slows down as it approaches the null point, causing the current density to accumulate there and rise rapidly. Ohmic dissipation will extract the energy in the wave at this point. This illustrates that null points play an important role in the rapid dissipation of fast magnetoacoustic waves and suggests the location where wave heating will occur in the corona. The Alfv\'en wave behaves in a different manner in that the wave energy is dissipated along the separatrices. For Alfv\'en waves that are decoupled from fast waves, the value of the plasma $\beta$ is unimportant. However, the phenomenon of dissipating the majority of the wave energy at a specific place is a feature of both wave types.

This paper is a demonstration of how the WKB approximation can be used to help solve the linearised 3D MHD equations. Using Charpit's Method and a Runge-Kutta numerical scheme, we have demonstrated this technique for a potential 3D magnetic null point, ${\bf{B}}=(x,\epsilon y -(\epsilon +1)z)$. Under our cold plasma assumption, we have considered two types of wave propagation: fast magnetoacoustic and Alfv\'en waves. We find that the fast magnetoacoustic wave experiences refraction towards the magnetic null point, and that the effect of this refraction depends upon the Alfv\'en speed profile. The wave, and thus the wave energy, accumulates at the null point. We have found that current build up is exponential and the exponent is dependent upon $\epsilon$. Thus, for the fast wave there is preferential heating at the null point. For the Alfv\'en wave, we find that the wave propagates along the fieldlines. For an Alfv\'en wave generated along the fan-plane, the wave accumulates along the spine. For an Alfv\'en wave generated across the spine, the value of $\epsilon$ determines where the wave accumulation will occur: fan-plane ($\epsilon=1$), along the $x-$axis ($0<\epsilon <1$) or along the $y-$axis ($\epsilon>1$). We have shown analytically that currents build up exponentially, leading to preferential heating in these areas. The work described here highlights the importance of understanding the magnetic topology of the coronal magnetic field for the location of wave heating.

The nature of fast magnetoacoustic and Alfv\'en waves is investigated in a zero $\beta$ plasma in the neighbourhood of a pair of two-dimensional null points. This gives an indication of wave propagation in the low $\beta$ solar corona, for a more complicated magnetic configuration than that looked at by McLaughlin & Hood (2004). It is found that the fast wave is attracted to the null points and that the front of the wave slows down as it approaches the null point pair, with the wave splitting and part of the wave accumulating at one null and the rest at the other. Current density will then accumulate at these points and ohmic dissipation will then extract the energy in the wave at these points. This suggests locations where wave heating will occur in the corona. The Alfv\'en wave behaves in a different manner in that the wave accumulates along the separatrices. Hence, the current density will accumulate at this part of the topology and this is where wave heating will occur. However, the phenomenon of wave accumulation at a specific place is a feature of both wave types, and illustrates the importance of studying the topology of the corona when considering MHD wave propagation.

This paper is the third in a series of investigations by the authors. The nature of fast magnetoacoustic and Alfv\'en waves is investigated in a 2D $\beta=0$ plasma in the neighbourhood of two dipoles. We use both numerical simulations (two-step Lax-Wendroff scheme) and analytical techniques (WKB approximation). It is found that the propagation of the linear fast wave is dictated by the Alfv\'en speed profile and that close to the null, the wave is attracted to the neutral point. However, it is also found that in this magnetic configuration some of the wave can escape the refraction effect; this had not been seen in previous investigations by the authors. The wave split occurs near the regions of very high Alfv\'en speed (found near the loci of the two dipoles). Also, for the set-up investigated it was found that 40% of the wave energy accumulates at the null. Ohmic dissipation will then extract the wave energy at this point. The Alfv\'en wave behaves in a different manner in that part of the wave accumulates along the separatrices and part escapes. Hence, the current density will accumulate at this part of the topology and this is where wave heating will occur. The phenomenon of wave accumulation at a specific place is a feature of both wave types, as is the result that a fraction of the wave can now escape the numerical box when propagating in this magnetic configuration.

Recent numerical investigations of wave propagation near coronal magnetic null points (McLaughlin and Hood: Astron. Astrophys. 459, 641,2006) have indicated how a fast MHD wave partially converts into a slow MHD wave as the disturbance passes from a low-beta plasma to a high-beta plasma. This is a complex process and a clear understanding of the conversion mechanism requires the detailed investigation of a simpler model. An investigation of mode conversion in a stratified, isothermal atmosphere, with a uniform, vertical magnetic field is carried out, both numerically and analytically. In contrast to previous investigations of upward-propagating waves (Zhugzhda and Dzhalilov: Astron. Astrophys. 112, 16, 1982a; Cally: Astrophys. J. 548, 473, 2001), this paper studies the downward propagation of waves from a low-beta to high-beta environment. A simple expression for the amplitude of the transmitted wave is compared with the numerical solution.

The emergence of magnetic flux from the convection zone into the corona is an important process for the dynamical evolution of the coronal magnetic field. In this paper we extend our previous numerical investigations, by looking at the process of flux interaction as an initially twisted flux tube emerges into a plane parallel, coronal magnetic field. Significant differences are found in the dynamical appearance and evolution of the emergence process depending on the relative orientation between the rising flux system and any preexisting coronal field. When the flux systems are nearly anti-parallel, the experiments show substantial reconnection and demonstrate clear signatures of a high temperature plasma located in the high velocity outflow regions extending from the reconnection region. However, the cases that have a more parallel orientation of the flux systems show very limited reconnection and none of the associated features. Despite the very different amount of reconnection between the two flux systems, it is found that the emerging flux that is still connected to the original tube, reaches the same height as a function of time. As a compensation for the loss of tube flux, a clear difference is found in the extent of the emerging loop in the direction perpendicular to the main axis of the initial flux tube. Increasing amounts of magnetic reconnection decrease the volume, which confines the remaining tube flux.

We present the results of a set of three-dimensional numerical simulations of magnetic flux emergence from below the photosphere into the corona that include a uniform and horizontal coronal magnetic field mimicking a pre-existing large-scale coronal magnetic system. Cases with different relative orientations of the upcoming and coronal fields are studied. Upon contact, a concentrated current sheet with the shape of an arch or bridge is formed at the interface which marks the positions of maximum jump in the field vector between the two systems. Relative angles above 90 degrees yield abundant magnetic reconnection and plasma heating. The reconnection is seen to be intrisincally three-dimensional in nature, except at singular positions along the current sheet. It drives collimated high-speed and high-temperature outflows only a short distance from the reconnection site that propagate along the ambient magnetic field lines as jets. Due to the low plasma density in the corona, these jets may propagate over large distances and, therefore help distribute high-density and high-temperature plasma along these newly reconnected field lines. The experiments permit to discern and visualize the three-dimensional shape and relative position of the upcoming plasma hill, high-speed jets and subphotospheric flux system. As a result of the reconnection, magnetic field lines from the magnetized plasma below the surface end up as coronal field lines, thus causing a profound change in the connectivity of the magnetic regions in the corona. The experiments presented here thus yield a number of features repeatedly observed with the TRACE satellite and the YOHKOH-SXT detector, like the establishment of connectivity between emergent and pre-existing active regions, local heating and high-velocity outflows.