The National Science Foundation’s Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand, and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities that will accompany full commissioning of the five facility instruments. With this Critical Science Plan (CSP) we attempt to anticipate some of what those capabilities will enable, providing a snapshot of some of the scientific pursuits that the DKIST hopes to engage as start-of-operations nears. The work builds on the combined contributions of the DKIST Science Working Group (SWG) and CSP Community members, who generously shared their experiences, plans, knowledge, and dreams. Discussion is primarily focused on those issues to which DKIST will uniquely contribute.

It is well known that photospheric flux emergence is an important process for stressing coronal fields and storing magnetic free energy, which may then be released during a flare. The Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) captured the entire emergence of NOAA AR 11158. This region emerged as two distinct bipoles, possibly connected underneath the photosphere, yet characterized by different photospheric field evolutions and fluxes. The combined active region complex produced 15 GOES C-class, 2 M-class, and the X2.2 Valentine's Day Flare during the four days after initial emergence on February 12th, 2011. The M and X class flares are of particular interest because they are nonhomologous, involving different subregions of the active region. We use a Magnetic Charge Topology together with the Minimum Current Corona model of the coronal field to model field evolution of the complex. Combining this with observations of flare ribbons in the 1600Å channel of the Atmospheric Imaging Assembly (AIA) onboard SDO, we propose a minimization algorithm for estimating the amount of reconnected flux and resulting drop in magnetic free energy during a flare. For the M6.6, M2.2, and X2.2 flares, we find a flux exchange of 4.2 × 10 20 Mx, 2.0 × 10 20 Mx, and 21.0 × 10 20 Mx, respectively, resulting in free energy drops of 3.89 × 10 30 ergs, 2.62 × 10 30 ergs, and 1.68 × 10 32 ergs.

Alfvénic 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.

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 produce 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én speed. Waves incident on the equipartition contour surrounding the null, where the sound and Alfvén 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.036E 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 regions of particular field orientation. We also model wavepacket propagation using a WKB method and find that it agrees qualitatively, though not quantitatively, with the results of the numerical simulation.

Using a simple two-dimensional, zero-β model, we explore the manner by which reconnection at a current sheet releases and dissipates free magnetic energy. We find that only a small fraction (3%-11% depending on currentsheet size) of the energy is stored close enough to the current sheet to be dissipated abruptly by the reconnection process. The remaining energy, stored in the larger-scale field, is converted to kinetic energy in a fast magnetosonic disturbance propagating away from the reconnection site, carrying the initial current and generating reconnectionassociated flows (inflow and outflow). Some of this reflects from the lower boundary (the photosphere) and refracts back to the X-point reconnection site. Most of this inward wave energy is reflected back again and continues to bounce between X-point and photosphere until it is gradually dissipated, over many transits. This phase of the energy dissipation process is thus global and lasts far longer than the initial purely local phase. In the process, a significant fraction of the energy (25%-60%) remains as undissipated fast magnetosonic waves propagating away from the reconnection site, primarily upward. This flare-generated wave is initiated by unbalanced Lorentz forces in the reconnection-disrupted current sheet, rather than by dissipation-generated pressure, as some previous models have assumed. Depending on the orientation of the initial current sheet, the wave front is either a rarefaction, with backward-directed flow, or a compression, with forward-directed flow.