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

We simulate a coronal mass ejection (CME) using a three-dimensional magnetohydrodynamic (MHD) code that includes coronal heating, thermal conduction, and radiative cooling in the energy equation. The magnetic flux distribution at 1 R$_s$ is produced by a localized subsurface dipole superimposed on a global dipole field, mimicking the presence of an active region within the global corona. Transverse electric fields are applied near the polarity inversion line to introduce a transverse magnetic field, followed by the imposition of a converging flow to form and destabilize a flux rope, producing an eruption. We examine the quantities responsible for plasma heating and cooling during the eruption, including thermal conduction, radiation, adiabatic effects, coronal heating, and ohmic heating. We find that ohmic heating is an important contributor to hot temperatures in the current sheet region early in the eruption, but in the late phase adiabatic compression plays an important role in heating the plasma there. Thermal conduction also plays an important role in the transport of thermal energy away from the current sheet region throughout the reconnection process, producing a ``thermal halo'' and widening the region of high temperatures. We simulate emission from solar telescopes for this eruption and find that there is evidence for emission from heated plasma above the flare loops late in the eruption, when the adiabatic heating is the dominant heating term. These results provide an explanation for hot supra-arcade plasma sheets that are often observed in X-rays and extreme ultraviolet wavelengths during the decay phase of large flares.

We present the observed local dispersion relations for magneto-acoustic-gravity waves in the Sun's atmosphere for different levels of magnetic field strength. We model these data with a theoretical local dispersion relation to produce spatial maps of the acoustic cut-off frequency in the Sun's photosphere. These maps have implications for the mechanical heating of the Sun's upper atmosphere, by magneto-acoustic-gravity waves, at different phases of the solar magnetic activity cycle.

(Abridged) In this white paper, we outline the major barriers to access within the educational and professional practice of astronomy. We present current best practices for inclusivity and accessibility, including classroom practices, institutional culture, support for infrastructure creation, hiring processes, and outreach initiatives. We present specific ways--beyond simple compliance with the ADA--that funding agencies, astronomers, and institutions can work together to make astronomy as a field more accessible, inclusive, and equitable. In particular, funding agencies should include the accessibility of institutions during proposal evaluation, hold institutions accountable for inaccessibility, and support efforts to gather data on the status and progress of astronomers and astronomy students with disabilities.

About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides in this intermediate region. We put forward that the next frontier for AMO science is to explore the AMO complexity that describes most of the Cosmos.

Magnetic reconnection is the most likely mechanism responsible for the high temperature events that are observed in strongly magnetized locations around the temperature minimum in the low solar chromosphere. This work improves upon our previous work ["Magnetic Reconnection in Strongly Magnetized Regions of the Low Solar Chromosphere", The Astrophysical Journal 852, 95 (2018) ] by using a more realistic radiative cooling model computed from the OPACITY project and the CHIANTI database. We find that the rate of ionization of the neutral component of the plasma is still faster than recombination within the current sheet region. For low $\beta$ plasmas, the ionized and neutral fluid flows are well-coupled throughout the reconnection region resembling the single-fluid Sweet-Parker model dynamics. Decoupling of the ion and neutral inflows appears in the higher $\beta$ case with $\beta_0=1.46$, which leads to a reconnection rate about three times faster than the rate predicted by the Sweet-Parker model. In all cases, the plasma temperature increases with time inside the current sheet, and the maximum value is above $2\times10^4$~K when the reconnection magnetic field strength is greater than $500$~G. While the more realistic radiative cooling model does not result in qualitative changes of the characteristics of magnetic reconnection, it is necessary for studying the variations of the plasma temperature and ionization fraction inside current sheets in strongly magnetized regions of the low solar atmosphere. It is also important for studying energy conversion during the magnetic reconnection process when the hydrogen-dominated plasma approaches full ionization.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

Magnetic reconnection in strongly magnetized regions around the temperature minimum region of the low solar atmosphere is studied by employing MHD-based simulations of a partially ionized plasma within a reactive 2.5D multi-fluid model. It is shown that in the absence of magnetic nulls in a low $\beta$ plasma the ionized and neutral fluid flows are well-coupled throughout the reconnection region. However, non-equilibrium ionization-recombination dynamics play a critical role in determining the structure of the reconnection region, lead to much lower temperature increases and a faster magnetic reconnection rate as compared to simulations that assume plasma to be in ionization-recombination equilibrium. The rate of ionization of the neutral component of the plasma is always faster than recombination within the current sheet region even when the initial plasma $\beta$ is as high as $\beta_0=1.46$. When the reconnecting magnetic field is in excess of a kilogauss and the plasma $\beta$ is lower than 0.0145, the initially weakly ionized plasmas can become fully ionized within the reconnection region and the current sheet can be strongly heated to above $2.5\times10^4$~K, even as most of the collisionally dissipated magnetic energy is radiated away. The Hall effect increases the reconnection rate slightly, but in the absence of magnetic nulls it does not result in significant asymmetries or change the characteristics of the reconnection current sheet down to meter scales.

