Evidence for the emergence of twisted flux tubes into the solar atmosphere has, so far, come from indirect signatures. In this work, we investigate the topological input of twisted flux tube emergence directly by studying helicity and winding fluxes. In magnetohydrodynamic simulations with domains spanning from the top of the convection zone to the lower corona, we simulate the emergence of twisted flux tubes with a range of different initial field strengths. One important feature of this work is the inclusion of a convectively-unstable layer beneath the photosphere. We find approximately self-similar behaviour in the helicity input for the different field strengths considered. As the tubes rise and reach the photosphere, there is a strong input of negative helicity since we consider left-handed twisted tubes. This phase is then followed by a reduction of the negative input and, for low initial field strengths, a net positive helicity input. This phase corresponds to the growing influence of convection on the field and the development of serpentine field structures during emergence. The winding flux can be used to detect when the twisted cores of the tubes reach the photosphere, giving clear information about the input of topologically complex magnetic field into the solar atmosphere. In short, the helicity and winding fluxes can provide much information about how a magnetic field emerges that is not directly available from other sources, such as magnetograms. In evaulating the helicity content of these simualtions we test numerous means for creating synthetic magnetograms, including methods which acount for both the evolving geometry and the finite extent of the photosphere. Whilst the general qualitative behaviours are same in each case, the different forms of averaging do affect the helicity and winding inputs quantitatively.

The kink instability of magnetohydrodynamics is believed to be fundamental to many aspects of the dynamic activity of the solar atmosphere, such as the initiation of flares and the heating of the solar corona. In this work, we investigate the importance of viscosity on the kink instability. In particular, we focus on two forms of viscosity; isotropic viscosity (independent of the magnetic field) and anisotropic viscosity (with a preferred direction following the magnetic field). Through the detailed analysis of magnetohydrodynamic simulations of the kink instability with both types of viscosity, we show that the form of viscosity has a significant effect on the nonlinear dynamics of the instability. The different viscosities allow for different flow and current structures to develop, thus affecting the behaviour of magnetic relaxation, the formation of secondary instabilities and the Ohmic and viscous heating produced. Our results have important consequences for the interpretation of solar observations of the kink instability.

Eruptive flares (EFs) are associated with erupting filaments and, in some models, filament eruption drives flare reconnection. Recently, however, observations of a confined flare (CF) have revealed all the hallmarks of an EF (impulsive phase, flare ribbons, etc.) without the filament eruption itself. Therefore, if the filament is not primarily responsible for impulsive flare reconnection, what is? In this Letter, we argue, based on mimimal requirements, that the plasmoid instability is a strong candidate for explaining the impulsive phase in the observed CF. We present magnetohydrodynamic simulation results of the nonlinear development of the plasmoid instability, in a model active region magnetic field geometry, to strengthen our claim. We also discuss how the ideas described in this Letter can be generalised to other situations, including EFs.

Magnetic helicity flux gives information about the topology of a magnetic field passing through a boundary. In solar physics applications, this boundary is the photosphere and magnetic helicity flux has become an important quantity in analysing magnetic fields emerging into the solar atmosphere. In this work we investigate the evolution of magnetic helicity flux in magnetohydrodynamic (MHD) simulations of solar flux emergence. We consider emerging magnetic fields with different topologies and investigate how the magnetic helicity flux patterns corresponds to the dynamics of emergence. To investigate how the helicity input is connected to the emergence process, we consider two forms of the helicity flux. The first is the standard form giving topological information weighted by magnetic flux. The second form represents the net winding and can be interpreted as the standard helicity flux less the magnetic flux. Both quantities provide important and distinct information about the structure of the emerging field and these quantities differ significantly for mixed sign helicity fields. A novel aspect of this study is that we account for the varying morphology of the photosphere due to the motion of the dense plasma lifted into the chromosphere. Our results will prove useful for the interpretation of magnetic helicity flux maps in solar observations.

We investigate the onset of the classical magnetohydrodynamic (MHD) tearing instability (TI) and focus on non-modal (transient) growth rather than the tearing mode. With the help of pseudospectral theory, the operators of the linear equations are shown to be highly non-normal, resulting in the possibility of significant transient growth at the onset of the TI. This possibility increases as the Lundquist number $S$ increases. In particular, we find evidence, numerically, that the maximum possible transient growth, measured in the $L_2$-norm, for the classical setup of current sheets unstable to the TI, scales as $O(S^{1/4})$ on time scales of $O(S^{1/4})$ for $S\gg 1$. This behaviour is much faster than the time scale $O(S^{1/2})$ when the solution behaviour is dominated by the tearing mode. The size of transient growth obtained is dependent on the form of the initial perturbation. Optimal initial conditions for the maximum possible transient growth are determined, which take the form of wave packets and can be thought of as noise concentrated at the current sheet. We also examine how the structure of the eigenvalue spectrum relates to physical quantities.

In this paper, we investigate the possibility of transient growth in the linear perturbation of current sheets. The resistive magnetohydrodynamic (MHD) operator for a background field consisting of a current sheet is non-normal, meaning that associated eigenvalues and eigenmodes can be very sensitive to perturbation. In a linear stability analysis of a tearing current sheet, we show that modes that are damped as $t\rightarrow\infty$ can produce transient energy growth, contributing faster growth rates and higher energy attainment (within a fixed finite time) than the unstable tearing mode found from normal-mode analysis. We determine the transient growth for tearing-stable and tearing-unstable regimes and discuss the consequences of our results for processes in the solar atmosphere, such as flares and coronal heating. Our results have significant potential impact on how fast current sheets can be disrupted. In particular, transient energy growth due to (asymptotically) damped modes may lead to accelerated current sheet thinning and, hence, a faster onset of the plasmoid instability, compared to the rate determined by the tearing mode alone.

