In order to efficiently analyse the vast amount of data generated by solar space missions and ground-base instruments, modern machine learning techniques such as decision trees, support vector machines (SVMs) and neural networks can be very useful. In this paper we present initial results from using a convolutional neural network (CNN) to analyse observations from the Atmospheric Imaging Assembly (AIA) in the 1600A wavelength. The data is pre-processed to locate flaring regions where flare ribbons are visible in the observations. The CNN is created and trained to automatically analyse the shape and position of the flare ribbons, by identifying whether each image belongs into one of four classes: two-ribbon flare, compact/circular ribbon flare, limb flare or quiet Sun, with the final class acting as a control for any data included in the training or test sets where flaring regions are not present. The network created can classify flare ribbon observations into any of the four classes with a final accuracy of 94%. Initial results show that most of the images are correctly classified with the compact flare class being the only class where accuracy drops below 90% ad some observations are wrongly classified as belonging to the limb class.

We discuss a family of Vlasov-Maxwell equilibrium distribution functions for current sheet equilibria that are intermediate cases between the Harris sheet and the force-free (or modified) Harris sheet. These equilibrium distribution functions have potential applications to space and astrophysical plasmas. The existence of these distribution function had been briefly discussed in by Harrison and Neukirch 2009, but here it is shown that their approach runs into problems in the limit where the guide field goes to zero. The nature of this problem will be discussed and an alternative approach will be suggested that avoids the problem. This is achieved by considering a slight variation of the magnetic field profile, which allows a smooth transition between the Harris and force-free Harris sheet cases.

Kinetic-scale current sheets observed in the solar wind are frequently approximately force-free despite the fact that their plasma $\beta$ is of the order of one. In-situ measurements have recently shown that plasma density and temperature often vary across the current sheets, while the plasma pressure is approximately uniform. In many cases these density and temperature variations are asymmetric with respect to the center of the current sheet. To model these observations theoretically we develop in this paper equilibria of kinetic-scale force-free current sheets that have plasma density and temperature gradients. The models can also be useful for analysis of stability and dissipation of the current sheets in the solar wind.

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

For the extrapolation of magnetic fields into the solar corona from measurements taken in the photosphere (or chromosphere) force-free magnetic fields are typically used. This does not take into account that the lower layers of the solar atmosphere are not force-free. While some numerical extrapolation methods using magnetohydrostatic magnetic fields have been suggested, a complementary and numerically comparatively cheap method is to use analytical magnetohydrostatic equilibria to extrapolate the magnetic field. In this paper, we present a new family of solutions for a special class of analytical three-dimensional magnetohydrostatic equilibria, which can be of use for such magnetic field extrapolation. The new solutions allow for the more flexible modelling of a transition from non-force-free to (linear) force-free magnetic fields. In particular, the height and width of the region where this transition takes place can be specified by choosing appropriate model parameters.

We present new analytical three-dimensional solutions of the magnetohydrostatic equations, which are applicable to the co-rotating frame of reference outside a rigidly rotating cylindrical body, and have potential applications to planetary magnetospheres and stellar coronae. We consider the case with centrifugal force only, and use a transformation method in which the governing equation for the "pseudo-potential" (from which the magnetic field can be calculated) becomes the Laplace partial differential equation. The new solutions extend the set of previously found solutions to those of a "fractional multipole" nature, and offer wider possibilities for modelling than before. We consider some special cases, and present example solutions.

We present a class of one-dimensional, strictly neutral, Vlasov-Maxwell equilibrium distribution functions for force-free current sheets, with magnetic fields defined in terms of Jacobian elliptic functions, extending the results of Abraham-Shrauner (Phys. Plasmas 20, 102117, 2013) to allow for non-uniform density and temperature profiles. To achieve this, we use an approach previously applied to the force-free Harris sheet by Kolotkov et al. (Phys. Plasmas 22, 112902, 2015). In one limit of the parameters, we recover the model of Kolotkov et al., while another limit gives a linear force-free field. We discuss conditions on the parameters such that the distribution functions are always positive, and give expressions for the pressure, density, temperature and bulk-flow velocities of the equilibrium, discussing differences from previous models. We also present some illustrative plots of the distribution function in velocity space.

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

