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

We present observations and magnetic field models of an intermediate filament present on the Sun in August 2012, associated with a polarity inversion line that extends from AR 11541 in the east into the quiet sun at its western end. A combination of SDO/AIA, SDO/HMI, and GONG H alpha data allow us to analyse the structure and evolution of the filament from 2012 August 4 23:00 UT to 2012 August 6 08:00 UT when the filament was in equilibrium. By applying the flux rope insertion method, nonlinear force-free field models of the filament are constructed using SDO/HMI line-of-sight magnetograms as the boundary condition at the two times given above. Guided by observed filament barbs, both modelled flux ropes are split into three sections each with a different value of axial flux to represent the non-uniform photospheric field distribution. The flux in the eastern section of the rope increases by 4$\times$10$^{20}$ Mx between the two models, which is in good agreement with the amount of flux cancelled along the internal PIL of AR 11541, calculated to be 3.2$\times$10$^{20}$ Mx. This suggests that flux cancellation builds flux into the filament's magnetic structure. Additionally, the number of field line dips increases between the two models in the locations where flux cancellation, the formation of new filament threads and growth of the filament is observed. This suggests that flux cancellation associated with magnetic reconnection forms concave-up magnetic field that lifts plasma into the filament. During this time, the free magnetic energy in the models increases by 0.2$\times$10$^{31}$ ergs.

Deriving the strength and direction of the three-dimensional (3D) magnetic field in the solar atmosphere is fundamental for understanding its dynamics. Volume information on the magnetic field mostly relies on coupling 3D reconstruction methods with photospheric and/or chromospheric surface vector magnetic fields. Infrared coronal polarimetry could provide additional information to better constrain magnetic field reconstructions. However, combining such data with reconstruction methods is challenging, e.g., because of the optical-thinness of the solar corona and the lack and limitations of stereoscopic polarimetry. To address these issues, we introduce the Data-Optimized Coronal Field Model (DOCFM) framework, a model-data fitting approach that combines a parametrized 3D generative model, e.g., a magnetic field extrapolation or a magnetohydrodynamic model, with forward modeling of coronal data. We test it with a parametrized flux rope insertion method and infrared coronal polarimetry where synthetic observations are created from a known "ground truth" physical state. We show that this framework allows us to accurately retrieve the ground truth 3D magnetic field of a set of force-free field solutions from the flux rope insertion method. In observational studies, the DOCFM will provide a means to force the solutions derived with different reconstruction methods to satisfy additional, common, coronal constraints. The DOCFM framework therefore opens new perspectives for the exploitation of coronal polarimetry in magnetic field reconstructions and for developing new techniques to more reliably infer the 3D magnetic fields that trigger solar flares and coronal mass ejections.

We study the properties of the local transverse deviations of magnetic field lines at a fixed moment in time. Those deviations "evolve" smoothly in a plane normal to the field-line direction as one moves that plane along the field line. Since the evolution can be described by a planar flow in the normal plane, we derive most of our results in the context of a toy model for planar fluid flow. We then generalize our results to include the effects of field-line curvature. We show that the type of flow is determined by the two non-zero eigenvalues of the gradient of the normalized magnetic field. The eigenvalue difference quantifies the local rate of squeezing or coiling of neighboring field lines, which we relate to standard notions of fluid vorticity and shear. The resulting squeezing rate can be used in the detection of null points, hyperbolic flux tubes and current sheets. Once integrated along field lines, that rate gives a squeeze factor, which is an approximation to the squashing factor, which is usually employed in locating quasi-separatrix layers (QSLs), which are possible sites for magnetic reconnection. Unlike the squeeze factor, the squashing factor can miss QSLs for which field lines are squeezed and then unsqueezed. In that regard, the squeeze factor is a better proxy for locating QSLs than the squashing factor. In another application of our analysis, we construct an approximation to the local rate of twist of neighboring field lines, which we refer to as the coiling rate. That rate can be integrated along a field line to give a coiling number, $\mathrm{N_c}$. We show that unlike the standard local twist number, $\mathrm{N_c}$ gives an unbiased approximation to the number of twists neighboring field lines make around one another. $\mathrm{N_c}$ can be useful for the study of flux rope instabilities, such as the kink instability, and can be used in the detection of flux ropes.

