Magnetic flux ropes (MFRs) are usually considered to be the magnetic structure that dominates the transport of helicity from the Sun into the heliosphere. They entrain a confined plasma within a helically organized magnetic structure and are able to cause geomagnetic activity. The formation, evolution and twist distribution of MFRs are issues subject to strong debate. Although different twist profiles have been suggested so far, none of them has been thoroughly explored yet. The aim of this paper is to present a theoretical study of the conditions under which MFRs with different twist profiles are kink stable and thereby shed some light on the aforementioned discussions. The magnetic field is modeled according to the circular-cylindrical analytical flux rope model in Nieves-Chinchilla et al. (2016) as well as the Lundquist and Gold-Hoyle models, and the kink stability is analyzed with a numerical method that has been developed based on Linton, Longcope, and Fisher (1996). The results are discussed in relation to MFR rotations, magnetic forces, the reversed chirality scenario and the expansion throughout the heliosphere, among others, providing a theoretical background to improve the current understanding of the internal magnetic configuration of CMEs. The data obtained by new missions like Parker Solar Probe or Solar Orbiter will give the opportunity to explore these results and ideas by observing MFRs closer than ever to the Sun.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

We study the propagation and dissipation of magnetohydrodynamic waves in a set of numerical models that each include a solar--like stratified atmosphere and a magnetic field with a null point. All simulations have the same magnetic field configuration but different transition region heights. Compressive wave packets introduced in the photospheric portion of the simulations refract towards the null and collapse it into a current sheet, which then undergoes reconnection. The collapsed null forms a current sheet due to a strong magnetic pressure gradient caused by the inability of magnetic perturbations to cross the null. Although the null current sheet undergoes multiple reconnection episodes due to repeated reflections off the lower boundary, we find no evidence of oscillatory reconnection arising from the dynamics of the null itself. Wave mode conversion around the null generates a series of slow mode shocks localized near each separatrix. The shock strength is asymmetric across each separatrix, and subsequent shock damping therefore creates a tangential discontinuity across each separatrix, with long--lived current densities. A parameter study of the injected wave energy to reach the null confirms our previous WKB estimates. Finally, using current estimates of the photospheric acoustic power, we estimate that the shock and Ohmic heating we describe may account for $\approx1-10\%$ of the radiative losses from coronal bright points with similar topologies, and are similarly insufficient to account for losses from larger structures such as ephemeral regions. At the same time, the dynamics are comparable to proposed mechanisms for generating type--II spicules.

It has been observationally well established that the magnetic configurations most favorable for producing energetic flaring events reside in delta-spots, a class of sunspots defined as having opposite polarity umbrae sharing a common penumbra. They are frequently characterized by extreme compactness, strong rotation and anti-Hale orientation. Numerous studies have shown that nearly all of the largest solar flares originate in delta-spots, making the understanding of these structures a fundamental step in predicting space weather. Despite their important influence on the space environment, surprisingly little is understood about the origin and behavior of delta-spots. In this paper, we perform a systematic study of the behavior of emerging flux ropes to test a theoretical model for the formation of delta-spots: the kink instability of emerging flux ropes. We simulated the emergence of highly twisted, kink-unstable flux ropes from the convection zone into the corona, and compared their photospheric properties to those of emerged weakly twisted, kink-stable flux ropes. We show that the photospheric manifestations of the emergence of highly twisted flux ropes closely match the observed properties of delta-spots, and we discuss the resulting implications for observations. Our results strongly support and extend previous theoretical work that suggested that the kink instability of emerging flux ropes is a promising candidate to explain delta-spot formation, as it reproduces their key characteristics very well.

We describe and analyze observations by the Solar Dynamics Observatory of the emergence of a small, bipolar active region within an area of unipolar magnetic flux that was surrounded by a circular, quiescent filament. Within only eight hours of the start of the emergence, a partial splitting of the filament and two consecutive coronal mass ejections took place. We argue that all three dynamic events occurred as a result of particular magnetic-reconnection episodes between the emerging bipole and the pre-existing coronal magnetic field. In order to substantiate our interpretation, we consider three-dimensional magnetohydrodynamic simulations that model the emergence of magnetic flux in the vicinity of a large-scale coronal flux rope. The simulations qualitatively reproduce most of the reconnection episodes suggested by the observations; as well as the filament-splitting, the first eruption, and the formation of sheared/twisted fields that may have played a role in the second eruption. Our results suggest that the position of emerging flux with respect to the background magnetic configuration is a crucial factor for the resulting evolution, while previous results suggest that parameters such as the orientation or the amount of emerging flux are important as well. This poses a challenge for predicting the onset of eruptions that are triggered by flux emergence, and it calls for a detailed survey of the relevant parameter space by means of numerical simulations.