Magneto-static models may overcome some of the issues facing force-free magnetic field extrapolations. So far they have seen limited use and have faced problems when applied to quiet-Sun data. Here we present a first application to an active region. We use solar vector magnetic field measurements gathered by the IMaX polarimeter during the flight of the \sunrise{} balloon-borne solar observatory in June 2013 as boundary condition for a magneto-static model of the higher solar atmosphere above an active region. The IMaX data are embedded in active region vector magnetograms observed with SDO/HMI. This work continues our magneto-static extrapolation approach, which has been applied earlier ({\it Paper I}) to a quiet Sun region observed with \sunrise{} I. In an active region the signal-to-noise-ratio in the measured Stokes parameters is considerably higher than in the quiet Sun and consequently the IMaX measurements of the horizontal photospheric magnetic field allow us to specify the free parameters of the model in a special class of linear magneto-static equilibria. The high spatial resolution of IMaX (110-130 km, pixel size 40 km) enables us to model the non-force-free layer between the photosphere and the mid chromosphere vertically by about 50 grid points. In our approach we can incorporate some aspects of the mixed beta layer of photosphere and chromosphere, e.g., taking a finite Lorentz force into account, which was not possible with lower resolution photospheric measurements in the past. The linear model does not, however, permit to model intrinsic nonlinear structures like strongly localized electric currents.

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

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

Results are presented of a first study of collisionless magnetic reconnection starting from a recently found exact nonlinear force-free Vlasov-Maxwell equilibrium. The initial state has a Harris sheet magnetic field profile in one direction and a non-uniform guide field in a second direction, resulting in a spatially constant magnetic field strength as well as a constant initial plasma density and plasma pressure. It is found that the reconnection process initially resembles guide field reconnection, but that a gradual transition to anti-parallel reconnection happens as the system evolves. The time evolution of a number of plasma parameters is investigated, and the results are compared with simulations starting from a Harris sheet equilibrium and a Harris sheet plus constant guide field equilibrium.

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.

Our aim is to model the 3D magnetic field structure of the upper solar atmosphere, including regions of non-negligible plasma beta. We use high-resolution photospheric magnetic field measurements from SUNRISE/IMaX as boundary condition for a magneto-static magnetic field model. The high resolution of IMaX allows us to resolve the interface region between photosphere and corona, but modelling this region is challenging for the following reasons. While the coronal magnetic field is thought to be force-free (the Lorentz-force vanishes), this is not the case in the mixed plasma $\beta$ environment in the photosphere and lower chromosphere. In our model, pressure gradients and gravity forces are taken self-consistently into account and compensate the non-vanishing Lorentz-force. Above a certain height (about 2 Mm) the non-magnetic forces become very weak and consequently the magnetic field becomes almost force-free. Here we apply a linear approach, where the electric current density consists of a superposition of a field-line parallel current and a current perpendicular to the Sun's gravity field. We illustrate the prospects and limitations of this approach and give an outlook for an extension towards a non-linear model.

The aim of this work is to investigate and characterise particle behaviour in a 3D magnetohydrodynamic (MHD) model of a reconnecting magnetic separator. We use a relativistic guiding-centre test-particle code to investigate electron and proton acceleration in snapshots from 3D MHD separator reconnection experiments, and compare the results with findings from an analytical separator reconnection model studied in a previous investigation. The behaviour and acceleration of large distributions of particles are examined in detail for both analytical and numerical separator reconnection models. Accelerated particle orbit trajectories are shown to follow the separator before leaving the system along the separatrix surface of one of the nulls (determined by particle species) in the system of both models. A sufficiently localised electric field about the separator causes the orbits to appear to follow the spine bounding the separatrix surface field lines instead. We analyse and discuss the locations and spread of accelerated particle orbit final positions, which are seen to change over time in the numerical separator reconnection model. We deduce a simple relationship between the final energy range of particle orbits and the model dimensions, and discuss its implications for observed magnetic separators in the solar corona.

