One of the greatest challenges in solar physics is understanding the heating of the Sun's corona. Most theories for coronal heating postulate that free energy in the form of magnetic twist/stress is injected by the photosphere into the corona where the free energy is converted into heat either through reconnection or wave dissipation. The magnetic helicity associated with the twist/stress, however, is expected to be conserved and appear in the corona. In previous work we showed that helicity associated with the small-scale twists undergoes an inverse cascade via stochastic reconnection in the corona, and ends up as the observed large-scale shear of filament channels. Our ``helicity condensation'' model accounts for both the formation of filament channels and the observed smooth, laminar structure of coronal loops. In this paper, we demonstrate, using helicity- and energy-conserving numerical simulations of a coronal system driven by photospheric motions, that the model also provides a natural mechanism for heating the corona. We show that the heat generated by the reconnection responsible for the helicity condensation process is sufficient to account for the observed coronal heating. We study the role that helicity injection plays in determining coronal heating and find that, crucially, the heating rate is only weakly dependent on the net helicity preference of the photospheric driving. Our calculations demonstrate that motions with 100\% helicity preference are least efficient at heating the corona; those with 0\% preference are most efficient. We discuss the physical origins of this result and its implications for the observed corona.

Solar flares are observed and classified according to their intensity measured with the GOES X-ray Sensors. We show that the duration of a flare, as measured by the full width at half maximum (FWHM) in GOES is not related to the size of the flare as measured by GOES intensity. The durations of X-class flares range from a few minutes to a few hours, and the same is true of M- and C-class flares. In this work, we therefore examine the statistical relationships between the basic properties of flares -- temperature, emission measure, energy, etc. -- in comparison to both their size and duration. We find that the size of the flare is directly related to all of these basic properties, as previously found by many authors. The duration is not so clear. When examining the whole data set, the duration appears to be independent of all of these properties. In larger flares, however, there are direct correlations between the GOES FWHM and magnetic reconnection flux and ribbon area. We discuss the possible explanations, finding that this discrepancy may be due to large uncertainties in small flares, though we cannot rule out the possibility that the driving physical processes are different in smaller flares than larger ones. We discuss the implications of this result and how it relates to the magnetic reconnection process that releases energy in flares.

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.

Numerous observations have revealed that power-law distributions are ubiquitous in energetic solar processes. Hard X-rays, soft X-rays, extreme ultraviolet radiation, and radio waves all display power-law frequency distributions. Since magnetic reconnection is the driving mechanism for many energetic solar phenomena, it is likely that reconnection events themselves display such power-law distributions. In this work, we perform numerical simulations of the solar corona driven by simple convective motions at the photospheric level. Using temperature changes, current distributions, and Poynting fluxes as proxies for heating, we demonstrate that energetic events occurring in our simulation display power-law frequency distributions, with slopes in good agreement with observations. We suggest that the braiding-associated reconnection in the corona can be understood in terms of a self-organized criticality model driven by convective rotational motions similar to those observed at the photosphere.

The underlying origin of solar eruptive events (SEEs), ranging from giant coronal mass ejections to small coronal-hole jets, is that the lowest-lying magnetic flux in the Sun's corona undergoes the continual buildup of stress and free energy. This magnetic stress has long been observed as the phenomenon of "filament channels:" strongly sheared magnetic field localized around photospheric polarity inversion lines. However, the mechanism for the stress buildup - the formation of filament channels - is still debated. We present magnetohydrodynamic simulations of a coronal volume that is driven by transient, cellular boundary flows designed to model the processes by which the photosphere drives the corona. The key feature of our simulations is that they accurately preserve magnetic helicity, the topological quantity that is conserved even in the presence of ubiquitous magnetic reconnection. Although small-scale random stress is injected everywhere at the photosphere, driving stochastic reconnection throughout the corona, the net result of the magnetic evolution is a coherent shearing of the lowest-lying field lines. This highly counter-intuitive result - magnetic stress builds up locally rather than spreading out to attain a minimum energy state - explains the formation of filament channels and is the fundamental mechanism underlying SEEs. Furthermore, this mechanism may be relevant to other astrophysical or laboratory plasmas.