We consider a self-consistent and physical approach to interacting dark energy models described by a Lagrangian, and identify a new class of models with variable dark energy sound speed. We show that if the interaction between dark energy in the form of quintessence and cold dark matter is purely momentum exchange this generally leads to a dark energy sound speed that deviates from unity. Choosing a specific sub-case, we study its phenomenology by investigating the effects of the interaction on the cosmic microwave background and linear matter power spectrum. We also perform a global fitting of cosmological parameters using CMB data, and compare our findings to $\Lambda$CDM.

This article completes and extends a recent study of the Grad-Shafranov (GS) reconstruction in toroidal geometry, as applied to a two and a half dimensional configurations in space plasmas with rotational symmetry. A further application to the benchmark study of an analytic solution to the toroidal GS equation with added noise shows deviations in the reconstructed geometry of the flux rope configuration, characterized by the orientation of the rotation axis, the major radius, and the impact parameter. On the other hand, the physical properties of the flux rope, including the axial field strength, and the toroidal and poloidal magnetic flux, agree between the numerical and exact GS solutions. We also present a real event study of a magnetic cloud flux rope from \textit{in situ} spacecraft measurements. The devised procedures for toroidal GS reconstruction are successfully executed. Various geometrical and physical parameters are obtained with associated uncertainty estimates. The overall configuration of the flux rope from the GS reconstruction is compared with the corresponding morphological reconstruction based on white-light images. The results show overall consistency, but also discrepancy in that the inclination angle of the flux rope central axis with respect to the ecliptic plane differs by about 20-30 degrees in the plane of the sky. We also compare the results with the original straight-cylinder GS reconstruction and discuss our findings.

We quantify the emergence and decay rates of preceder (p) and follower (f) sunspots within ten active regions from 2010-2014 using Space-weather Helioseismic Magnetic Imager Active Region Patch data. The sunspots are small- to mid-sized regions and contain a signed flux within a single polarity sunspot of $(1.1-6.5)\times 10^{21}$ Mx. The net unsigned flux within the regions, including plage, ranges from $(5.1-20)\times 10^{21}$ Mx. Rates are calculated with and without intensity contours to differentiate between sunspot formation and flux emergence. Signed flux emergence rates, calculated with intensity contours, for the p (f) spots average $6.8$ (4.9) $\times 10^{19}$ Mx h$^{-1}$, while decay rates are $-1.9 (-3.4)\times 10^{19}$ Mx h$^{-1}$. The mean, signed flux emergence rate of the regions, including plage, is $7.1 \times 10^{19}$ Mx h$^{-1}$ for a mean peak flux of $5.9 \times 10^{21}$ Mx. Using a synthesis of these results and others reported previously, there is a clear trend for larger flux regions to emerge faster than smaller ones. Observed emergence rates ($d{\phi}/dt$, Mx h$^{-1}$) scale with total signed peak flux, $\tilde{\phi}_{max}$, as a power law with an exponent of 0.36, i.e., $d{\phi}/dt = A \tilde{\phi}_{max}^{0.36}$. The observed rates may assist in constraining the boundary and initial conditions in simulations which already demonstrate increased rates for flux tubes with higher buoyancy and twist, or in the presence of a strong upflow. Overall, the observed emergence rates are smaller than those in simulations, which may indicate a slower rise of the flux in the interior than captured in simulations.