Quasi-Separatrix Layers (QSLs) are a useful proxy for the locations where current sheets can develop in the solar corona, and give valuable information about the connectivity in complicated magnetic field configurations. However, calculating QSL maps even for 2-dimensional slices through 3-dimensional models of coronal magnetic fields is a non-trivial task as it usually involves tracing out millions of magnetic field lines with immense precision. Thus, extending QSL calculations to three dimensions has rarely been done until now. In order to address this challenge, we present QSL Squasher -- a public, open-source code, which is optimized for calculating QSL maps in both two and three dimensions on GPUs. The code achieves large processing speeds for three reasons, each of which results in an order-of-magnitude speed-up. 1) The code is parallelized using OpenCL. 2) The precision requirements for the QSL calculation are drastically reduced by using perturbation theory. 3) A new boundary detection criterion between quasi-connectivity domains is used, which quickly identifies possible QSL locations which need to be finely sampled by the code. That boundary detection criterion relies on finding the locations of abrupt field-line length changes, which we do by introducing a new Field-line Length Edge (FLEDGE) map. We find FLEDGE maps useful on their own as a quick-and-dirty substitute for QSL maps. QSL Squasher allows constructing high-resolution 3D FLEDGE maps in a matter of minutes, which is two orders of magnitude faster than calculating the corresponding 3D QSL maps. We include a sample of calculations done using QSL Squasher to demonstrate its capabilities as a QSL calculator, as well as to compare QSL and FLEDGE maps.

Magnetic reconnection, the central engine that powers explosive phenomena throughout the Universe, is also perceived as one of the principal mechanisms for accelerating particles to high energies. Although various signatures of magnetic reconnection have been frequently reported, observational evidence that links particle acceleration directly to the reconnection site has been rare, especially for space plasma environments currently inaccessible to $\textit{in situ}$ measurements. Here we utilize broadband radio dynamic imaging spectroscopy available from the Karl G. Jansky Very Large Array to observe decimetric type III radio bursts in a solar jet with high angular ($\sim$20$''$), spectral ($\sim$1 %), and temporal resolution (50 milliseconds). These observations allow us to derive detailed trajectories of semi-relativistic (tens of keV) electron beams in the low solar corona with unprecedentedly high angular precision ($<0''.65$). We found that each group of electron beams, which corresponds to a cluster of type III bursts with 1-2-second duration, diverges from an extremely compact region ($\sim$600 km$^2$) in the low solar corona. The beam-diverging sites are located behind the erupting jet spire and above the closed arcades, coinciding with the presumed location of magnetic reconnection in the jet eruption picture supported by extreme ultraviolet/X-ray data and magnetic modeling. We interpret each beam-diverging site as a reconnection null point where multitudes of magnetic flux tubes join and reconnect. Our data suggest that the null points likely consist of a high level of density inhomogeneities possibly down to 10-km scales. These results, at least in the present case, strongly favor a reconnection-driven electron acceleration scenario.

Magnetic helicity is a quantity of great importance in solar studies because it is conserved in ideal magneto-hydrodynamics. While many methods to compute magnetic helicity in Cartesian finite volumes exist, in spherical coordinates, the natural coordinate system for solar applications, helicity is only treated approximately. We present here a method to properly compute relative magnetic helicity in spherical geometry. The volumes considered are finite, of shell or wedge shape, and the three-dimensional magnetic field is considered fully known throughout the studied domain. Testing of the method with well-known, semi-analytic, force-free magnetic-field models reveals that it has excellent accuracy. Further application to a set of nonlinear force-free reconstructions of the magnetic field of solar active regions, and comparison with an approximate method used in the past, indicates that the proposed methodology can be significantly more accurate, thus making our method a promising tool in helicity studies that employ the spherical geometry. Additionally, the range of applicability of the approximate method is determined and discussed.