Two of the most widely observed and yet most puzzling features of the Sun's magnetic field are coronal loops that are smooth and laminar and prominences/filaments that are strongly sheared. These two features would seem to be quite unrelated in that the loops are near their minimum-energy current-free state, whereas filaments are regions of high magnetic stress and intense electric currents. We argue that, in fact, these two features are inextricably linked in that both are due to a single process: the injection of magnetic helicity into the corona by photospheric motions and the subsequent evolution of this helicity by coronal reconnection. In this paper, we present numerical simulations of the response of a \citet{Parker72} corona to photospheric driving motions that have varying degrees of helicity preference. We obtain four main conclusions: 1) in agreement with the helicity condensation model of \citet{Antiochos13}, the inverse cascade of helicity by magnetic reconnection results in the formation of prominences/filaments localized about polarity inversion lines (PILs); 2) this same process removes most structure from the rest of the corona, resulting in smooth and laminar coronal loops; 3) the amount of remnant tangling in coronal loops is inversely dependent on the net helicity injected by the driving motions; and 4) the structure of the solar corona depends only on the helicity preference of the driving motions and not on their detailed time dependence. We discuss the implications of our results for high-resolution observations of the corona.

Two of the most important features of the solar atmosphere are its hot, smooth coronal loops and the concentrations of magnetic shear, known as filament channels, that reside above photospheric polarity inversion lines (PILs). The shear observed in filament channels represents magnetic helicity, while the smoothness of the coronal loops indicates an apparent lack of magnetic helicity in the rest of the corona. At the same time, models that attempt to explain the high temperatures observed in these coronal loops require magnetic energy, in the form of twist, to be injected at the photosphere. In addition to magnetic energy, this twist also represents magnetic helicity. Unlike magnetic energy, magnetic helicity is conserved under reconnection, and is consequently expected to accumulate and be observed in the corona. However, filament channels, rather than the coronal loops, are the locations in the corona where magnetic helicity is observed, and it manifests itself in the form of shear, rather than twist. This naturally raises the question: if magnetic helicity needs to be injected to heat coronal loops, why is it only observed in filament channels, while coronal loops are observed to be laminar and smooth? This thesis addresses this question using a series of numerical simulations that demonstrate that magnetic helicity is transported throughout the solar corona by magnetic reconnection in such a way that it accumulates above PILs, forming filament channels, and leaving the rest of the corona generally smooth. In the process, it converts magnetic energy into heat, accounting for the large observed temperatures. This thesis presents a model for the formation of filament channels in the solar corona and the presence of smooth, hot coronal loops, and shows how the transport of magnetic helicity throughout the solar corona by magnetic reconnection is responsible for both of these phenomena.

We study nonlinear waves in a prominence foot using 2.5D MHD model motivated by recent high-resolution observations with Hinode/SOT in Ca~II emission of a prominence on October 10, 2012 showing highly dynamic small-scale motions in the prominence material. Observations of H$\alpha$ intensities and of Doppler shifts show similar propagating fluctuations. However the optically thick nature of the emission lines inhibits unique quantitative interpretation in terms of density. Nevertheless, we find evidence of nonlinear wave activity in the prominence foot by examining the relative magnitude of the fluctuation intensity ($\delta I/I\sim \delta n/n$). The waves are evident as significant density fluctuations that vary with height, and apparently travel upward from the chromosphere into the prominence material with quasi-periodic fluctuations with typical period in the range of 5-11 minutes, and wavelengths $\sim <$2000 km. Recent Doppler shift observations show the transverse displacement of the propagating waves. The magnetic field was measured with THEMIS instrument and was found to be 5-14 G. For the typical prominence density the corresponding fast magnetosonic speed is $\sim$20 km s$^{-1}$, in qualitative agreement with the propagation speed of the detected waves. The 2.5D MHD numerical model is constrained with the typical parameters of the prominence waves seen in observations. Our numerical results reproduce the nonlinear fast magnetosonic waves and provide strong support for the presence of these waves in the prominence foot. We also explore gravitational MHD oscillations of the heavy prominence foot material supported by dipped magnetic field structure.

A major unexplained feature of the solar atmosphere is the accumulation of magnetic shear, in the form of filament channels, at photospheric polarity inversion lines (PILs). In addition to free energy, this shear also represents magnetic helicity, which is conserved under reconnection. In this paper, we address the problem of filament channel formation and show how they acquire their shear and magnetic helicity. The results of 3D simulations using the Adaptively Refined Magnetohydrodynamics Solver (ARMS) are presented that support the model of filament channel formation by magnetic helicity condensation developed by \citet{Antiochos13}. We consider the supergranular twisting of a quasi-potential flux system that is bounded by a PIL and contains a coronal hole (CH). The magnetic helicity injected by the small-scale photospheric motions is shown to inverse-cascade up to the largest allowable scales that define the closed flux system: the PIL and the CH. This process produces field lines that are both sheared and smooth, and are sheared in opposite senses at the PIL and the CH. The accumulated helicity and shear flux are shown to be in excellent quantitative agreement with the helicity-condensation model. We present a detailed analysis of the simulations, including comparisons of our analytical and numerical results, and discuss their implications for observations.