The discovery of clear criteria that can deterministically describe the eruptive state of a solar active region would lead to major improvements on space weather predictions. Using series of numerical simulations of the emergence of a magnetic flux rope in a magnetized coronal, leading either to eruptions or to stable configurations, we test several global scalar quantities for the ability to discriminate between the eruptive and the non-eruptive simulations. From the magnetic field generated by the three-dimensional magnetohydrodynamical simulations, we compute and analyse the evolution of the magnetic flux, of the magnetic energy and its decomposition into potential and free energies, and of the relative magnetic helicity and its decomposition. Unlike the magnetic flux and magnetic energies, magnetic helicities are able to markedly distinguish the eruptive from the non-eruptive simulations. We find that the ratio of the magnetic helicity of the current-carrying magnetic field to the total relative helicity presents the highest values for the eruptive simulations, in the pre-eruptive phase only. We observe that the eruptive simulations do not possess the highest value of total magnetic helicity. In the framework of our numerical study, the magnetic energies and the total relative helicity do not correspond to good eruptivity proxies. Our study highlights that the ratio of magnetic helicities diagnoses very clearly the eruptive potential of our parametric simulations. Our study shows that magnetic-helicity-based quantities may be very efficient for the prediction of solar eruptions.

Models for the evolution of the solar coronal magnetic field are vital for understanding solar activity, yet the best measurements of magnetic field lie at the photosphere, necessitating the development of coronal models which are \textit{"data-driven"} at the photosphere. We present an investigation to determine the feasibility and accuracy of such methods. Our validation framework uses a simulation of active region (AR) formation, modeling the emergence of magnetic flux from the convection zone to the corona, as a ground-truth dataset, to supply both the photospheric information, and to perform the validation of the data-driven method. We focus our investigation on how the data-driven model accuracy depends on the temporal frequency of the driving data. The Helioseismic and Magnetic Imager on NASA's Solar {\jl Dynamics} Observatory produces full-disk vector magnetic field measurements at a 12 minute cadence. Using our framework we show that ARs that emerge over 25 hours can be modeled by the data-driving method with only $\sim$1\% error in the free magnetic energy, assuming the photospheric information is specified every 12 minutes. However, for rapidly evolving features, under-sampling of the dynamics at this cadence leads to a strobe effect, generating large electric currents and incorrect coronal morphology and energies. We derive a sampling condition for the driving cadence based on the evolution of these small-scale features, and show that higher-cadence driving can lead to acceptable errors. Future work will investigate the source of errors associated with deriving plasma variables from the photospheric magnetograms as well as other sources of errors such as reduced resolution and instrument bias and noise.

We identify coronal mass ejections (CMEs) associated with magnetic clouds (MCs) observed near Earth by the Wind spacecraft from 2008 to mid-2012, a time period when the two STEREO spacecraft were well positioned to study Earth-directed CMEs. We find 31 out of 48 Wind MCs during this period can be clearly connected with a CME that is trackable in STEREO imagery all the way from the Sun to near 1 AU. For these events, we perform full 3-D reconstructions of the CME structure and kinematics, assuming a flux rope morphology for the CME shape, considering the full complement of STEREO and SOHO imaging constraints. We find that the flux rope orientations and sizes inferred from imaging are not well correlated with MC orientations and sizes inferred from the Wind data. However, velocities within the MC region are reproduced reasonably well by the image-based reconstruction. Our kinematic measurements are used to provide simple prescriptions for predicting CME arrival times at Earth, provided for a range of distances from the Sun where CME velocity measurements might be made. Finally, we discuss the differences in the morphology and kinematics of CME flux ropes associated with different surface phenomena (flares, filament eruptions, or no surface activity).

We perform nonlinear MHD simulations to study the propagation of magnetoacoustic waves from the photosphere to the low corona. We focus on a 2D system with a gravitationally stratified atmosphere and three photospheric concentrations of magnetic flux that give rise to a magnetic null point with a magnetic dome topology. We find that a single wavepacket introduced at the lower boundary splits into multiple secondary wavepackets. A portion of the packet refracts towards the null due to the varying Alfv\'en speed. Waves incident on the equipartition contour surrounding the null, where the sound and Alfv\'en speeds coincide, partially transmit, reflect, and mode convert between branches of the local dispersion relation. Approximately $15.5\%$ of the wavepacket's initial energy ($E_{input}$) converges on the null, mostly as a fast magnetoacoustic wave. Conversion is very efficient: $70\%$ of the energy incident on the null is converted to slow modes propagating away from the null, $7\%$ leaves as a fast wave, and the remaining $23\%$ (0.036$E_{input}$) is locally dissipated. The acoustic energy leaving the null is strongly concentrated along field lines near each of the null's four separatrices. The portion of the wavepacket that refracts towards the null, and the amount of current accumulation, depends on the vertical and horizontal wavenumbers and the centroid position of the wavepacket as it crosses the photosphere. Regions that refract towards or away from the null do not simply coincide with regions of open versus closed magnetic field or the local field orientation. We also modeled wavepacket propagation using a WKB method and found that it agrees qualitatively, though not quantitatively, with the results of the numerical simulation.