We perform two-dimensional resistive magnetohydrodynamic simulations of coronal jets driven by flux emergence along the lower boundary. The reconnection layers are susceptible the formation of blobs that are ejected in the jet. Our simulation with low plasma $\beta$ (Case I) shows that magnetic islands form easily and propagate upwards in the jet. These islands are multithermal and thus are predicted to show up in hot channels (335 \AA\, and 211 \AA) and the cool channel (304 \AA) in observations by the Atmospheric Imaging Assembly (AIA) on the \emph{Solar Dynamics Observatory}. The islands have maximum temperatures of 8 MK, lifetimes of 120 s, diameters of 6 Mm, and velocities of 200 km s$^{-1}$. These parameters are similar to the properties of blobs observed in EUV jets by AIA\@. The Kelvin-Helmholtz instability develops in our simulation with moderately high plasma $\beta$ (Case II), and leads to the formation of bright vortex-like blobs above the multiple high magneto-sonic Mach number regions that appear along the jet. These vortex-like blobs can also be identified in the AIA channels. However, they eventually move downward and disappear after the high magneto-sonic Mach number regions disappear. In the lower plasma $\beta$ case, the lifetime for the jet is shorter, the jet and magnetic islands are formed with higher velocities and temperatures, the current sheet fragments are more chaotic, and more magnetic islands are generated. Our results show that the plasmoid instability and Kelvin-Helmholtz instability along the jet are both possible causes of the formation of blobs observed at extreme ultraviolet (EUV) wavelengths.

Starting in December 2014, Kepler K2 observed Neptune continuously for 49 days at a 1-minute cadence. The goals consisted of studying its atmospheric dynamics (Simon et al. 2016), detecting its global acoustic oscillations (Rowe et al., submitted), and those of the Sun, which we report on here. We present the first indirect detection of solar oscillations in intensity measurements. Beyond the remarkable technical performance, it indicates how Kepler would see a star like the Sun. The result from the global asteroseismic approach, which consists of measuring the oscillation frequency at maximum amplitude "nu_max" and the mean frequency separation between mode overtones "Delta nu", is surprising as the nu_max measured from Neptune photometry is larger than the accepted value. Compared to the usual reference nu_max_sun = 3100 muHz, the asteroseismic scaling relations therefore make the solar mass and radius appear larger by 13.8 +/- 5.8 % and 4.3 +/- 1.9 % respectively. The higher nu_max is caused by a combination of the value of nu_max_sun, being larger at the time of observations than the usual reference from SOHO/VIRGO/SPM data (3160 +/- 10 muHz), and the noise level of the K2 time series, being ten times larger than VIRGO's. The peak-bagging method provides more consistent results: despite a low signal-to-noise ratio (SNR), we model ten overtones for degrees l=0,1,2. We compare the K2 data with simultaneous SOHO/VIRGO/SPM photometry and BiSON velocity measurements. The individual frequencies, widths, and amplitudes mostly match those from VIRGO and BiSON within 1 sigma, except for the few peaks with lowest SNR.

While theoretical models and simulations of magnetic reconnection often assume symmetry such that the magnetic null point when present is co-located with a flow stagnation point, the introduction of asymmetry typically leads to non-ideal flows across the null point. To understand this behavior, we present exact expressions for the motion of three-dimensional linear null points. The most general expression shows that linear null points move in the direction along which the vector field and its time derivative are antiparallel. Null point motion in resistive magnetohydrodynamics results from advection by the bulk plasma flow and resistive diffusion of the magnetic field, which allows non-ideal flows across topological boundaries. Null point motion is described intrinsically by parameters evaluated locally; however, global dynamics help set the local conditions at the null point. During a bifurcation of a degenerate null point into a null-null pair or the reverse, the instantaneous velocity of separation or convergence of the null-null pair will typically be infinite along the null space of the Jacobian matrix of the magnetic field, but with finite components in the directions orthogonal to the null space. Not all bifurcating null-null pairs are connected by a separator. Furthermore, except under special circumstances, there will not exist a straight line separator connecting a bifurcating null-null pair. The motion of separators cannot be described using solely local parameters, because the identification of a particular field line as a separator may change as a result of non-ideal behavior elsewhere along the field line.