On 20 August 2010 an energetic disturbance triggered large-amplitude longitudinal oscillations in a nearby filament. The triggering mechanism appears to be episodic jets connecting the energetic event with the filament threads. In the present work we analyze this periodic motion in a large fraction of the filament to characterize the underlying physics of the oscillation as well as the filament properties. The results support our previous theoretical conclusions that the restoring force of large-amplitude longitudinal oscillations is solar gravity, and the damping mechanism is the ongoing accumulation of mass onto the oscillating threads. Based on our previous work, we used the fitted parameters to determine the magnitude and radius of curvature of the dipped magnetic field along the filament, as well as the mass accretion rate onto the filament threads. These derived properties are nearly uniform along the filament, indicating a remarkable degree of cohesiveness throughout the filament channel. Moreover, the estimated mass accretion rate implies that the footpoint heating responsible for the thread formation, according to the thermal nonequilibrium model, agrees with previous coronal heating estimates. We estimate the magnitude of the energy released in the nearby event by studying the dynamic response of the filament threads, and discuss the implications of our study for filament structure and heating.

On 20 August 2010 an energetic disturbance triggered damped large-amplitude longitudinal (LAL) oscillations in almost an entire filament. In the present work we analyze this periodic motion in the filament to characterize the damping and restoring mechanism of the oscillation. Our method involves placing slits along the axis of the filament at different angles with respect to the spine of the filament, finding the angle at which the oscillation is clearest, and fitting the resulting oscillation pattern to decaying sinusoidal and Bessel functions. These functions represent the equations of motion of a pendulum damped by mass accretion. With this method we determine the period and the decaying time of the oscillation. Our preliminary results support the theory presented by Luna and Karpen (2012) that the restoring force of LAL oscillations is solar gravity in the tubes where the threads oscillate, and the damping mechanism is the ongoing accumulation of mass onto the oscillating threads. Following an earlier paper, we have determined the magnitude and radius of curvature of the dipped magnetic flux tubes hosting a thread along the filament, as well as the mass accretion rate of the filament threads, via the fitted parameters.

We report an unusual set of observations of waves in a large prominence pillar which consist of pulses propagating perpendicular to the prominence magnetic field. We observe a huge quiescent prominence with the Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA) in EUV on 2012 October 10 and only a part of it, the pillar, which is a foot or barb of the prominence, with the Hinode Solar Optical Telescope (SOT) (in Ca II and H\alpha lines), Sac Peak (in H\alpha, H\beta\ and Na-D lines), THEMIS ("T\'elescope H\'eliographique pour l' Etude du Magn\'etisme et des Instabilit\'es Solaires") with the MTR (MulTi-Raies) spectropolarimeter (in He D_3 line). The THEMIS/MTR data indicates that the magnetic field in the pillar is essentially horizontal and the observations in the optical domain show a large number of horizontally aligned features on a much smaller scale than the pillar as a whole. The data is consistent with a model of cool prominence plasma trapped in the dips of horizontal field lines. The SOT and Sac Peak data over the 4 hour observing period show vertical oscillations appearing as wave pulses. These pulses, which include a Doppler signature, move vertically, perpendicular to the field direction, along thin quasi-vertical columns in the much broader pillar. The pulses have a velocity of propagation of about 10 km/s, a period about 300 sec, and a wavelength around 2000 km. We interpret these waves in terms of fast magneto-sonic waves and discuss possible wave drivers.

The acceleration of solar flare ions during magnetic reconnection is explored via particle-in-cell simulations that self-consistently follow the motions of both protons and $\alpha$ particles. We demonstrate that the dominant ion heating during reconnection with a guide field (a magnetic component perpendicular to the reconnection plane) results from pickup behavior during the entry into reconnection exhausts. In contrast with anti-parallel reconnection, the temperature increment is dominantly transverse, rather than parallel, to the local magnetic field. The comparison of protons and alphas reveals a mass-to-charge ($M/Q$) threshold in pickup behavior that favors heating of high $M/Q$ ions over protons, which is consistent with impulsive flare observations.