We present the analysis of an unusual failed eruption captured in high cadence and in many wavelengths during the observing campaign in support of the VAULT2.0 sounding rocket launch. The refurbished Very high Angular resolution Ultraviolet Telescope (VAULT2.0) is a Ly$\alpha$ ($\lambda$ 1216 {\AA}) spectroheliograph launched on September 30, 2014. The campaign targeted active region NOAA AR 12172 and was closely coordinated with the Hinode and IRIS missions and several ground-based observatories (NSO/IBIS, SOLIS, and BBSO). A filament eruption accompanied by a low level flaring event (at the GOES C-class level) occurred around the VAULT2.0 launch. No Coronal Mass Ejection (CME) was observed. The eruption and its source region, however, were recorded by the campaign instruments in many atmospheric heights ranging from the photosphere to the corona in high cadence and spatial resolution. This is a rare occasion which enables us to perform a comprehensive investigation on a failed eruption. We find that a rising Magnetic Flux Rope-like (MFR) structure was destroyed during its interaction with the ambient magnetic field creating downflows of cool plasma and diffuse hot coronal structures reminiscent of "cusps". We employ magnetofrictional simulations to show that the magnetic topology of the ambient field is responsible for the destruction of the MFR. Our unique observations suggest that the magnetic topology of the corona is a key ingredient for a successful eruption.

Coronal jets represent important manifestations of ubiquitous solar transients, which may be the source of significant mass and energy input to the upper solar atmosphere and the solar wind. While the energy involved in a jet-like event is smaller than that of "nominal" solar flares and coronal mass ejections (CMEs), jets share many common properties with these phenomena, in particular, the explosive magnetically driven dynamics. Studies of jets could, therefore, provide critical insight for understanding the larger, more complex drivers of the solar activity. On the other side of the size-spectrum, the study of jets could also supply important clues on the physics of transients close or at the limit of the current spatial resolution such as spicules. Furthermore, jet phenomena may hint to basic process for heating the corona and accelerating the solar wind; consequently their study gives us the opportunity to attack a broad range of solar-heliospheric problems.

The standard model for eruptive flares has in the past few years been extended to 3D. It predicts typical J-shaped photospheric footprints of the coronal current layer, forming at similar locations as the Quasi-Separatrix Layers (QSLs). Such a morphology is also found for flare ribbons observed in the EUV band, as well as in non-linear force-free field (NLFFF) magnetic field extrapolations and models. We study the evolution of the photospheric traces of the current density and flare ribbons, both obtained with the SDO instruments. We investigate the photospheric current evolution during the 6 September 2011 X-class flare (SOL2011-09-06T22:20) from observational data of the magnetic field obtained with HMI. This evolution is compared with that of the flare ribbons observed in the EUV filters of the AIA. We also compare the observed electric current density and the flare ribbon morphology with that of the QSLs computed from the flux rope insertion method/NLFFF model. The NLFFF model shows the presence of a fan-spine configuration of overlying field lines, due to the presence of a parasitic polarity, embedding an elongated flux rope that appears in the observations as two parts of a filament. The QSLs, evolved via a magnetofrictional method, also show similar morphology and evolution as both the current ribbons and the EUV flare ribbons obtained at several times during the flare. For the first time, we propose a combined analysis of the photospheric traces of an eruptive flare, in a complex topology, with direct measurements of electric currents and QSLs from observational data and a magnetic field model. The results, obtained by two different and independent approaches, 1) confirm previous results of current increase during the impulsive phase of the flare, 2) show how NLFFF models can capture the essential physical signatures of flares even in a complex magnetic field topology.

In this paper we present a topological magnetic field investigation of seven two-ribbon flares in sigmoidal active regions observed with Hinode, STEREO, and SDO. We first derive the 3D coronal magnetic field structure of all regions using marginally unstable 3D coronal magnetic field models created with the flux rope insertion method. The unstable models have been shown to be a good model of the flaring magnetic field configurations. Regions are selected based on their pre-flare configurations along with the appearance and observational coverage of flare ribbons, and the model is constrained using pre-flare features observed in extreme ultraviolet and X-ray passbands. We perform a topology analysis of the models by computing the squashing factor, Q, in order to determine the locations of prominent quasi-separatrix layers (QSLs). QSLs from these maps are compared to flare ribbons at their full extents. We show that in all cases the straight segments of the two J-shaped ribbons are matched very well by the flux-rope-related QSLs, and the matches to the hooked segments are less consistent but still good for most cases. In addition, we show that these QSLs overlay ridges in the electric current density maps. This study is the largest sample of regions with QSLs derived from 3D coronal magnetic field models, and it shows that the magnetofrictional modeling technique that we employ gives a very good representation of flaring regions, with the power to predict flare ribbon locations in the event of a flare following the time of the model.