Realistic models of magnetic reconnection in the solar chromosphere must take into account that the plasma is partially ionized and that plasma conditions within any two magnetic flux bundles undergoing reconnection may not be the same. Asymmetric reconnection in the chromosphere may occur when newly emerged flux interacts with pre-existing, overlying flux. We present 2.5D simulations of asymmetric reconnection in weakly ionized, reacting plasmas where the magnetic field strengths, ion and neutral densities, and temperatures are different in each upstream region. The plasma and neutral components are evolved separately to allow non-equilibrium ionization. As in previous simulations of chromospheric reconnection, the current sheet thins to the scale of the neutral-ion mean free path and the ion and neutral outflows are strongly coupled. However, the ion and neutral inflows are asymmetrically decoupled. In cases with magnetic asymmetry, a net flow of neutrals through the current sheet from the weak field (high density) upstream region into the strong field upstream region results from a neutral pressure gradient. Consequently, neutrals dragged along with the outflow are more likely to originate from the weak field region. The Hall effect leads to the development of a characteristic quadrupole magnetic field modified by asymmetry, but the X-point geometry expected during Hall reconnection does not occur. All simulations show the development of plasmoids after an initial laminar phase.

As the interface between the Sun's photosphere and corona, the chromosphere and transition region play a key role in the formation and acceleration of the solar wind. Observations from the Interface Region Imaging Spectrograph reveal the prevalence of intermittent small-scale jets with speeds of 80-250 km/s from the narrow bright network lanes of this interface region. These jets have lifetimes of 20-80 seconds and widths of 300 km or less. They originate from small-scale bright regions, often preceded by footpoint brightenings and accompanied by transverse waves with ~20 km/s amplitudes. Many jets reach temperatures of at least ~100000 K and constitute an important element of the transition region structures. They are likely an intermittent but persistent source of mass and energy for the solar wind.

Magnetic reconnection is believed to be the dominant energy release mechanism in solar flares. The standard flare model predicts both downward and upward outflow plasmas with speeds close to the coronal Alfv\'{e}n speed. Yet, spectroscopic observations of such outflows, especially the downflows, are extremely rare. With observations of the newly launched Interface Region Imaging Spectrograph (IRIS), we report the detection of greatly redshifted ($\sim$125 km s$^{-1}$ along line of sight) Fe {\sc{xxi}} 1354.08\AA{} emission line with a $\sim$100 km s$^{-1}$ nonthermal width at the reconnection site of a flare. The redshifted Fe {\sc{xxi}} feature coincides spatially with the loop-top X-Ray source observed by the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI). We interpret this large redshift as the signature of downward-moving reconnection outflow/hot retracting loops. Imaging observations from both IRIS and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) also reveal the eruption and reconnection processes. Fast downward-propagating blobs along these loops are also found from cool emission lines (e.g., Si {\sc{iv}}, O {\sc{iv}}, C {\sc{ii}}, Mg {\sc{ii}}) and images of AIA and IRIS. Furthermore, the entire Fe {\sc{xxi}} line is blueshifted by $\sim$260 km s$^{-1}$ at the loop footpoints, where the cool lines mentioned above all exhibit obvious redshift, a result that is consistent with the scenario of chromospheric evaporation induced by downward-propagating nonthermal electrons from the reconnection site.

We propose a revolutionary way of studying the sur-face of Mars using a wind-driven network of mobile sensors- Gone with the Wind ON_Mars (GOWON). GOWON is envisioned to be a scalable, 100% self energy-generating and distributed system that allows in-situ mapping of a wide range of phenomena in a much larger portion of the surface of Mars compared to earlier missions. It could radically improve the possibility of finding rare phenomena like bio signatures through random wind-driven search. It could explore difficult terrains that were beyond the reach of previous missions, such as regions with very steep slopes, cluttered surfaces and/or sand dunes; GOWON is envisioned as an on going mission with a long life span. It could achieve any of NASA's scientific objectives on Mars in a cost-effective way, leaving a long lasting sensing and searching infrastructure on Mars. GOWON is a 2012 Step B invitee for NASA Innovative Advanced Concept (NIAC). It addresses the challenge area of the Mars Surface System Capabilities area. We believe the challenge to be near-term, i.e., 2018-2024.