While the exact acceleration mechanism of energetic particles during solar flares is (as yet) unknown, magnetic reconnection plays a key role both in the release of stored magnetic energy of the solar corona and the magnetic restructuring during a flare. Recent work has shown that special field lines, called separators, are common sites of reconnection in 3D numerical experiments. To date, 3D separator reconnection sites have received little attention as particle accelerators. We investigate the effectiveness of separator reconnection as a particle acceleration mechanism for electrons and protons. We study the particle acceleration using a relativistic guiding-centre particle code in a time-dependent kinematic model of magnetic reconnection at a separator. The effect upon particle behaviour of initial position, pitch angle and initial kinetic energy are examined in detail, both for specific (single) particle examples and for large distributions of initial conditions. The separator reconnection model contains several free parameters and we study the effect of changing these parameters upon particle acceleration, in particular in view of the final particle energy ranges which agree with observed energy spectra.

In solar flares, a large number of charged particles is accelerated to high energies. By which physical processes this is achieved is one of the main open problems in solar physics. It has been suggested that during a flare, regions of the rapidly relaxing magnetic field can form a collapsing magnetic trap (CMT) and that this trap may contribute to particle energisation.} In this Research Note we focus on a particular analytical CMT model based on kinematic magnetohydrodynamics. Previous investigations of particle acceleration for this CMT model focused on the non-relativistic energy regime. It is the specific aim of this Research Note to extend the previous work to relativistic particle energies. Particle orbits were calculated numerically using the relativistic guiding centre equations. We also calculated particle orbits using the non-relativistic guiding centre equations for comparison. For mildly relativistic energies the relativistic and non-relativistic particle orbits mainly agree well, but clear deviations are seen for higher energies. In particular, the final particle energies obtained from the relativistic calculations are systematically lower than the energies reached from the corresponding non-relativistic calculations, and the mirror points of the relativistic orbits are systematically higher than for the corresponding non-relativistic orbits. While the overall behaviour of particle orbits in CMTs does not differ qualitatively when using the relativistic guiding centre equations, there are a few systematic quantitative differences between relativistic and non-relativistic particle dynamics.

Collapsing magnetic traps (CMTs) have been suggested as one possible mechanism responsible for the acceleration of high-energy particles during solar flares. An important question regarding the CMT acceleration mechanism is which particle orbits escape and which are trapped during the time evolution of a CMT. While some models predict the escape of the majority of particle orbits, other more sophisticated CMT models show that, in particular, the highest-energy particles remain trapped at all times. The exact prediction is not straightforward because both the loss cone angle and the particle orbit pitch angle evolve in time in a CMT. Our aim is to gain a better understanding of the conditions leading to either particle orbit escape or trapping in CMTs.

Context. It has been suggested that collapsing magnetic traps may contribute to accelerating particles to high energies during solar flares. Aims. We present a detailed investigation of the energization processes of particles in collapsing magnetic traps, using a specific model. We also compare for the first time the energization processes in a symmetric and an asymmetric trap model. Methods. Particle orbits are calculated using guiding centre theory. We systematically investigate the dependence of the energization process on initial position, initial energy and initial pitch angle. Results. We find that in our symmetric trap model particles can gain up to about 50 times their initial energy, but that for most initial conditions the energy gain is more moderate. Particles with an initial position in the weak field region of the collapsing trap and with pitch angles around 90 degrees achieve the highest energy gain, with betatron acceleration of the perpendicular energy the dominant energization mechanism. For particles with smaller initial pitch angle, but still outside the loss cone, we find the possibility of a significant increase in parallel energy. This increase in parallel energy can be attributed to the curvature term in the parallel equation of motion and the associated energy gain happens in the center of the trap where the field line curvature has its maximum. We find qualitatively similar results for the asymmetric trap model, but with smaller energy gains and a larger number of particles escaping from the trap.