Penumbrae are the manifestation of magnetoconvection in highly inclined (to the vertical direction) magnetic field. The penumbra of a sunspot tends to form, initially, along the arc of the umbra antipodal to the main region of flux emergence. The question of how highly inclined magnetic field can concentrate along the antipodal curves of umbrae, at least initially, remains to be answered. Previous observational studies have suggested the existence of some form of overlying magnetic canopy which acts as the progenitor for penumbrae. We propose that such overlying magnetic canopies are a consequence of how the magnetic field emerges into the atmosphere and are, therefore, part of the emerging region. We show, through simulations of twisted flux tube emergence, that canopies of highly inclined magnetic field form preferentially at the required locations above the photosphere.

We introduce to the stellar physics community a method of modelling stellar coronae that can be considered to be an extension of the potential field. In this approach, the magnetic field is coupled to the background atmosphere. The model is magnetohydrostatic (MHS) and is a balance between the Lorentz force, the pressure gradient and gravity. Analytical solutions are possible and we consider a particular class of equilibria in this paper. The model contains two free parameters and the effects of these on both the geometry and topology of the coronal magnetic field are investigated. A demonstration of the approach is given using a magnetogram derived from Zeeman-Doppler imaging of the 0.75 M$_{\odot}$ M-dwarf star GJ 182.

Aims: To investigate the relationship between surges and magnetic reconnection during the emergence of small-scale active regions. In particular, to examine how the large-scale geometry of the magnetic field, shaped by different phases of reconnection, guides the flowing of surges. Methods: We present three flux emergence models. The first model, and the simplest, consists of a region emerging into a horizontal ambient field that is initially parallel to the top of the emerging region. The second model is the same as the first but with an extra smaller emerging region which perturbs the main region. This is added to create a more complex magnetic topology and to test how this complicates the development of surges compared to the first model. The last model has a non-uniform ambient magnetic field to model the effects of emergence near a sunspot field and impose asymmetry on the system through the ambient magnetic field. At each stage, we trace the magnetic topology to identify the locations of reconnection. This allows for field lines to be plotted from different topological regions, highlighting how their geometry affects the development of surges. Results: In the first model, we identify distinct phases of reconnection. Each phase is associated with a particular geometry for the magnetic field and this determines the paths of the surges. The second model follows a similar pattern to the first but with a more complex magnetic topology and extra eruptions. The third model highlights how an asymmetric ambient field can result in preferred locations for reconnection, subsequently guiding the direction of surges.

We present an analysis of the formation of atmospheric flux ropes in a magnetohydrodynamic (MHD) solar flux emergence simulation. The simulation domain ranges from the top of the solar interior to the low corona. A twisted magnetic flux tube emerges from the solar interior and into the atmosphere where it interacts with the ambient magnetic field. By studying the connectivity of the evolving magnetic field, we are able to better understand the process of flux rope formation in the solar atmosphere. In the simulation, two flux ropes are produced as a result of flux emergence. Each has a different evolution resulting in different topological structures. These are determined by plasma flows and magnetic reconnection. As the flux rope is the basic structure of the coronal mass ejection (CME), we discuss the implications of our findings for solar eruptions.

Aims. Linear magnetohydrostatic (MHS) models of solar magnetic fields balance plasma pressure gradients, gravity and Lorentz forces where the current density is composed of a linear force-free component and a cross-field component that depends on gravitational stratification. In this paper, we investigate an efficient numerical procedure for calculating such equilibria. Methods. The MHS equations are reduced to two scalar elliptic equations - one on the lower boundary and the other within the interior of the computational domain. The normal component of the magnetic field is prescribed on the lower boundary and a multigrid method is applied on both this boundary and within the domain to find the poloidal scalar potential. Once solved to a desired accuracy, the magnetic field, plasma pressure and density are found using a finite difference method. Results. We investigate the effects of the cross-field currents on the linear MHS equilibria. Force-free and non-force-free examples are given to demonstrate the numerical scheme and an analysis of speed-up due to parallelization on a graphics processing unit (GPU) is presented. It is shown that speed-ups of x30 are readily achievable.

Oscillatory reconnection is a time-dependent magnetic reconnection mechanism that naturally produces periodic outputs from aperiodic drivers. This paper aims to quantify and measure the periodic nature of oscillatory reconnection for the first time. We solve the compressible, resistive, nonlinear MHD equations using 2.5D numerical simulations. We identify two distinct periodic regimes: the impulsive and stationary phases. In the impulsive phase, we find the greater the amplitude of the initial velocity driver, the longer the resultant current sheet and the earlier its formation. In the stationary phase, we find that the oscillations are exponentially decaying and for driving amplitudes 6.3 - 126.2 km/s, we measure stationary-phase periods in the range 56.3 - 78.9 s, i.e. these are high frequency (0.01 - 0.02 Hz) oscillations. In both phases, we find that the greater the amplitude of the initial velocity driver, the shorter the resultant period, but note that different physical processes and periods are associated with both phases. We conclude that the oscillatory reconnection mechanism behaves akin to a damped harmonic oscillator.

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.

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.