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.

When magnetic flux emerges from beneath the photosphere it displaces the preexisting field in the corona, and a current sheet generally forms at the boundary between the old and new magnetic domains. Reconnection in the current sheet relaxes this highly stressed configuration to a lower energy state. This scenario is most familiar, and most often studied, in flares, where the flux transfer is rapid. We present here a study of steady, quiescent flux transfer occurring at a rate three orders of magnitude below that in a large flare. In particular we quantify the reconnection rate, and related energy release, occurring as new polarity emerges to form Active Region 11112 (SOL16 October 2010T00:00:00L205C117) within a region of preexisting flux. A bright, low lying kernel of coronal loops above the emerging polarity, observed with the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory and the X-ray Telescope onboard Hinode, originally shows magnetic connectivity only between regions of newly emerged flux when overlaid on magnetograms from the Helioseisimic and Magnetic Imager. Over the course of several days, this bright kernel advances into the preexisting flux. The advancement of an easily visible boundary into the old flux regions allows measurement of the rate of reconnection between old and new magnetic domains. We compare the reconnection rate to the inferred heating of the coronal plasma. To our knowledge, this is the first measurement of steady, quiescent heating related to reconnection. We determine that the newly emerged flux reconnects at a fairly steady rate of 0.38 × 10 16 Mx s −1 over two days, while the radiated power varies between (2 ∼ 8) × 10 25 erg s −1 over the same time. We find that as much as 40% of the total emerged flux at any given time may have reconnected. The total amount of transferred flux (∼ 1 × 10 21 Mx) and radiated energy (∼ 7.2 × 10 30 ergs) are comparable to that of a large M-or small X-class flare, but are stretched out over 45 hours.

The Minimum Current Corona (MCC) model provides a way to estimate stored coronal energy using the number of field lines connecting regions of positive and negative photospheric flux. This information is quantified by the net flux connecting pairs of opposing regions in a connectivity matrix. Changes in the coronal magnetic field, due to processes such as magnetic reconnection, manifest themselves as changes in the connectivity matrix. However, the connectivity matrix will also change when flux sources emerge or submerge through the photosphere, as often happens in active regions. We have developed an algorithm to estimate the changes in flux due to emergence and submergence of magnetic flux sources. These estimated changes must be accounted for in order to quantify storage and release of magnetic energy in the corona. To perform this calculation over extended periods of time, we must additionally have a consistently labeled connectivity matrix over the entire observational time span. We have therefore developed an automated tracking algorithm to generate a consistent connectivity matrix as the photospheric source regions evolve over time. We have applied this method to NOAA Active Region 11112, which underwent a GOES M-2.9 class flare around 19:00 on Oct.16 th , 2010, and calculated a lower bound on the free magnetic energy buildup of ∼ 8.25 × 10 30 ergs over 3 days.

In order to investigate the relation between magnetic structures and the signatures of heating in plage regions, we observed a plage region with the He i 1083.0 nm and Si i 1082.7 nm lines on 2018 October 3 using the integral field unit mode of the GREGOR Infrared Spectrograph (GRIS) installed at the GREGOR telescope. During the GRIS observation, the Interface Region Imaging Spectrograph obtained spectra of the ultraviolet Mg ii doublet emitted from the same region. In the periphery of the plage region, within the limited field of view seen by GRIS, we find that the Mg ii radiative flux increases with the magnetic field in the chromosphere with a factor of proportionality of 2.38 × 104 erg cm−2 s−1 G−1. The positive correlation implies that magnetic flux tubes can be heated by Alfvén wave turbulence or by collisions between ions and neutral atoms relating to Alfvén waves. Within the plage region itself, the radiative flux was large between patches of strong magnetic field strength in the photosphere or at the edges of magnetic patches. On the other hand, we do not find any significant spatial correlation between the enhanced radiative flux and the chromospheric magnetic field strength or the electric current. In addition to the Alfvén wave turbulence or collisions between ions and neutral atoms relating to Alfvén waves, other heating mechanisms related to magnetic field perturbations produced by interactions of magnetic flux tubes could be at work in the plage chromosphere.