Magnetic reconnection, a fundamentally important process in many aspects of astrophysics, is believed to be initiated by the tearing instability of an electric current sheet, a region where magnetic field abruptly changes direction and electric currents build up. Recent studies have suggested that the amount of magnetic shear in these structures is a critical parameter for the switch-on nature of magnetic reconnection in the solar atmosphere, at fluid spatial scales much larger than kinetic scales. We present results of simulations of reconnection in 3D current sheets with conditions appropriate to the solar corona. Using high-fidelity simulations, we follow the evolution of the linear and non-linear 3D tearing instability, leading to reconnection. We find that, depending on the parameter space, magnetic shear can play a vital role in the onset of significant energy release and heating via non-linear tearing. Two regimes in our study exist, dependent on whether the current sheet is longer or shorter than the wavelength of the fastest growing parallel mode (in the corresponding infinite system), thus determining whether sub-harmonics are present in the actual system. In one regime, where the fastest growing parallel mode has sub-harmonics, the non-linear interaction of these sub-harmonics and the coalescence of 3D plasmoids dominates the non-linear evolution, with magnetic shear playing only a weak role in the amount of energy released. In the second regime, where the fastest growing parallel mode has no-sub-harmonics, then only strongly sheared current sheets, where oblique mode are strong enough to compete with the dominant parallel mode, show any significant energy release. We expect both regimes to exist on the Sun, and so our results have important consequences for the the question of reconnection onset in different solar physics applications.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

Relative magnetic helicity is a gauge invariant quantity suitable for the study of the magnetic helicity content of heliospheric plasmas. Relative magnetic helicity can be decomposed uniquely into two gauge invariant quantities, the magnetic helicity of the non-potential component of the field, and a complementary volume-threading helicity. Recent analysis of numerical experiments simulating the generation of solar eruptions have shown that the ratio of the non-potential helicity to the total relative helicity is a clear marker of the eruptivity of the magnetic system, and that the high value of that quantity could be a sufficient condition for the onset of the instability generating the eruptions. The present study introduces the first analytical examination of the time variations of these non-potential and volume-threading helicities. The validity of the analytical formulas derived are confirmed with analysis of three-dimensional (3D) magnetohydrodynamics (MHD) simulations of solar coronal dynamics. Both the analytical investigation, and the numerical application show that, unlike magnetic helicity, the non-potential and the volume-threading helicities are not conserved quantities, even in the ideal MHD regime. A term corresponding to the transformation between the non-potential and volume-threading helicities frequently dominates their dynamics. This finding has an important consequence for their estimation in the solar corona: unlike with relative helicity, their volume coronal evolution cannot be ascertained by the flux of these quantities through the volume's boundaries. Only techniques extrapolating the 3D coronal field will enable both the proper study of the non-potential and volume-threading helicities, and the observational analysis of helicity-based solar-eruptivity proxies.

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.

It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years.

Alfv\'enic waves have been proposed as an important energy transport mechanism in coronal loops, capable of delivering energy to both the corona and chromosphere and giving rise to many observed features, of flaring and quiescent regions. In previous work, we established that resistive dissipation of waves (ambipolar diffusion) can drive strong chromospheric heating and evaporation, capable of producing flaring signatures. However, that model was based on a simplified assumption that the waves propagate instantly to the chromosphere, an assumption which the current work removes. Via a ray tracing method, we have implemented traveling waves in a field-aligned hydrodynamic simulation that dissipate locally as they propagate along the field line. We compare this method to and validate against the magnetohydrodynamics code Lare3D. We then examine the importance of travel times to the dynamics of the loop evolution, finding that (1) the ionization level of the plasma plays a critical role in determining the location and rate at which waves dissipate; (2) long duration waves effectively bore a hole into the chromosphere, allowing subsequent waves to penetrate deeper than previously expected, unlike an electron beam whose energy deposition rises in height as evaporation reduces the mean-free paths of the electrons; (3) the dissipation of these waves drives a pressure front that propagates to deeper depths, unlike energy deposition by an electron beam.

The physical conditions that determine whether or not solar active regions (ARs) produce strong flares and coronal mass ejections (CMEs) are not yet well understood. Here we investigate the association between electric-current neutralization, magnetic shear along polarity inversion lines (PILs), and eruptive activity in four ARs; two emerging and two well-developed ones. We find that the CME-producing ARs are characterized by a strongly non-neutralized total current, while the total current in the ARs that did not produce CMEs is almost perfectly neutralized. The difference in the PIL-shear between these two groups is much less pronounced, which suggests that the degree of current-neutralization may serve as a better proxy for assessing the ability of ARs to produce CMEs.

Solar flares and coronal mass ejections (CMEs) are among the most energetic events in the solar system, impacting the near-Earth environment. Flare productivity is empirically known to be correlated with the size and complexity of active regions. Several indicators, based on magnetic-field data from active regions, have been tested for flare forecasting in recent years. None of these indicators, or combinations thereof, have yet demonstrated an unambiguous eruption or flare criterion. Furthermore, numerical simulations have been only barely used to test the predictability of these parameters. In this context, we used the 3D parametric MHD numerical simulations of the self-consistent formation of the flux emergence of a twisted flux tube, inducing the formation of stable and unstable magnetic flux ropes of Leake (2013, 2014). We use these numerical simulations to investigate the eruptive signatures observable in various magnetic scalar parameters and provide highlights on data analysis processing. Time series of 2D photospheric-like magnetograms are used from parametric simulations of stable and unstable flux emergence, to compute a list of about 100 different indicators. This list includes parameters previously used for operational forecasting, physical parameters used for the first time, as well as new quantities specifically developed for this purpose. Our results indicate that only parameters measuring the total non-potentiality of active regions associated with magnetic inversion line properties, such as the Falconer parameters $L_{ss}$, $WL_{ss}$, $L_{sg}$ and $WL_{sg}$, as well as the new current integral $WL_{sc}$ and length $L_{sc}$ parameters, present a significant ability to distinguish the eruptive cases of the model from the non-eruptive cases, possibly indicating that they are promising flare and eruption predictors.

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.