We present two-dimensional resistive magnetohydrodynamic simulations of line-tied asymmetric magnetic reconnection in the context of solar flare and coronal mass ejection current sheets. The reconnection process is made asymmetric along the inflow direction by allowing the initial upstream magnetic field strengths and densities to differ, and along the outflow direction by placing the initial perturbation near a conducting wall boundary that represents the photosphere. When the upstream magnetic fields are asymmetric, the post-flare loop structure is distorted into a characteristic skewed candle flame shape. The simulations can thus be used to provide constraints on the reconnection asymmetry in post-flare loops. More hard X-ray emission is expected to occur at the footpoint on the weak magnetic field side because energetic particles are more likely to escape the magnetic mirror there than at the strong magnetic field footpoint. The footpoint on the weak magnetic field side is predicted to move more quickly because of the requirement in two dimensions that equal amounts of flux must be reconnected from each upstream region. The X-line drifts away from the conducting wall in all simulations with asymmetric outflow and into the strong magnetic field region during most of the simulations with asymmetric inflow. There is net plasma flow across the X-line for both the inflow and outflow directions. The reconnection exhaust directed away from the obstructing wall is significantly faster than the exhaust directed towards it. The asymmetric inflow condition allows net vorticity in the rising outflow plasmoid which would appear as rolling motions about the flux rope axis.

There are relatively few observations of UV emission during the impulsive phases of solar flares, so the nature of that emission is poorly known. Photons produced by solar flares can resonantly scatter off atoms and ions in the corona. Based on off-limb measurements by SOHO/UVCS, we derive the O VI $\lambda$1032 luminosities for 29 flares during the impulsive phase and the Ly$\alpha$ luminosities of 5 flares, and we compare them with X-ray luminosities from GOES measurements. The upper transition region and lower transition region luminosities of the events observed are comparable. They are also comparable to the luminosity of the X-ray emitting gas at the beginning of the flare, but after 10-15 minutes the X-ray luminosity usually dominates. In some cases we can use Doppler dimming to estimate flow speeds of the O VI emitting gas, and 5 events show speeds in the 40 to 80 $\rm km s^{-1}$ range. The O VI emission could originate in gas evaporating to fill the X-ray flare loops, in heated chromospheric gas at the footpoints, or in heated prominence material in the coronal mass ejection. All three sources may contribute in different events or even in a single event, and the relative timing of UV and X-ray brightness peaks, the flow speeds, and the total O VI luminosity favor each source in one or more events.

We perform a time-dependent ionization analysis to constrain plasma heating requirements during a fast partial halo coronal mass ejection (CME) observed on 2000 June 28 by the Ultraviolet Coronagraph Spectrometer (UVCS) aboard the Solar and Heliospheric Observatory (SOHO). We use two methods to derive densities from the UVCS measurements, including a density sensitive O V line ratio at 1213.85 and 1218.35 Angstroms, and radiative pumping of the O VI 1032,1038 doublet by chromospheric emission lines. The most strongly constrained feature shows cumulative plasma heating comparable to or greater than the kinetic energy, while features observed earlier during the event show cumulative plasma heating of order or less than the kinetic energy. SOHO Michelson Doppler Imager (MDI) observations are used to estimate the active region magnetic energy. We consider candidate plasma heating mechanisms and provide constraints when possible. Because this CME was associated with a relatively weak flare, the contribution by flare energy (e.g., through thermal conduction or energetic particles) is probably small; however, the flare may have been partially behind the limb. Wave heating by photospheric motions requires heating rates significantly larger than those previously inferred for coronal holes, but the eruption itself could drive waves which heat the plasma. Heating by small-scale reconnection in the flux rope or by the CME current sheet is not significantly constrained. UVCS line widths suggest that turbulence must be replenished continually and dissipated on time scales shorter than the propagation time in order to be an intermediate step in CME heating.

Understanding the physical processes responsible for accelerating the solar wind requires detailed measurements of the collisionless plasma in the extended solar corona. Some key clues about these processes have come from instruments that combine the power of an ultraviolet (UV) spectrometer with an occulted telescope. This combination enables measurements of ion emission lines far from the bright solar disk, where most of the solar wind acceleration occurs. Although the UVCS instrument on SOHO made several key discoveries, many questions remain unanswered because its capabilities were limited. This white paper summarizes these past achievements and also describes what can be accomplished with next-generation instrumentation of this kind.