Transient collimated plasma eruptions in the solar corona, commonly known as coronal (or X-ray) jets, are among the most interesting manifestations of solar activity. It has been suggested that these events contribute to the mass and energy content of the corona and solar wind, but the extent of these contributions remains uncertain. We have recently modeled the formation and evolution of coronal jets using a three-dimensional (3D) magnetohydrodynamic (MHD) code with thermodynamics in a large spherical domain that includes the solar wind. Our model is coupled to 3D MHD flux-emergence simulations, i.e, we use boundary conditions provided by such simulations to drive a time-dependent coronal evolution. The model includes parametric coronal heating, radiative losses, and thermal conduction, which enables us to simulate the dynamics and plasma properties of coronal jets in a more realistic manner than done so far. Here we employ these simulations to calculate the amount of mass and energy transported by coronal jets into the outer corona and inner heliosphere. Based on observed jet-occurrence rates, we then estimate the total contribution of coronal jets to the mass and energy content of the solar wind to (0.4-3.0) % and (0.3-1.0) %, respectively. Our results are largely consistent with the few previous rough estimates obtained from observations, supporting the conjecture that coronal jets provide only a small amount of mass and energy to the solar wind. We emphasize, however, that more advanced observations and simulations are needed to substantiate this conjecture.

We review our understanding of ionized plasma and neutral gas coupling in the weakly ionized, stratified, electromagnetically-permeated regions of the Sun's chromosphere and Earth's ionosphere/thermosphere. Using representative models for each environment we derive fundamental descriptions of the coupling of the constituent parts to each other and to the electric and magnetic fields, and we examine the variation in magnetization of the ionized component. Using these descriptions we compare related phenomena in the two environments, and discuss electric currents, energy transfer and dissipation. We present a coupled theoretical and numerical study of plasma instabilities in the two environments that serves as an example of how the chromospheric and ionospheric communities can further collaborate. We also suggest future collaborative studies that will help improve our understanding of these two different atmospheres which share many similarities, but have large disparities in key quantities.

Views of two bright prominence eruptions trackable all the way to 1AU are here presented, using the heliospheric imagers on the Solar TErrestrial RElations Observatory (STEREO) spacecraft. The two events first erupted from the Sun on 2011 June 7 and 2012 August 31, respectively. Only these two examples of clear prominence eruptions observable this far from the Sun could be found in the STEREO image database, emphasizing the rarity of prominence eruptions this persistently bright. For the 2011 June event, a time-dependent 3-D reconstruction of the prominence structure is made using point-by-point triangulation. This is not possible for the August event due to a poor viewing geometry. Unlike the coronal mass ejection (CME) that accompanies it, the 2011 June prominence exhibits little deceleration from the Sun to 1 AU, as a consequence moving upwards within the CME. This demonstrates that prominences are not necessarily tied to the CME's magnetic structure far from the Sun. A mathematical framework is developed for describing the degree of self-similarity for the prominence's expansion away from the Sun. This analysis suggests only modest deviations from self-similar expansion, but close to the Sun the prominence expands radially somewhat more rapidly than self-similarity would predict.

Coronal jets are transient, collimated eruptions that occur in regions of predominantly open magnetic field in the solar corona. Our understanding of these events has greatly evolved in recent years but several open questions, such as the contribution of coronal jets to the solar wind, remain. Here we present an overview of the observations and numerical modeling of coronal jets, followed by a brief description of "next-generation" simulations that include an advanced description of the energy transfer in the corona ("thermodynamic MHD"), large spherical computational domains, and the solar wind. These new models will allow us to address some of the open questions.