Solar flares are driven by the release of magnetic energy from reconnection events in the solar corona, whereafter energy is transported to the chromosphere, heating the plasma and causing the characteristic radiative losses. In the collisional thick-target model, electrons accelerated to energies exceeding 10 keV traverse the corona and impact the chromosphere, where they deposit their energy through collisions with the much denser plasma in the lower atmosphere. While there are undoubtedly high energy non-thermal electrons accelerated in flares, it is unclear whether these electron beams are the sole mechanism of energy transport, or whether they only dominate in certain phases of the flare's evolution. Alfv\'enic waves are generated during the post-reconnection relaxation of magnetic field lines, so it is important to examine their role in energy transport.

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.

Sunquakes are the surface signatures of acoustic waves in the Sun's interior that are produced by some but not all flares and coronal mass ejections (CMEs). This paper explores a mechanism for sunquake generation by the magnetic field changes that occur during flares and CMEs, using MHD simulations with a semiempirical FAL-C atmosphere to demonstrate the generation of acoustic waves in the interior in response to changing magnetic tilt in the corona. We find that Alfven-sound resonance combined with the ponderomotive force produces acoustic waves in the interior with sufficient energy to match sunquake observations when the magnetic field angle changes by the order of 10 degrees in a region where the coronal field strength is a few hundred gauss or more. The most energetic sunquakes are produced when the coronal field is strong, while the variation of magnetic field strength with height and the time scale of the tilt change are of secondary importance.

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.

We present the first magnetohydrodynamic model of the stellar chromospheric heating and acceleration of the outer atmospheres of cool evolved stars, using alpha Tau as a case study. We used a 1.5D MHD code with a generalized Ohm's law that accounts for the effects of partial ionization in the stellar atmosphere to study Alfven wave dissipation and wave reflection. We have demonstrated that due to inclusion of the effects of ion-neutral collisions in magnetized weakly ionized chromospheric plasma on resistivity and the appropriate grid resolution, the numerical resistivity becomes 1-2 orders of magnitude smaller than the physical resistivity. The motions introduced by non-linear transverse Alfven waves can explain non-thermally broadened and non-Gaussian profiles of optically thin UV lines forming in the stellar chromosphere of alpha Tau and other late-type giant and supergiant stars. The calculated heating rates in the stellar chromosphere due to resistive (Joule) dissipation of electric currents, induced by upward propagating non-linear Alfven waves, are consistent with observational constraints on the net radiative losses in UV lines and the continuum from alpha Tau. At the top of the chromosphere, Alfven waves experience significant reflection, producing downward propagating transverse waves that interact with upward propagating waves and produce velocity shear in the chromosphere. Our simulations also suggest that momentum deposition by non-linear Alfven waves becomes significant in the outer chromosphere at 1 stellar radius from the photosphere. The calculated terminal velocity and the mass loss rate are consistent with the observationally derived wind properties in alpha Tau.

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.

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 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.

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.

For a magnetic flux tube, or indeed any flux, to emerge into the Solar corona from the convection zone it must pass through the partially ionised layers of the lower atmosphere: the photosphere and the chromosphere. In such regions the ion-neutral collisions lead to an increased resistivity for currents flowing across magnetic field lines. This Cowling resistivity can exceed the Spitzer resistivity by orders of magnitude and in 2.5D simulations has been shown to be sufficient to remove all cross field current from emerging flux. Here we extend this modelling into 3D. Once again it is found that the Cowling resistivity removes perpendicular current. However the presence of 3D structure prevents the simple comparison possible in 2.5D simulations. With a fully ionised atmosphere the flux emergence leads to an unphysically low temperature region in the overlying corona, lifting of chromospheric material and the subsequent onset of the Rayleigh-Taylor instability. Including neutrals removes the low temperature region, lifts less chromospheric matter and shows no signs of the Rayleigh-Taylor instability. Simulations of flux emergence therefore should include such a neutral layer in order to obtain the correct perpendicular current, remove the Rayleigh-Taylor instability and get the correct temperature profile. In situations when the temperature is not important, i.e. when no simulated spectral emission is required, a simple model for the neutral layer is demonstrated to adequately reproduce the results of fully consistent simulations.

Certain regions of the solar atmosphere are at sufficiently low temperatures to be only partially ionised. The lower chromosphere contains neutral atoms, the existence of which greatly increases the efficiency of the damping of waves due to collisional friction momentum transfer. More specifically the Cowling conductivity can be up to 12 orders of magnitude smaller than the Spitzer value, so that the main damping mechanism in this region is due to the collisions between neutrals and positive ions. Using values for the gas density and temperature as functions of height taken from the VAL C model of the quiet Sun, an estimate is made for the dependance of the Cowling conductivity on height and strength of magnetic field. Using both analytic and numerical approaches the passage of Alfven waves over a wide spectrum through this partially ionised region is investigated. Estimates of the efficiency of this region in the damping of Alfven waves are made and compared for both approaches. We find that Alfven waves with frequencies above 0.6Hz are completely damped and frequencies below 0.01 Hz unaffected.