For axisymmetric models for coronal loops the relationship between the bifurcation points of magnetohydrodynamic (MHD) equilibrium sequences and the points of linear ideal MHD instability is investigated imposing line-tied boundary conditions. Using a well-studied example based on the Gold-Hoyle equilibrium, it is demonstrated that if the equilibrium sequence is calculated using the Grad-Shafranov equation, the instability corresponds to the second bifurcation point and not the first bifurcation point because the equilibrium boundary conditions allow for modes which are excluded from the linear ideal stability analysis. This is shown by calculating the bifurcating equilibrium branches and comparing the spatial structure of the solutions close to the bifurcation point with the spatial structure of the unstable mode. If the equilibrium sequence is calculated using Euler potentials the first bifurcation point of the Grad-Shafranov case is not found, and the first bifurcation point of the Euler potential description coincides with the ideal instability threshold. An explanation of this results in terms of linear bifurcation theory is given and the implications for the use of MHD equilibrium bifurcations to explain eruptive phenomena is briefly discussed.

Context: During solar flares a large number of charged particles are accelerated to high energies, but the exact mechanism responsible for this is, so far, still unclear. Acceleration in collapsing magnetic traps is one of the mechanisms proposed. Aims: In the present paper we want to extend previous 2D models for collapsing magnetic traps to 3D models and to 2D models with shear flow. Methods: We use analytic solutions of the kinematic magnetohydrodynamic (MHD) equations to construct the models. Particle orbits are calculated using the guiding centre approximation. Results: We present a general theoretical framework for constructing kinematic MHD models of collapsing magnetic traps in 3D and in 2D with shear flow. A few illustrative examples of collapsing trap models are presented, together with some preliminary studies of particle orbits. For these example orbits, the energy increases roughly by a factor of 5 or 6, which is consistent with the energy increase found in previous 2D models.

We present three-dimensional solutions of the magnetohydrostatic equations in the co-rotating frame of reference outside a magnetized rigidly rotating cylinder. We make no symmetry assumption for the magnetic field, but to be able to make analytical progress we neglect outflows and specify a particular form for the current density. The magnetohydrostatic equations can then be reduced to a single linear partial differential equation for a pseudo-potential $U$, from which the magnetic field can be calculated by differentiation. The equation for $U$ can be solved by standard methods. The solutions can also be used to determine the plasma pressure, density and temperature as functions of all three spatial coordinates. Despite the obvious limitations of this approach, it can for example be used as a simple tool to create three-dimensional models for the closed field line regions of rotating magnetospheres without rotational symmetry.

We present a method to include stereoscopic information about the three dimensional structure of flux tubes into the reconstruction of the coronal magnetic field. Due to the low plasma beta in the corona we can assume a force free magnetic field, with the current density parallel to the magnetic field lines. Here we use linear force free fields for simplicity. The method uses the line of sight magnetic field on the photosphere as observational input. The value of $\alpha$ is determined iteratively by comparing the reconstructed magnetic field with the observed structures. The final configuration is the optimal linear force solution constrained by both the photospheric magnetogram and the observed plasma structures. As an example we apply our method to SOHO MDI/EIT data of an active region. In the future it is planned to apply the method to analyse data from the SECCHI instrument aboard the STEREO mission.

Recent observations of the solar corona with the LASCO coronagraph on board of the SOHO spacecraft have revealed the occurrence of triple helmet streamers even during solar minimum, which occasionally go unstable and give rise to large coronal mass ejections. There are also indications that the slow solar wind is either a combination of a quasi-stationary flow and a highly fluctuating component or may even be caused completely by many small eruptions or instabilities. As a first step we recently presented an analytical method to calculate simple two-dimensional stationary models of triple helmet streamer configurations. In the present contribution we use the equations of time- dependent resistive magnetohydrodynamics to investigate the stability and the dynamical behaviour of these configurations. We particularly focus on the possible differences between the dynamics of single isolated streamers and triple streamers and on the way in which magnetic reconnection initiates both small scale and large scale dynamical behaviour of the streamers. Our results indicate that small eruptions at the helmet streamer cusp may incessantly accelerate small amounts of plasma without significant changes of the equilibrium configuration and might thus contribute to the non-stationary slow solar wind. On larger time and length scales, large coronal eruptions can occur as a consequence of large scale magnetic reconnection events inside the streamer configuration. Our results also show that triple streamers are usually more stable than a single streamer.