Flux emergence is widely recognized to play an important role in the initiation of coronal mass ejections. The Chen-Shibata (2000) model, which addresses the connection between emerging flux and flux rope eruptions, can be implemented numerically to study how emerging flux through the photosphere can impact the eruption of a pre-existing coronal flux rope. The model's sensitivity to the initial conditions and reconnection micro-physics is investigated with a parameter study. In particular, we aim to understand the stability of the coronal flux rope in the context of X-point collapse and the effects of boundary driving in both unstratified and stratified atmospheres. In the absence of driving, we assess the behavior of waves in the vicinity of the X-point. With boundary driving applied, we study the effects of reconnection micro-physics and atmospheric stratification on the eruption. We find that the Chen-Shibata equilibrium can be unstable to an X-point collapse even in the absence of driving due to wave accumulation at the X-point. However, the equilibrium can be stabilized by reducing the compressibility of the plasma, which allows small-amplitude waves to pass through the X-point without accumulation. Simulations with the photospheric boundary driving evaluate the impact of reconnection micro-physics and atmospheric stratification on the resulting dynamics: we show the evolution of the system to be determined primarily by the structure of the global magnetic fields with little sensitivity to the micro-physics of magnetic reconnection; and in a stratified atmosphere, we identify a novel mechanism for producing quasi-periodic behavior at the reconnection site behind a rising flux rope as a possible explanation of similar phenomena observed in solar and stellar flares.

We present results from 3D magnetohydrodynamic (MHD) simulations of the emergence of a twisted convection zone flux tube into a pre-existing coronal dipole field. As in previous simulations, following the partial emergence of the sub-surface flux into the corona, a combination of vortical motions and internal magnetic reconnection forms a coronal flux rope. Then, in the simulations presented here, external reconnection between the emerging field and the pre-existing dipole coronal field allows further expansion of the coronal flux rope into the corona. After sufficient expansion, internal reconnection occurs beneath the coronal flux rope axis, and the flux rope erupts up to the top boundary of the simulation domain ($\sim$ 36 Mm above the surface). We find that the presence of a pre-existing field, orientated in a direction to facilitate reconnection with the emerging field, is vital to the fast rise of the coronal flux rope. The simulations shown in this paper are able to self-consistently create many of the surface and coronal signatures used by coronal mass ejection (CME) models. These signatures include: surface shearing and rotational motions; quadrupolar geometry above the surface; central sheared arcades reconnecting with oppositely orientated overlying dipole fields; the formation of coronal flux ropes underlying potential coronal field; and internal reconnection which resembles the classical flare reconnection scenario. This suggests that proposed mechanisms for the initiation of a CME, such as "magnetic breakout", are operating during the emergence of new active regions.

The helical kink instability of a twisted magnetic flux tube has been suggested as a trigger mechanism for solar filament eruptions and coronal mass ejections (CMEs). In order to investigate if estimations of the pre-eruptive twist can be obtained from observations of writhe in such events, we quantitatively analyze the conversion of twist into writhe in the course of the instability, using numerical simulations. We consider the line tied, cylindrically symmetric Gold-Hoyle flux rope model and measure the writhe using the formulae by Berger and Prior which express the quantity as a single integral in space. We find that the amount of twist converted into writhe does not simply scale with the initial flux rope twist, but depends mainly on the growth rates of the instability eigenmodes of higher longitudinal order than the basic mode. The saturation levels of the writhe, as well as the shapes of the kinked flux ropes, are very similar for considerable ranges of initial flux rope twists, which essentially precludes estimations of pre-eruptive twist from measurements of writhe. However, our simulations suggest an upper twist limit of $\sim 6\pi$ for the majority of filaments prior to their eruption.

There has been a long-lasting debate on the question of whether or not electric currents in solar active regions are neutralized. That is, whether or not the main (or direct) coronal currents connecting the active region polarities are surrounded by shielding (or return) currents of equal total value and opposite direction. Both theory and observations are not yet fully conclusive regarding this question, and numerical simulations have, surprisingly, barely been used to address it. Here we quantify the evolution of electric currents during the formation of a bipolar active region by considering a three-dimensional magnetohydrodynamic simulation of the emergence of a sub-photospheric, current-neutralized magnetic flux rope into the solar atmosphere. We find that a strong deviation from current neutralization develops simultaneously with the onset of significant flux emergence into the corona, accompanied by the development of substantial magnetic shear along the active region's polarity inversion line. After the region has formed and flux emergence has ceased, the strong magnetic fields in the region's center are connected solely by direct currents, and the total direct current is several times larger than the total return current. These results suggest that active regions, the main sources of coronal mass ejections and flares, are born with substantial net currents, in agreement with recent observations. Furthermore, they support eruption models that employ pre-eruption magnetic fields containing such currents.