We present a new sample of M subdwarfs compiled from the 7th data release of the Sloan Digital Sky Survey. With 3517 new subdwarfs, this new sample significantly increases the number of spectroscopically confirmed low-mass subdwarfs. This catalog also includes 905 extreme and 534 ultra sudwarfs. We present the entire catalog including observed and derived quantities, and template spectra created from co-added subdwarf spectra. We show color-color and reduced proper motion diagrams of the three metallicity classes, which are shown to separate from the disk dwarf population. The extreme and ultra subdwarfs are seen at larger values of reduced proper motion as expected for more dynamically heated populations. We determine 3D kinematics for all of the stars with proper motions. The color-magnitude diagrams show a clear separation of the three metallicity classes with the ultra and extreme subdwarfs being significantly closer to the main sequence than the ordinary subdwarfs. All subdwarfs lie below (fainter) and to the left (bluer) of the main sequence. Based on the average $(U,V,W)$ velocities and their dispersions, the extreme and ultra subdwarfs likely belong to the Galactic halo, while the ordinary subdwarfs are likely part of the old Galactic (or thick) disk. An extensive activity analysis of subdwarfs is performed using H$\alpha$ emission and 208 active subdwarfs are found. We show that while the activity fraction of subdwarfs rises with spectral class and levels off at the latest spectral classes, consistent with the behavior of M dwarfs, the extreme and ultra subdwarfs are basically flat.

We present a statistical parallax study of nearly 2,000 M subdwarfs with photometry and spectroscopy from the Sloan Digital Sky Survey. Statistical parallax analysis yields the mean absolute magnitudes, mean velocities and velocity ellipsoids for homogenous samples of stars. We selected homogeneous groups of subdwarfs based on their photometric colors and spectral appearance. We examined the color-magnitude relations of low-mass subdwarfs and quantified their dependence on the newly-refined metallicity parameter, zeta. We also developed a photometric metallicity parameter, delta(g-r), based on the g-r and r-z colors of low-mass stars and used it to select stars with similar metallicities. The kinematics of low-mass subdwarfs as a function of color and metallicity were also examined and compared to main sequence M dwarfs. We find that the SDSS subdwarfs share similar kinematics to the inner halo and thick disk. The color-magnitude relations derived in this analysis will be a powerful tool for identifying and characterizing low-mass metal-poor subdwarfs in future surveys such as GAIA and LSST, making them important and plentiful tracers of the stellar halo.

We have observed the persistent but optically unidentified X-ray source X1908+075 with the PCA and HEXTE instruments on RXTE. The binary nature of this source was established by Wen, Remillard, & Bradt (2000) who found a 4.4-day orbital period in results from the RXTE ASM. We report the discovery of 605 s pulsations in the X-ray flux. The Doppler delay curve is measured and provides a mass function of 6.1 Msun which is a lower limit to the mass of the binary companion of the neutron star. The degree of attenuation of the low-energy end of the spectrum is found to be a strong function of orbital phase. A simple model of absorption in a stellar wind from the companion star fits the orbital phase dependence reasonably well and limits the orbital inclination angle to the range 38 to 72 degrees. These measured parameters lead to an orbital separation of 60 to 80 lt-s, a mass for the companion star in the range 9-31 Msun, and an upper limit to the size of the companion of ~22 Rsun. From our analysis we also infer a wind mass loss rate from the companion star of >~ 1.3 x 10^-6 Msun/yr and, when the properties of the companion star and the effects of photoionization are considered, likely >~ 4 x 10^-6 Msun/yr. Such a high rate is inconsistent with the allowed masses and radii that we find for a main sequence or modestly evolved star unless the mass loss rate is enhanced in the binary system relative to that of an isolated star. We discuss the possibility that the companion might be a Wolf-Rayet star that could evolve to become a black hole in 10^4 to 10^5 yr. If so, this would be the first identified progenitor of a neutron star--black hole binary.