To investigate the impact of current sheet motion on the reconnection process, we perform resistive magnetohydrodynamic (MHD) simulations of two closely located reconnection sites which move apart from each other as reconnection develops. This simulation develops less quickly than an otherwise equivalent single perturbation simulation but eventually exhibits a higher reconnection rate. The unobstructed outflow jets are faster and longer than the outflow jets directed towards the magnetic island that forms between the two current sheets. The X-line and flow stagnation point are located near the trailing end of each current sheet very close to the obstructed exit. The speed of X-line retreat ranges from ~0.02-0.06 while the speed of stagnation point retreat ranges from ~0.03-0.07, in units of the initial upstream Alfven velocity. Early in time, the flow stagnation point is located closer to the center of the current sheet than the X-line, but later on the relative positions of these two points switch. Consequently, late in time there is significant plasma flow across the X-line in the opposite direction of X-line retreat. Throughout the simulation, the velocity at the X-line does not equal the velocity of the X-line. Motivated by these results, an expression for the rate of X-line retreat is derived in terms of local parameters at the X-point. This expression shows that X-line retreat is due to both advection by the bulk plasma flow and diffusion of the normal component of the magnetic field.

The POLAR Investigation of the Sun (POLARIS) mission uses a combination of a gravity assist and solar sail propulsion to place a spacecraft in a 0.48 AU circular orbit around the Sun with an inclination of 75 degrees with respect to solar equator. This challenging orbit is made possible by the challenging development of solar sail propulsion. This first extended view of the high-latitude regions of the Sun will enable crucial observations not possible from the ecliptic viewpoint or from Solar Orbiter. While Solar Orbiter would give the first glimpse of the high latitude magnetic field and flows to probe the solar dynamo, it does not have sufficient viewing of the polar regions to achieve POLARIS' primary objective : determining the relation between the magnetism and dynamics of the Sun's polar regions and the solar cycle.

It is now well established that stellar winds of hot stars are fragmentary and that the X-ray emission from stellar winds has a strong contribution from shocks in winds. Chandra high spectral resolution observations of line profiles of O and B stars have shown numerous properties that had not been expected. Here we suggest explanations by considering the X-rays as arising from bow shocks that occur where the stellar wind impacts on spherical clumps in the winds. We use an accurate and stable numerical hydrodynamical code to obtain steady-state physical conditions for the temperature and density structure in a bow shock. We use these solutions plus analytic approximations to interpret some major X-ray features: the simple power-law distribution of the observed emission measure derived from many hot star X-ray spectra and the wide range of ionization stages that appear to be present in X-ray sources throughout the winds. Also associated with the adiabatic cooling of the gas around a clump is a significant transverse velocity for the hot plasma flow around the clumps, and this can help to understand anomalies associated with observed line widths, and the differences in widths seen in stars with high and low mass-loss rates. The differences between bow shocks and the planar shocks that are often used for hot stars are discussed. We introduce an ``on the shock'' (OTSh) approximation that is useful for interpreting the X-rays and the consequences of clumps in hot star winds and elsewhere in astronomy.

We report the results of a high angular resolution near-infrared survey of dusty Wolf-Rayet stars using the Keck-1 Telescope, including new multi-wavelength images of the pinwheel nebulae WR 98a, WR 104, and WR 112. Angular sizes were measured for an additional 8 dusty WR stars using aperture masking interferometry, allowing us to probe characteristics sizes down to ~20 milliarcseconds (~40 AU for typical sources). With angular sizes and specific fluxes, we can directly measure the wavelength-dependent surface brightness and size relations for our sample. We discovered tight correlations of these properties within our sample which could not be explained by simple spherically-symmetric dust shells or even the more realistic ``pinwheel nebula'' (3-D) radiative transfer model, when using optical constants of Zubko. While the tightly-correlated surface brightness relations we uncovered offer compelling indirect evidence of a shared and distinctive dust shell geometry amongst our sample, long-baseline interferometers should target the marginally-resolved objects in our sample in order to conclusively establish the presence or absence of the putative underyling colliding wind binaries thought to produce the dust shells around WC Wolf-Rayets.

Brown dwarfs are classified as objects which are not massive enough to sustain nuclear fusion of hydrogen, and are distinguished from planets by their ability to burn deuterium. Old (>10 Myr) brown dwarfs are expected to possess short-lived magnetic fields and, since they no longer generate energy from collapse and accretion, weak radio and X-ray emitting coronae. Several efforts have been undertaken in the past to detect chromospheric activity from the brown dwarf LP944-20 at X-ray and optical wavelengths, but only recently an X-ray flare from this object was detected. Here we report on the discovery of quiescent and flaring radio emission from this source, which represents the first detection of persistent radio emission from a brown dwarf, with luminosities that are several orders of magnitude larger than predicted from an empirical relation between the X-ray and radio luminosities of many stellar types. We show in the context of synchrotron emission, that LP944-20 possesses an unusually weak magnetic field in comparison to active dwarf M stars, which might explain the null results from previous optical and X-ray observations of this source, and the deviation from the empirical relations.