We present a 2.5D MHD simulation of a magnetic flux rope (FR) propagating in the heliosphere and investigate the cause of the observed sharp plasma beta transition. Specifically, we consider a strong internal magnetic field and an explosive fast start, such that the plasma beta is significantly lower in the FR than the sheath region that is formed ahead. This leads to an unusual FR morphology in the first stage of propagation, while the more traditional view (e.g. from space weather simulations like Enlil) of a `pancake' shaped FR is observed as it approaches 1 AU. We investigate how an equipartition line, defined by a magnetic Weber number, surrounding a core region of a propagating FR can demarcate a boundary layer where there is a sharp transition in the plasma beta. The substructure affects the distribution of toroidal flux, with the majority of the flux remaining in a small core region which maintains a quasi-cylindrical structure. Quantitatively, we investigate a locus of points where the kinetic energy density of the relative inflow field is equal to the energy density of the transverse magnetic field (i.e. effective tension force). The simulation provides compelling evidence that at all heliocentric distances the distribution of toroidal magnetic flux away from the FR axis is not linear; with 80% of the toroidal flux occurring within 40% of the distance from the FR axis. Thus our simulation displays evidence that the competing ideas of a pancaking structure observed remotely can coexist with a quasi-cylindrical magnetic structure seen in situ.

We present results from 3D visco-resistive magnetohydrodynamic (MHD) simulations of the emergence of a convection zone magnetic flux tube into a solar atmosphere containing a pre-existing dipole coronal field, which is orientated to minimize reconnection with the emerging field. We observe that the emergence process is capable of producing a coronal flux rope by the transfer of twist from the convection zone as found in previous simulations. We find that this flux rope is stable, with no evidence of a fast rise, and that its ultimate height in the corona is determined by the strength of the pre-existing dipole field. We also find that although the electric currents in the initial convection zone flux tube are almost perfectly neutralized, the resultant coronal flux rope carries a significant net current. These results suggest that flux tube emergence is capable of creating non-current-neutralized stable flux ropes in the corona, tethered by overlying potential fields, a magnetic configuration that is believed to be the source of coronal mass ejections.

Magnetic reconnection in partially ionized plasmas is a ubiquitous phenomenon spanning the range from laboratory to intergalactic scales, yet it remains poorly understood and relatively little studied. Here, we present results from a self-consistent multi-fluid simulation of magnetic reconnection in a weakly ionized reacting plasma with a particular focus on the parameter regime of the solar chromosphere. The numerical model includes collisional transport, interaction and reactions between the species, and optically thin radiative losses. This model improves upon our previous work in Leake et al. 2012 \cite{Leake2012} by considering realistic chromospheric transport coefficients, and by solving a generalized Ohm's law that accounts for finite ion-inertia and electron-neutral drag. We find that during the two dimensional reconnection of a Harris current sheet with an initial width larger than the neutral-ion collisional coupling scale, the current sheet thins until its width becomes less than this coupling scale, and the neutral and ion fluids decouple upstream from the reconnection site. During this process of decoupling, we observe reconnection faster than the single-fluid Sweet-Parker prediction, with recombination and plasma outflow both playing a role in determining the reconnection rate. As the current sheet thins further and elongates it becomes unstable to the secondary tearing instability, and plasmoids are seen. The reconnection rate, outflows and plasmoids observed in this simulation provide evidence that magnetic reconnection in the chromosphere could be responsible for jet-like transient phenomena such as spicules and chromospheric jets.

We present results of 2.5D numerical simulations of the emergence of sub-surface magnetic flux into the solar atmosphere, with emerging flux regions ranging from $10^{18}$ to $10^{21}$ Mx, representing both ephemeral and active regions. We include the presence of neutral Hydrogen in the governing equations, improve upon previous models by including the ionization in the equation of state, and use a more realistic convection zone model. We find that ionization and recombination of plasma during the rise of a convection zone flux tube reduces the rise speed of the tube's axis. The presence of neutral Hydrogen allows the effective flow of mass across fieldlines, by the addition of a Pedersen resistivity to the generalized Ohm's law, which dissipates current perpendicular to the magnetic field. This causes an increase of up to 10% in the amount of magnetic in-plane flux supplied to the corona and a reduction of up to 89% in the amount of sub-surface plasma brought up into the corona. However, it also reduces the amount of free magnetic energy supplied to the corona, and thus does not positively affect the likelihood of creating unstable coronal structures.

We present results from the first self-consistent multi-fluid simulations of chromospheric magnetic reconnection in a weakly ionized reacting plasma. We simulate two dimensional magnetic reconnection in a Harris current sheet with a numerical model which includes ion-neutral scattering collisions, ionization, recombination, optically thin radiative loss, collisional heating, and thermal conduction. In the resulting tearing mode reconnection the neutral and ion fluids become decoupled upstream from the reconnection site, creating an excess of ions in the reconnection region and therefore an ionization imbalance. Ion recombination in the reconnection region, combined with Alfv\'{e}nic outflows, quickly removes ions from the reconnection site, leading to a fast reconnection rate independent of Lundquist number. In addition to allowing fast reconnection, we find that these non-equilibria partial ionization effects lead to the onset of the nonlinear secondary tearing instability at lower values of the Lundquist number than has been found in fully ionized plasmas.These simulations provide evidence that magnetic reconnection in the chromosphere could be responsible for jet-like transient phenomena such as spicules and chromospheric jets.

We present the science case for a broadband X-ray imager with high-resolution spectroscopy, including simulations of X-ray spectral diagnostics of both active regions and solar flares. This is part of a trilogy of white papers discussing science, instrument (Bandler et al. 2010), and missions (Bookbinder et al. 2010) to exploit major advances recently made in transition-edge sensor (TES) detector technology that enable resolution better than 2 eV in an array that can handle high count rates. Combined with a modest X-ray mirror, this instrument would combine arcsecondscale imaging with high-resolution spectra over a field of view sufficiently large for the study of active regions and flares, enabling a wide range of studies such as the detection of microheating in active regions, ion-resolved velocity flows, and the presence of non-thermal electrons in hot plasmas. It would also enable more direct comparisons between solar and stellar soft X-ray spectra, a waveband in which (unusually) we currently have much better stellar data than we do of the Sun.

Current coronal mass ejection (CME) models set their lower boundary to be in the lower corona. They do not calculate accurately the transfer of free magnetic energy from the convection zone to the magnetically dominated corona because they model the effects of flux emergence using kinematic boundary conditions or simply assume the appearance of flux at these heights. We test the importance of including dynamical flux emergence in CME modeling by simulating, in 2.5D, the emergence of sub-surface flux tubes into different coronal magnetic field configurations. We investigate how much free magnetic energy, in the form of shear magnetic field, is transported from the convection zone to the corona, and whether dynamical flux emergence can drive CMEs. We find that multiple coronal flux ropes can be formed during flux emergence, and although they carry some shear field into the corona, the majority of shear field is confined to the lower atmosphere. Less than 10% of the magnetic energy in the corona is in the shear field, and this, combined with the fact that the coronal flux ropes bring up significant dense material, means that they do not erupt. Our results have significant implications for all CME models which rely on the transfer of free magnetic energy from the lower atmosphere into the corona but which do not explicitly model this transfer. Such studies of flux emergence and CMEs are timely, as we have new capabilities to observe this with Hinode and SDO, and therefore to test the models against observations.

Numerical simulations of the helical ($m\!=\!1$) kink instability of an arched, line-tied flux rope demonstrate that the helical deformation enforces reconnection between the legs of the rope if modes with two helical turns are dominant as a result of high initial twist in the range $\Phi\gtrsim6\pi$. Such reconnection is complex, involving also the ambient field. In addition to breaking up the original rope, it can form a new, low-lying, less twisted flux rope. The new flux rope is pushed downward by the reconnection outflow, which typically forces it to break as well by reconnecting with the ambient field. The top part of the original rope, largely rooted in the sources of the ambient flux after the break-up, can fully erupt or be halted at low heights, producing a "failed eruption." The helical current sheet associated with the instability is squeezed between the approaching legs, temporarily forming a double current sheet. The leg-leg reconnection proceeds at a high rate, producing sufficiently strong electric fields that it would be able to accelerate particles. It may also form plasmoids, or plasmoid-like structures, which trap energetic particles and propagate out of the reconnection region up to the top of the erupting flux rope along the helical current sheet. The kinking of a highly twisted flux rope involving leg-leg reconnection can explain key features of an eruptive but partially occulted solar flare on 18 April 2001, which ejected a relatively compact hard X-ray and microwave source and was associated with a fast coronal mass ejection.

We study the evolution of the magnetic field in a Y-type current sheet subject to a brief, localized magnetic reconnection event. The reconnection produces up- and down-flowing reconnected flux tubes which rapidly decelerate when they hit the Y-lines and underlying magnetic arcade loops at the ends of the current sheet. This localized reconnection outflow followed by a rapid deceleration reproduces the observed behavior of post-CME downflowing coronal voids. These simulations support the hypothesis that these observed coronal downflows are the retraction of magnetic fields reconnected in localized patches in the high corona.

We present a novel model in which shortening of a magnetic flux tube following localized, three-dimensional reconnection generates strong gas-dynamic shocks around its apex. The shortening releases magnetic energy by progressing away from the reconnection site at the Alfven speed. This launches inward flows along the field lines whose collision creates a pair of gas-dynamic shocks. The shocks raise both the mass density and temperature inside the newly shortened flux tube. Reconnecting field lines whose initial directions differ by more that 100 degrees can produce a concentrated knot of plasma hotter that 20 MK, consistent with observations. In spite of these high temperatures, the shocks convert less than 10% of the liberated magnetic energy into heat - the rest remains as kinetic energy of bulk motion. These gas-dynamic shocks arise only when the reconnection is impulsive and localized in all three dimensions; they are distinct from the slow magnetosonic shocks of the Petschek steady-state reconnection model.

We show, theoretically and via MHD simulations, how a short burst of reconnection localized in three dimensions on a one-dimensional current sheet creates a pair of reconnected flux tubes. We focus on the post-reconnection evolution of these flux tubes, studying their velocities and shapes. We find that slow-mode shocks propagate along these reconnected flux tubes, releasing magnetic energy as in steady-state Petschek reconnection. The geometry of these three-dimensional shocks, however, differs dramatically from the classical two-dimensional geometry. They propagate along the flux tube legs in four isolated fronts, whereas in the two-dimensional Petschek model, they form a continuous, stationary pair of V-shaped fronts. We find that the cross sections of these reconnected flux tubes appear as teardrop shaped bundles of flux propagating away from the reconnection site. Based on this, we argue that the descending coronal voids seen by Yohkoh SXT, LASCO, and TRACE are reconnected flux tubes descending from a flare site in the high corona, for example after a coronal mass ejection. In this model, these flux tubes would then settle into equilibrium in the low corona, forming an arcade of post-flare coronal loops.

The discrete nature of the solar magnetic field as it emerges into the corona through the photosphere indicates that it exists as isolated flux tubes in the convection zone, and will remain as discrete flux tubes in the corona until it collides and reconnects with other coronal fields. Collisions of these flux tubes will in general be three dimensional, and will often lead to reconnection, both rearranging the magnetic field topology in fundamental ways, and releasing magnetic energy. With the goal of better understanding these dynamics, we carry out a set of numerical experiments exploring fundamental characteristics of three dimensional magnetic flux tube reconnection. We first show that reconnecting flux tubes at opposite extremes of twist behave very differently: in some configurations, low twist tubes slingshot while high twist tubes tunnel. We then discuss a theory explaining these differences: by assuming helicity conservation during the reconnection one can show that at high twist, tunneled tubes reach a lower magnetic energy state than slingshot tubes, whereas at low twist the opposite holds. We test three predictions made by this theory. 1) We find that the level of twist at which the transition from slingshot to tunnel occurs is about two to three times higher than predicted on the basis of energetics and helicity conservation alone, probably because the dynamics of the reconnection play a large role as well. 2) We find that the tunnel occurs at all flux tube collision angles predicted by the theory. 3) We find that the amount of magnetic energy a slingshot or a tunnel reconnection releases agrees reasonably well with the theory, though at the high resistivities we have to use for numerical stability, a significant amount of magnetic energy is lost to diffusion, independent of reconnection.