Knowledge of the structure of the coronal magnetic field is important for our understanding of many solar activity phenomena, e.g. flares and CMEs. However, the direct measurement of coronal magnetic fields is not possible with present methods, and therefore the coronal field has to be extrapolated from photospheric measurements. Due to the low plasma beta the coronal magnetic field can usually be assumed to be approximately force free, with electric currents flowing along the magnetic field lines. There are both observational and theoretical reasons which suggest that at least prior to an eruption the coronal magnetic field is in a nonlinear force free state. Unfortunately the computation of nonlinear force free fields is way more difficult than potential or linear force free fields and analytic solutions are not generally available. We discuss several methods which have been proposed to compute nonlinear force free fields and focus particularly on an optimization method which has been suggested recently. We compare the numerical performance of a newly developed numerical code based on the optimization method with the performance of another code based on an MHD relaxation method if both codes are applied to the reconstruction of a semi-analytic nonlinear force-free solution. The optimization method has also been tested for cases where we add random noise to the perfect boundary conditions of the analytic solution, in this way mimicking the more realistic case where the boundary conditions are given by vector magnetogram data. We find that the convergence properties of the optimization method are affected by adding noise to the boundary data and we discuss possibilities to overcome this difficulty.

Context: This paper presents a method which can be used to calculate models of the global solar corona from observational data. Aims: We present an optimization method for computing nonlinear magnetohydrostatic equilibria in spherical geometry with the aim to obtain self-consistent solutions for the coronal magnetic field, the coronal plasma density and plasma pressure using observational data as input. Methods: Our code for the self-consistent computation of the coronal magnetic fields and the coronal plasma solves the non-force-free magnetohydrostatic equilibria using an optimization method. Previous versions of the code have been used to compute non-linear force-free coronal magnetic fields from photospheric measurements in Cartesian and spherical geometry, and magnetostatic-equilibria in Cartesian geometry. We test our code with the help of a known analytic 3D equilibrium solution of the magnetohydrostatic equations. The detailed comparison between the numerical calculations and the exact equilibrium solutions is made by using magnetic field line plots, plots of density and pressure and some of the usual quantitative numerical comparison measures. Results: We find that the method reconstructs the equilibrium accurately, with residual forces of the order of the discretisation error of the analytic solution. The correlation with the reference solution is better than 99.9% and the magnetic energy is computed accurately with an error of <0.1%. Conclusions: We applied the method so far to an analytic test case. We are planning to use this method with real observational data as input as soon as possible.

CONTEXT: As the coronal magnetic field can usually not be measured directly, it has to be extrapolated from photospheric measurements into the corona. AIMS: We test the quality of a non-linear force-free coronal magnetic field extrapolation code with the help of a known analytical solution. METHODS: The non-linear force-free equations are numerically solved with the help of an optimization principle. The method minimizes an integral over the force-free and solenoidal condition. As boundary condition we use either the magnetic field components on all six sides of the computational box in Case I or only on the bottom boundary in Case II. We check the quality of the reconstruction by computing how well force-freeness and divergence-freeness are fulfilled and by comparing the numerical solution with the analytical solution. The comparison is done with magnetic field line plots and several quantitative measures, like the vector correlation, Cauchy Schwarz, normalized vector error, mean vector error and magnetic energy. RESULTS: For Case I the reconstructed magnetic field shows good agreement with the original magnetic field topology, whereas in Case II there are considerable deviations from the exact solution. This is corroborated by the quantitative measures, which are significantly better for Case I. CONCLUSIONS: Despite the strong nonlinearity of the considered force-free equilibrium, the optimization method of extrapolation is able to reconstruct it; however, the quality of reconstruction depends significantly on the consistency of the input data, which is given only if the known solution is provided also at the lateral and top boundaries, and on the presence or absence of flux concentrations near the boundaries of the magnetogram.

AIMS: We develop an optimization principle for computing stationary MHD equilibria. METHODS: Our code for the self-consistent computation of the coronal magnetic fields and the coronal plasma uses non-force-free MHD equilibria. Previous versions of the code have been used to compute non-linear force-free coronal magnetic fields from photospheric measurements. The program uses photospheric vector magnetograms and coronal EUV images as input. We tested our reconstruction code with the help of a semi-analytic MHD-equilibrium. The quality of the reconstruction was judged by comparing the exact and reconstructed solution qualitatively by magnetic field-line plots and EUV-images and quantitatively by several different numerical criteria. RESULTS: Our code is able to reconstruct the semi-analytic test equilibrium with high accuracy. The stationary MHD optimization code developed here has about the same accuracy as its predecessor, a non-linear force-free optimization code. The computing time for MHD-equilibria is, however, longer than for force-free magnetic fields. We also extended a well-known class of nonlinear force-free equilibria to the non-force-free regime for purposes of testing the code. CONCLUSIONS: We demonstrate that the code works in principle using tests with analytical equilibria, but it still needs to be applied to real data.

An analytical MHD model of a normal-polarity prominence with compressible flow is presented. The exact solution is constructed via a systematic nonlinear separation of variables method used to calculate several classes of MHD equilibria in Cartesian geometry and uniform gravity. Although the model is 2D, a third magnetic/velocity vector field component is included and the highly sheared fields observed in prominences are reproduced. A description is given of the balance of gas pressure gradient with gravity and the Lorentz or inertial forces acting along and across the prominence. It is found that the flow does not significantly influence the heating profile. The analyzed model has dimensions, plasma density, temperature and velocity profiles which agree with those in the observations literature.

The concept of hyperbolic flux tubes (HFTs) is a generalization of the concept of separator field lines for coronal magnetic fields with a trivial magnetic topology. An effective mechanism of a current layer formation in HFTs is proposed. This mechanism is called magnetic pinching and it is caused by large-scale shearing motions applied to the photospheric feet of HFTs in a way as if trying to twist the HFT. It is shown that in the middle of an HFT such motions produce a hyperbolic flow that causes an exponentially fast growth of the current density in a thin force-free current layer. The magnetic energy associated with the current layer that is built up over a few hours is sufficient for a large flare. Other implications of HFT pinching for solar flares are discussed as well.

A self-consistent mathematical description of the magnetic field of fluted sunspot penumbrae is presented. This description is based on an expansion of the nonlinear force-free magnetohydrostatic equations written in cylindrical coordinates. The lowest order solutions are mathematically equivalent to laminated force-free equilibria in Cartesian geometry. The lowest order solutions have no toroidal component of the magnetic field and the magnetic pressure does not vary with azimuth but the solutions allow arbitrary variations of the magnetic field components with azimuth. Explicit solutions are presented which have a realistic radial profile of the magnetic field strength and reproduce the basic features of the observations.

Many observations indicate the occurrence of ionised gas in the distant halos of galaxies (including our own). Since photoionisation by stars (mainly O stars, young or evolved low-mass stars depending on the kind of galaxy) does not seem to be exclusively responsible for the ionisation of the hydrogen filaments that should otherwise cool fast and recombine quickly, the question arises which extra energy source can produce the quasistationary ionisation. We show that stationary localized magnetic reconnection in current filaments may contribute to the ionisation of the extraplanar halo gas. In these filaments magnetic energy is dissipated. Consequently, the ionised as well as the neutral component are heated and re-ionised on a time scale significantly shorter than the recombination time scale. The amount of energy required for efficient re-ionisation can in principle easily be provided by the free magnetic energy. We present quasi-static models that are characterized by plasma temperatures and densities that agree well with the observed values for the diffuse ionised gas component of the interstellar medium. Plasma-neutral gas fluid simulations are made to show that the recombination induced dynamical reconnection process indeed works in a self-regulatory way.

Recent observations of the solar corona with the LASCO coronagraph on board of the SOHO spacecraft have revealed the occurrence of triple helmet streamers even during solar minimum, which occasionally go unstable and give rise to particularly huge coronal mass ejections. We present a method to calculate (semi-)analytically self-consistent stationary configurations of triple helmet streamers which can serve as input for stability considerations and dynamical calculations. The method is based on an asymptotic expansion procedure using the elongated structure of the streamers. The method is very flexible and can be used in both Cartesian and spherical geometry. We discuss the effects of magnetic shear, gravity and field-aligned flow on open field lines. Example solutions illustrating the influence of each of these features on the solution are presented.