Flare ribbons are bright manifestations of flare energy dissipation in the lower solar atmosphere. For the first time, we report on high-resolution imaging spectroscopy observations of flare ribbons situated off-limb in the H$\beta$ and Ca II 8542 {\AA} lines and make a detailed comparison with radiative hydrodynamic simulations. Observations of the X8.2-class solar flare SOL2017-09-10T16:06 UT obtained with the Swedish Solar Telescope reveal bright horizontal emission layers in H$\beta$ line wing images located near the footpoints of the flare loops. The apparent separation between the ribbon observed in the H$\beta$ wing and the nominal photospheric limb is about 300 - 500 km. The Ca II 8542 {\AA} line wing images show much fainter ribbon emissions located right on the edge of the limb, without clear separation from the limb. RADYN models are used to investigate synthetic spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. The simulations show that, towards the limb, where the line of sight is substantially oblique with respect to the vertical direction, the flaring atmosphere model reproduces the high contrast of the off-limb H$\beta$ ribbons and their significant elevation above the photosphere. The ribbons in the Ca II 8542 {\AA} line wing images are located deeper in the lower solar atmosphere with a lower contrast. A comparison of the height deposition of electron beam energy and the intensity contribution function shows that the H$\beta$ line wing intensities can be an useful tracer of flare energy deposition in the lower solar atmosphere

We report multi-wavelength monitoring observations of an M-dwarf flare star AD Leonis with Seimei Telescope (6150--7930 {\AA}), SCAT (Spectroscopic Chuo-university Astronomical Telescope; 3700--7500 {\AA}), NICER (Neutron Star Interior Composition Explorer; 0.2--12.0 keV), and collaborations of OISTER (Optical and Infrared Synergetic Telescopes for Education and Research) program. Twelve flares are detected in total which include ten H$\alpha$, four X-ray, and four optical-continuum flares; one of them is a superflare with the total energy of $\sim$ 2.0$\times$10$^{33}$ erg. We found that (1) during the superflare, the H$\alpha$ emission line full width at 1/8 maximum dramatically increases to 14 {\AA} from 8 {\AA} in the low-resolution spectra (R$\sim$ 2000) accompanied with the large white-light flares, (2) some weak H$\alpha$/X-ray flares are not accompanied with white-light emissions, and (3) the non-flaring emissions show clear rotational modulations in X-ray and H$\alpha$ intensity in the same phase. To understand these observational features, one-dimensional hydrodynamic flare simulations are performed by using the RADYN code. As a result of simulations, we found the simulated H$\alpha$ line profiles with hard and high-energy non-thermal electron beams are consistent with that of the initial phase line profiles of the superflares, while those with more soft- and/or weak-energy beam are consistent with those in decay phases, indicating the changes in the energy fluxes injected to the lower atmosphere. Also, we found that the relation between optical continuum and H$\alpha$ intensity is nonlinear, which can be one cause of the non-white-light flares. The flare energy budget exhibits diversity in the observations and models, and more observations of stellar flares are necessary for constraining the occurrence of various emission line phenomena in stellar flares.

We study the evolution of chromospheric line and continuum emission during the impulsive phase of the X-class SOL2014-09-10T17:45 solar flare. We extend previous analyses of this flare to multiple chromospheric lines of Fe I, Fe II, Mg II, C I, and Si II, observed with IRIS, combined with radiative-hydrodynamical (RHD) modeling. For multiple flaring kernels, the lines all show a rapidly evolving double-component structure: an enhanced, emission component at rest, and a broad, highly red-shifted component of comparable intensity. The red-shifted components migrate from 25-50 km s$^{-1}$ towards the rest wavelength within $\sim$30 seconds. Using Fermi hard X-ray observations, we derive the parameters of an accelerated electron beam impacting the dense chromosphere, using them to drive a RHD simulation with the RADYN code. As in Kowalski et al. 2017a, our simulations show that the most energetic electrons penetrate into the deep chromosphere, heating it to T$\sim$10,000 K, while the bulk of the electrons dissipate their energy higher, driving an explosive evaporation, and its counterpart condensation -- a very dense (n$_e \sim 2 \times 10^{14}$ cm$^{-3}$), thin layer (30--40 km thickness), heated to 8--12,000 K, moving towards the stationary chromosphere at up to 50 km s$^{-1}$. The synthetic Fe II 2814.45A profiles closely resemble the observational data, including a continuum enhancement, and both a stationary and a highly red-shifted component, rapidly moving towards the rest wavelength. Importantly, the absolute continuum intensity, ratio of component intensities, relative time of appearance, and red-shift amplitude, are sensitive to the model input parameters, showing great potential as diagnostics.

To extract the information that the Mg II NUV spectra (observed by the Interface Region Imaging Spectrograph; IRIS), carries about the chromosphere during solar flares, and to validate models of energy transport via model-data comparison, forward modelling is required. The assumption of statistical equilibrium is typically used to obtain the atomic level populations from snapshots of flare atmospheres, due to computational necessity. However it is possible that relying on statistical equilibrium could lead to spurious results. We compare solving the atomic level populations via statistical equilibrium versus a non-equilibrium time-dependent approach. This was achieved using flare simulations from RADYN alongside the minority species version, MS_RADYN, from which the time-dependent Mg II atomic level populations and radiation transfer were computed in complete frequency redistribution. The impacts on the emergent profiles, lightcurves, line ratios, and formation heights are discussed. In summary we note that non-equilibrium effects during flares are typically important only in the initial stages and for a short period following the cessation of the energy injection. An analysis of the timescales of ionisation equilibrium reveals that for most of the duration of the flare, when the temperatures and densities are sufficiently enhanced, the relaxation timescales are short ($\tau_{\mathrm{relax}}<0.1$ s), so that the equilibrium solution is an adequate approximation. These effects vary with the size of the flare, however. In weaker flares effects can be more pronounced. We recommend that non-equilibrium effects be considered when possible, but that statistical equilibrium is sufficient at most stages of the flare.

The Interface Region Imaging Spectrograph (IRIS) has routinely observed the flaring Mg II NUV spectrum, offering excellent diagnostic potential and a window into the location of energy deposition. A number of studies have forward modelled both the general properties of these lines and specific flare observations. Generally these have forward modelled radiation via post-processing of snapshots from hydrodynamic flare simulations through radiation transfer codes. There has, however, not been a study of how the physics included in these radiation transport codes affects the solution. A baseline setup for forward modelling MgII in flares is presented and contrasted with approaches that add or remove complexity. It is shown for Mg II: (1) PRD is still required during flare simulations despite the increased densities, (2) using full angle-dependent PRD affects the solution but takes significantly longer to process a snapshot, (3) including Mg I in NLTE results in negligible differences to the Mg II lines but does affect the NUV quasi-continuum, (4) only hydrogen and Mg II need to be included in NLTE, (5) ideally the non-equilibrium hydrogen populations, with non-thermal collisional rates, should be used rather than the statistical equilibrium populations, (6) an atom consisting of only the ground state, h & k upper levels, and continuum level is insufficient to model the resonance lines, and (7) irradiation from a hot, dense flaring transition region can affect the formation of Mg II. We discuss modifications to the RH code allowing straightforward inclusion of transition region and coronal irradiation in flares.

The GOES X1 flare SOL2014-10-25T17:08:00 was a three-ribbon solar flare observed with IRIS in the near and far ultraviolet. One of the flare ribbons crossed a sunspot umbra, producing a dramatic, $\sim1000$% increase in the near-ultraviolet (NUV) continuum radiation. We comprehensively analyze the ultraviolet spectral data of the umbral flare brightenings, which provide new challenges for radiative-hydrodynamic modeling of the chromospheric velocity field and the white-light continuum radiation. The emission line profiles in the umbral flare brightenings exhibit redshifts and profile asymmetries, but these are significantly smaller than in another, well-studied X-class solar flare. We present a ratio of the NUV continuum intensity to the Fe II 2814.45 Ang intensity. This continuum-to-line ratio is a new spectral diagnostic of significant heating at high column mass (log $m/$[g cm$^{-2}] >-2$) during solar flares because the continuum and emission line radiation originate from relatively similar temperatures but moderately different optical depths. The full spectral readout of these IRIS data also allow for a comprehensive survey of the flaring NUV landscape: in addition to many lines of Fe II and Cr II, we identify a new solar flare emission line, He I $\lambda2829.91$ (as previously identified in laboratory and early-type stellar spectra). The Fermi/GBM hard X-ray data provide inputs to radiative-hydrodynamic models (which will be presented in Paper II) in order to better understand the large continuum-to-line ratios, the origin of the white-light continuum radiation, and the role of electron beam heating in the low atmosphere.

Observations from the \textit{Interface Region Imaging Spectrograph} (\textsl{IRIS}) often reveal significantly broadened and non-reversed profiles of the Mg II h, k and triplet lines at flare ribbons. To understand the formation of these optically thick Mg II lines, we perform plane parallel radiative hydrodynamics modeling with the RADYN code, and then recalculate the Mg II line profiles from RADYN atmosphere snapshots using the radiative transfer code RH. We find that the current RH code significantly underestimates the Mg II h \& k Stark widths. By implementing semi-classical perturbation approximation results of quadratic Stark broadening from the STARK-B database in the RH code, the Stark broadenings are found to be one order of magnitude larger than those calculated from the current RH code. However, the improved Stark widths are still too small, and another factor of 30 has to be multiplied to reproduce the significantly broadened lines and adjacent continuum seen in observations. Non-thermal electrons, magnetic fields, three-dimensional effects or electron density effect may account for this factor. Without modifying the RADYN atmosphere, we have also reproduced non-reversed Mg II h \& k profiles, which appear when the electron beam energy flux is decreasing. These profiles are formed at an electron density of $\sim 8\times10^{14}\ \mathrm{cm}^{-3}$ and a temperature of $\sim1.4\times10^4$ K, where the source function slightly deviates from the Planck function. Our investigation also demonstrates that at flare ribbons the triplet lines are formed in the upper chromosphere, close to the formation heights of the h \& k lines.

The EUV (100-912 {\AA}) is a spectral region notoriously difficult to observe due to attenuation by neutral hydrogen gas in the interstellar medium. Despite this, hundreds to thousands of nearby stars of different spectral types and magnetic activity levels are accessible in the EUV range. The EUV probes interesting and complicated regions in the stellar atmosphere like the lower corona and transition region that are inaccessible from other spectral regions. In this white paper we describe how direct EUV observations, which require a dedicated grazing-incidence observatory, cannot yet be accurately substituted with models and theory. Exploring EUV emission from cool dwarf stars in the time domain can make a major contribution to understanding stellar outer atmospheres and magnetism, and offers the clearest path toward detecting coronal mass ejections on stars other than the Sun.

The planetary effective surface temperature alone is insufficient to characterize exoplanet atmospheres and their stability or evolution. Considering the star-planet system as a whole is necessary, and a critical component of the system is the photoionizing stellar extreme ultraviolet emission (EUV; 100-912 {\AA}). EUV photons drive atmospheric mass loss through thermal and nonthermal processes, and an accurate accounting of the EUV energy deposition in a planet's energy budget is essential, especially for terrestrial habitable zone planets and close-in gaseous planets. Direct EUV observations of exoplanet host stars would require a new, dedicated observatory. Archival observations from the $\textit{EUVE}$ satellite, models, and theory alone are insufficient to accurately characterize EUV spectra of the majority of exoplanet host stars, especially for low-mass stars.

The Interface Region Imaging Spectrograph (IRIS) routinely observes the Si IV resonance lines. When analyzing observations of these lines it has typically been assumed they form under optically thin conditions. This is likely valid for the quiescent Sun, but this assumption has also been applied to the more extreme flaring scenario. We used 36 electron beam driven radiation hydrodynamic solar flare simulations, computed using the RADYN code, to probe the validity of this assumption. Using these simulated atmospheres we solved the radiation transfer equations to obtain the non-LTE, non-equilibrium populations, line profiles, and opacities for a model Silicon atom, including charge exchange processes. This was achieved using the `minority species' version of RADYN. The inclusion of charge exchange resulted in a substantial fraction of Si IV at cooler temperatures than those predicted by ionisation equilibrium. All simulations with an injected energy flux $F>5\times10^{10}$ erg cm$^{-2}$ s$^{-1}$ resulted in optical depth effects on the Si IV emission, with differences in both intensity and line shape compared to the optically thin calculation. Weaker flares (down to $F\approx5\times10^{9}$ erg cm$^{-2}$ s$^{-1}$) also resulted in Si IV emission forming under optically thick conditions, depending on the other beam parameters. When opacity was significant, the atmospheres generally had column masses in excess of $5\times10^{-6}$ g cm$^{-2}$ over the temperature range $40$ to $100$ kK, and the Si IV formation temperatures were between $30$ and $60$ kK. We urge caution when analyzing Si IV flare observations, or when computing synthetic emission without performing a full radiation transfer calculation.

The Extreme Ultraviolet Normal Incidence Spectrograph (EUNIS) and the Hinode/ EUV Imaging Spectrometer (EIS) observed AR 11726 on 2013 April 23. We present intensity images in numerous atomic lines constructed from these observations. These lines are formed over a wide range of temperatures, and we use their relative intensities to constrain a parameterization of nanoflare heating. We construct a 3D model of the magnetic field in this active region by extrapolating the surface magnetic field into the corona and using SDO/AIA images of coronal loops to ensure that extrapolated magnetic field lines co-align with observed coronal loops. We trace 2848 magnetic field lines within the volume of this active region and model how they fill with hot plasma in response to nanoflare heating. We perform a parameter study to determine how the frequency and energy released in nanoflares scale with magnetic field strength and loop length. From our 3D model, we construct synthetic images of the lines observed by EUNIS and EIS and constrain the parameter study by minimizing the difference between the synthetic and observed images.

We investigate the cause of the suppressed Balmer series and the origin of the white-light continuum emission in the X1.0 class solar flare on 2014 June 11. We use radiative hydrodynamic simulations to model the response of the flaring atmosphere to both electron and proton beams which are energetically constrained using RHESSI and Fermi observations. A comparison of synthetic spectra with the observations allow us to narrow the range of beam fluxes and low energy cut-off that may be applicable to this event. We conclude that the electron and proton beams that can reproduce the observed spectral features are those that have relatively low fluxes and high values for the low energy cut-off. While electron beams shift the upper chromosphere and transition region to greater geometrical heights, proton beams with a similar flux leave these areas of the atmosphere relatively undisturbed. It is easier for proton beams to penetrate to the deeper layers and not deposit their energy in the upper chromosphere where the Balmer lines are formed. The relatively weak particle beams that are applicable to this flare do not cause a significant shift of the $\tau=1$ surface and the observed excess WL emission is optically thin.

We present the results of 1D hydrodynamic simulations of coronal loops which are subject to nanoflares, caused by either in-situ thermal heating, or non-thermal electrons (NTE) beams. The synthesized intensity and Doppler shifts can be directly compared with IRIS and AIA observations of rapid variability in the transition region (TR) of coronal loops, associated with transient coronal heating. We find that NTE with high enough low-energy cutoff (E$_\textrm{C}$) deposit energy in the lower TR and chromosphere causing blueshifts (up to~$\sim$~20 km/s) in the \emph{IRIS} \siiv~lines, which thermal conduction cannot reproduce. The E$_\textrm{C}$ threshold value for the blueshifts depends on the total energy of the events ($\approx$~5 keV for 10$^{24}$ ergs, up to 15 keV for 10$^{25}$ ergs). The observed footpoint emission intensity and flows, combined with the simulations, can provide constraints on both the energy of the heating event and E$_\textrm{C}$. The response of the loop plasma to nanoflares depends crucially on the electron density: significant \siiv~intensity enhancements and flows are observed only for initially low-density loops ($<$~10$^{9}$~cm$^{-3}$). This provides a possible explanation of the relative scarcity of observations of significant moss variability. While the TR response to single heating episodes can be clearly observed, the predicted coronal emission (AIA 94\AA) for single strands is below current detectability, and can only be observed when several strands are heated closely in time. Finally, we show that the analysis of the IRIS \mgii~chromospheric lines can help further constrain the properties of the heating mechanisms.

The origin of the near-ultraviolet and optical continuum radiation in flares is critical for understanding particle acceleration and impulsive heating in stellar atmospheres. Radiative-hydrodynamic simulations in 1D have shown that high energy deposition rates from electron beams produce two flaring layers at T~10^4 K that develop in the chromosphere: a cooling condensation (downflowing compression) and heated non-moving (stationary) flare layers just below the condensation. These atmospheres reproduce several observed phenomena in flare spectra, such as the red wing asymmetry of the emission lines in solar flares and a small Balmer jump ratio in M dwarf flares. The high beam flux simulations are computationally expensive in 1D, and the (human) timescales for completing NLTE models with adaptive grids in 3D will likely be unwieldy for a time to come. We have developed a prescription for predicting the approximate evolved states, continuum optical depth, and the emergent continuum flux spectra of radiative-hydrodynamic model flare atmospheres. These approximate prescriptions are based on an important atmospheric parameter: the column mass (m_ref) at which hydrogen becomes nearly completely ionized at the depths that are approximately in steady state with the electron beam heating. Using this new modeling approach, we find that high energy flux density (>F11) electron beams are needed to reproduce the brightest observed continuum intensity in IRIS data of the 2014-Mar-29 X1 solar flare and that variation in m_ref from 0.001 to 0.02 g/cm2 reproduces most of the observed range of the optical continuum flux ratios at the peaks of M dwarf flares.

The broadening of the hydrogen lines during flares is thought to result from increased charge (electron, proton) density in the flare chromosphere. However, disagreements between theory and modeling prescriptions have precluded an accurate diagnostic of the degree of ionization and compression resulting from flare heating in the chromosphere. To resolve this issue, we have incorporated the unified theory of electric pressure broadening of the hydrogen lines into the non-LTE radiative transfer code RH. This broadening prescription produces a much more realistic spectrum of the quiescent, A0 star Vega compared to the analytic approximations used as a damping parameter in the Voigt profiles. We test recent radiative-hydrodynamic (RHD) simulations of the atmospheric response to high nonthermal electron beam fluxes with the new broadening prescription and find that the Balmer lines are over-broadened at the densest times in the simulations. Adding many simultaneously heated and cooling model loops as a "multithread" model improves the agreement with the observations. We revisit the three-component phenomenological flare model of the YZ CMi Megaflare using recent and new RHD models. The evolution of the broadening, line flux ratios, and continuum flux ratios are well-reproduced by a multithread model with high-flux nonthermal electron beam heating, an extended decay phase model, and a "hot spot" atmosphere heated by an ultrarelativistic electron beam with reasonable filling factors: 0.1%, 1%, and 0.1% of the visible stellar hemisphere, respectively. The new modeling motivates future work to understand the origin of the extended gradual phase emission.

We present unique NUV observations of a well-observed X-class flare from NOAA 12087 obtained at Ond\v{r}ejov Observatory. The flare shows a strong white-light continuum but no detectable emission in the higher Balmer and Lyman lines. RHESSI and Fermi observations indicate an extremely hard X-ray spectrum and gamma-ray emission. We use the RADYN radiative hydrodynamic code to perform two type of simulations. One where an energy of 3 x 10^11 erg/cm^2/s is deposited by an electron beam with a spectral index of ~3 and a second where the same energy is applied directly to the photosphere. The combination of observations and simulations allow us to conclude that the white-light emission and the suppression or complete lack of hydrogen emission lines is best explained by a model where the dominant energy deposition layer is located in the lower layers of the solar atmosphere rather than the chromosphere.

Recent observations from the Interface Region Imaging Spectrograph (IRIS) appear to show impulsive brightenings in high temperature lines, which when combined with simultaneous ground based observations in H$\alpha$, appear co-spatial to Ellerman Bombs (EBs). We use the RADYN 1-dimensional radiative transfer code in an attempt to try and reproduce the observed line profiles and simulate the atmospheric conditions of these events. Combined with the MULTI/RH line synthesis codes, we compute the H$\alpha$, Ca II 8542~\AA, and Mg II h \& k lines for these simulated events and compare them to previous observations. Our findings hint that the presence of superheated regions in the photosphere ($>$10,000 K) is not a plausible explanation for the production of EB signatures. While we are able to recreate EB-like line profiles in H$\alpha$, Ca II 8542~\AA, and Mg II h \& k, we cannot achieve agreement with all of these simultaneously.

Even in the absence of resolved flares, the corona is heated to several million degrees. However, despite its importance for the structure, dynamics, and evolution of the solar atmosphere, the origin of this heating remains poorly understood. Several observational and theoretical considerations suggest that the heating is driven by small, impulsive energy bursts which could be Parker-style "nanoflares" (Parker 1988) that arise via reconnection within the tangled and twisted coronal magnetic field. The classical "smoking gun" (Klimchuk 2009; Cargill et al. 2013) for impulsive heating is the direct detection of widespread hot plasma (T > 6 MK) with a low emission measure. In recent years there has been great progress in the development of Transition Edge Sensor (TES) X-ray microcalorimeters that make them more ideal for studying the Sun. When combined with grazing-incidence focusing optics, they provide direct spectroscopic imaging over a broad energy band (0.5 to 10 keV) combined with extremely impressive energy resolution in small pixels, as low as 0.7 eV (FWHM) at 1.5 keV (Lee 2015), and 1.56 eV (FWHM) at 6 keV (Smith 2012), two orders of magnitude better than the current best traditional solid state photon-counting spectrometers. Decisive observations of the hot plasma associated with nanoflare models of coronal heating can be provided by new solar microcalorimeters. These measurements will cover the most important part of the coronal spectrum for searching for the nanoflare-related hot plasma and will characterize how much nanoflares can heat the corona both in active regions and the quiet Sun. Finally, microcalorimeters will enable to study all of this as a function of time and space in each pixel simultaneously a capability never before available.

How impulsive magnetic energy release leads to solar eruptions and how those eruptions are energized and evolve are vital unsolved problems in Heliophysics. The standard model for solar eruptions summarizes our current understanding of these events. Magnetic energy in the corona is released through drastic restructuring of the magnetic field via reconnection. Electrons and ions are then accelerated by poorly understood processes. Theories include contracting loops, merging magnetic islands, stochastic acceleration, and turbulence at shocks, among others. Although this basic model is well established, the fundamental physics is poorly understood. HXR observations using grazing-incidence focusing optics can now probe all of the key regions of the standard model. These include two above-the-looptop (ALT) sources which bookend the reconnection region and are likely the sites of particle acceleration and direct heating. The science achievable by a direct HXR imaging instrument can be summarized by the following science questions and objectives which are some of the most outstanding issues in solar physics (1) How are particles accelerated at the Sun? (1a) Where are electrons accelerated and on what time scales? (1b) What fraction of electrons is accelerated out of the ambient medium? (2) How does magnetic energy release on the Sun lead to flares and eruptions? A Focusing Optics X-ray Solar Imager (FOXSI) instrument, which can be built now using proven technology and at modest cost, would enable revolutionary advancements in our understanding of impulsive magnetic energy release and particle acceleration, a process which is known to occur at the Sun but also throughout the Universe.

We study the temporal evolution of the Na I D1 line profiles in the M3.9 flare SOL2014-06-11T21:03 UT, using high spectral resolution observations obtained with the IBIS instrument on the Dunn Solar Telescope combined with radiative hydrodynamic simulations. Our results show a significant increase in line core and wing intensities during the flare. The analysis of the line profiles from the flare ribbons reveal that the Na I D1 line has a central reversal with excess emission in the blue wing (blue asymmetry). We combine RADYN and RH simulations to synthesise Na I D1 line profiles of the flaring atmosphere and find good agreement with the observations. Heating with a beam of electrons modifies the radiation field in the flaring atmosphere and excites electrons from the ground state $\mathrm{3s~^2S}$ to the first excited state $\mathrm{3p~^2P}$, which in turn modifies relative population of the two states. The change in temperature and the population density of the energy states make the sodium line profile revert from absorption into emission. Analysis of the simulated spectra also reveals that the Na I D1 flare profile asymmetries are produced by the velocity gradients generated %and opacity effects in the lower solar atmosphere.

The 2014 March 29 X1 solar flare (SOL20140329T17:48) produced bright continuum emission in the far- and near-ultraviolet (NUV) and highly asymmetric chromospheric emission lines, providing long-sought constraints on the heating mechanisms of the lower atmosphere in solar flares. We analyze the continuum and emission line data from the Interface Region Imaging Spectrograph (IRIS) of the brightest flaring magnetic footpoints in this flare. We compare the NUV spectra of the brightest pixels to new radiative-hydrodynamic predictions calculated with the RADYN code using constraints on a nonthermal electron beam inferred from the collisional thick-target modeling of hard X-ray data from RHESSI. We show that the atmospheric response to a high beam flux density satisfactorily achieves the observed continuum brightness in the NUV. The NUV continuum emission in this flare is consistent with hydrogen (Balmer) recombination radiation that originates from low optical depth in a dense chromospheric condensation and from the stationary beam-heated layers just below the condensation. A model producing two flaring regions (a condensation and stationary layers) in the lower atmosphere is also consistent with the asymmetric Fe II chromospheric emission line profiles observed in the impulsive phase.

Spectroscopic observations of solar flares provide critical diagnostics of the physical conditions in the flaring atmosphere. Some key features in observed spectra have not yet been accounted for in existing flare models. Here we report a data-driven simulation of the well-observed X1.0 flare on 2014 March 29 that can reconcile some well-known spectral discrepancies. We analyzed spectra of the flaring region from the Interface Region Imaging Spectrograph (IRIS) in MgII h&k, the Interferometric BIdimensional Spectropolarimeter at the Dunn Solar Telescope (DST/IBIS) in H$\alpha$ 6563 \AA\ and CaII 8542 \AA, and the Reuven Ramaty High Energy Solar Spectroscope Imager (RHESSI) in hard X-rays. We constructed a multi-threaded flare loop model and used the electron flux inferred from RHESSI data as the input to the radiative hydrodynamic code RADYN to simulate the atmospheric response. We then synthesized various chromospheric emission lines and compared them with the IRIS and IBIS observations. In general, the synthetic intensities agree with the observed ones, especially near the northern footpoint of the flare. The simulated MgII line profile has narrower wings than the observed one. This discrepancy can be reduced by using a higher microturbulent velocity (27 km/s) in a narrow chromospheric layer. In addition, we found that an increase of electron density in the upper chromosphere within a narrow height range of $\approx$800 km below the transition region can turn the simulated MgII line core into emission and thus reproduce the single peaked profile, which is a common feature in all IRIS flares.

We use radiation hydrodynamic simulations to examine two models of solar flare chromospheric heating: Alfv\'en wave dissipation and electron beam collisional losses. Both mechanisms are capable of strong chromospheric heating, and we show that the distinctive atmospheric evolution in the mid-to-upper chromosphere results in Mg II k-line emission that should be observably different between wave-heated and beam-heated simulations. We also present Ca II 8542A profiles which are formed slightly deeper in the chromosphere. The Mg II k-line profiles from our wave-heated simulation are quite different from those from a beam-heated model and are more consistent with IRIS observations. The predicted differences between the Ca II 8542A in the two models are small. We conclude that careful observational and theoretical study of lines formed in the mid-to-upper chromosphere holds genuine promise for distinguishing between competing models for chromospheric heating in flares.

The asymmetries observed in the line profiles of solar flares can provide important diagnostics of the properties and dynamics of the flaring atmosphere. In this paper the evolution of the Halpha and Ca II 8542 {\AA} lines are studied using high spatial, temporal and spectral resolution ground-based observations of an M1.1 flare obtained with the Swedish 1-m Solar Telescope. The temporal evolution of the Halpha line profiles from the flare kernel shows excess emission in the red wing (red asymmetry) before flare maximum, and excess in the blue wing (blue asymmetry) after maximum. However, the Ca II 8542 {\AA} line does not follow the same pattern, showing only a weak red asymmetry during the flare. RADYN simulations are used to synthesise spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. We show that the red asymmetry observed in Halpha is not necessarily associated with plasma downflows, and the blue asymmetry may not be related to plasma upflows. Indeed, we conclude that the steep velocity gradients in the flaring chromosphere modifies the wavelength of the central reversal in the Halpha line profile. The shift in the wavelength of maximum opacity to shorter and longer wavelengths generates the red and blue asymmetries, respectively.

We present a unified computational framework which can be used to describe impulsive flares on the Sun and on dMe stars. The models assume that the flare impulsive phase is caused by a beam of charged particles that is accelerated in the corona and propagates downward depositing energy and momentum along the way. This rapidly heats the lower stellar atmosphere causing it to explosively expand and dramatically brighten. Our models consist of flux tubes that extend from the sub-photosphere into the corona. We simulate how flare-accelerated charged particles propagate down one-dimensional flux tubes and heat the stellar atmosphere using the Fokker-Planck kinetic theory. Detailed radiative transfer is included so that model predictions can be directly compared with observations. The flux of flare-accelerated particles drives return currents which additionally heat the stellar atmosphere. These effects are also included in our models. We examine the impact of the flare-accelerated particle beams on model solar and dMe stellar atmospheres and perform parameter studies varying the injected particle energy spectra. We find the atmospheric response is strongly dependent on the accelerated particle cutoff energy and spectral index.

(abridged) The heating mechanism at high densities during M dwarf flares is poorly understood. Spectra of M dwarf flares in the optical and near-ultraviolet wavelength regimes have revealed three continuum components during the impulsive phase: 1) an energetically dominant blackbody component with a color temperature of T $\sim$ 10,000 K in the blue-optical, 2) a smaller amount of Balmer continuum emission in the near-ultraviolet at lambda $<$ 3646 Angstroms and 3) an apparent pseudo-continuum of blended high-order Balmer lines. These properties are not reproduced by models that employ a typical "solar-type" flare heating level in nonthermal electrons, and therefore our understanding of these spectra is limited to a phenomenological interpretation. We present a new 1D radiative-hydrodynamic model of an M dwarf flare from precipitating nonthermal electrons with a large energy flux of $10^{13}$ erg cm$^{-2}$ s$^{-1}$. The simulation produces bright continuum emission from a dense, hot chromospheric condensation. For the first time, the observed color temperature and Balmer jump ratio are produced self-consistently in a radiative-hydrodynamic flare model. We find that a T $\sim$ 10,000 K blackbody-like continuum component and a small Balmer jump ratio result from optically thick Balmer and Paschen recombination radiation, and thus the properties of the flux spectrum are caused by blue light escaping over a larger physical depth range compared to red and near-ultraviolet light. To model the near-ultraviolet pseudo-continuum previously attributed to overlapping Balmer lines, we include the extra Balmer continuum opacity from Landau-Zener transitions that result from merged, high order energy levels of hydrogen in a dense, partially ionized atmosphere. This reveals a new diagnostic of ambient charge density in the densest regions of the atmosphere that are heated during dMe and solar flares.

We investigated the response of the solar atmosphere to non-thermal electron beam heating using the radiative transfer and hydrodynamics modelling code RADYN. The temporal evolution of the parameters that describe the non-thermal electron energy distribution were derived from hard X-ray observations of a particular flare, and we compared the modelled and observed parameters. The evolution of the non-thermal electron beam parameters during the X1.5 solar flare on 2011 March 9 were obtained from analysis of RHESSI X-ray spectra. The RADYN flare model was allowed to evolve for 110 seconds, after which the electron beam heating was ended, and was then allowed to continue evolving for a further 300s. The modelled flare parameters were compared to the observed parameters determined from extreme-ultraviolet spectroscopy. The model produced a hotter and denser flare loop than that observed and also cooled more rapidly, suggesting that additional energy input in the decay phase of the flare is required. In the explosive evaporation phase a region of high-density cool material propagated upward through the corona. This material underwent a rapid increase in temperature as it was unable to radiate away all of the energy deposited across it by the non-thermal electron beam and via thermal conduction. A narrow and high-density ($n_{e} \le 10^{15}$ cm$^{-3}$) region at the base of the flare transition region was the source of optical line emission in the model atmosphere. The collision-stopping depth of electrons was calculated throughout the evolution of the flare, and it was found that the compression of the lower atmosphere may permit electrons to penetrate farther into a flaring atmosphere compared to a quiet Sun atmosphere.

The physical processes causing energy exchange between the Sun's hot corona and its cool lower atmosphere remain poorly understood. The chromosphere and transition region (TR) form an interface region between the surface and the corona that is highly sensitive to the coronal heating mechanism. High resolution observations with the Interface Region Imaging Spectrograph (IRIS) reveal rapid variability (about 20 to 60 seconds) of intensity and velocity on small spatial scales at the footpoints of hot dynamic coronal loops. The observations are consistent with numerical simulations of heating by beams of non-thermal electrons, which are generated in small impulsive heating events called "coronal nanoflares". The accelerated electrons deposit a sizable fraction of their energy in the chromosphere and TR. Our analysis provides tight constraints on the properties of such electron beams and new diagnostics for their presence in the nonflaring corona.

Recent progress in obtaining high spatial resolution images of the solar corona in the extreme-ultraviolet (EUV) with Hinode, TRACE, SDO and recent Hi-C missions and soft X-ray (SXR) bands opened a new avenue in understanding the solar coronal heating, the major goal of solar physics. The data from EUV/SXR missions suggest that solar corona is a non-uniform environment structured into active regions (AR) represented by bundles magnetic loops heated to temperatures exceeding 5 MK. Any viable coronal heating model should be capable of reproducing EUV and SXR emission from coronal active regions well as dynamic activity. Measurements of emission measures (EM) for ARs provide clues to time dependence of the heating mechanism: static versus impulsive. While static equilibrium coronal loop models are successful in reproducing SXR emission within an AR, they cannot adequately predict the bright EUV loops. Meantime, impulsive heating is capable in reproducing both EUV and SXR loop emission. The major goal of this paper is to construct realistic synthetic EM images of specific solar corona active region, AR 11117 by using our 1D fully non-linear time-dependent single-fluid hydrodynamic code. We first construct a magnetic skeleton for the entire active region using the HMI/SDO magnetogram for AR 11117 and populate magnetic field lines with plasma. We then parametrically specify impulsive heating of individual strands (flux tubes) comprising coronal loops. Next, we simulated the response of the entire active region (with LOS projection effects) to the heating function (volumetric heating rate) scaled with magnetic field and spatial scale parameters and find the best match between synthetic and actual (reconstructed) DEMs obtained by SDO.

This paper presents measurements of the energy radiated by the lower solar atmosphere, at optical, UV, and EUV wavelengths, during an X-class solar flare (SOL2011-02-15T01:56) in response to an injection of energy assumed to be in the form of nonthermal electrons. Hard X-ray observations from RHESSI were used to track the evolution of the parameters of the nonthermal electron distribution to reveal the total power contained in flare accelerated electrons. By integrating over the duration of the impulsive phase, the total energy contained in the nonthermal electrons was found to be $>2\times10^{31}$ erg. The response of the lower solar atmosphere was measured in the free-bound EUV continua of H I (Lyman), He I, and He II, plus the emission lines of He II at 304\AA\ and H I (Ly$\alpha$) at 1216\AA\ by SDO/EVE, the UV continua at 1600\AA\ and 1700\AA\ by SDO/AIA, and the WL continuum at 4504\AA, 5550\AA, and 6684\AA, along with the Ca II H line at 3968\AA\ using Hinode/SOT. The summed energy detected by these instruments amounted to $\sim3\times10^{30}$ erg; about 15% of the total nonthermal energy. The Ly$\alpha$ line was found to dominate the measured radiative losses. Parameters of both the driving electron distribution and the resulting chromospheric response are presented in detail to encourage the numerical modelling of flare heating for this event, to determine the depth of the solar atmosphere at which these line and continuum processes originate, and the mechanism(s) responsible for their generation.

Near-ultraviolet spectroscopic data obtained with the HST STIS instrument on the dMe flare star YZ Canis Minoris (YZ CMi) were analyzed. Flare and quiet intervals were identified from the broadband near-UV light curve, and the spectrum of each flare was separately extracted. Mg II and Fe II line profiles show similar behavior during the flares. Two large flares allowed time-resolved spectra to be analyzed, revealing a very broad component to the Mg II k line profile in at least one flare spectrum (F9b). If interpreted as a velocity, this component requires chromospheric material to be moving with FWHM ~ 250 km/sec, implying kinetic energy far in excess of the radiative energy. The Mg II k flare line profiles were compared to recent radiative hydrodynamic models of flare atmospheres undergoing electron beam heating. The models successfully predict red enhancements in the line profile with typical velocity of a few km/sec, but do not reproduce the flares showing blue enhancements, or the strongly broadened line observed in flare F9b. A more complete calculation of redistribution into the line wings, including the effect of collisions with the electron beam, may resolve the origin of the excess line broadening.

We report on multi-wavelength observations spanning radio to X-ray wavelengths of the M dwarf flare star, EV Lacertae, probing the characteristics of the outer atmospheric plasma from the upper chromosphere to the corona. We detect the star at a wavelength of 2 cm (15 GHz) for the first time. UV and FUV line profiles show evidence of nonthermal broadening, and the velocity width appear to peak at lower temperatures than in the Sun; this trend is confirmed in another active M dwarf flare star. Electron density measurements indicate nearly constant electron pressures between $\log T=$5.2 and 6.4. At higher coronal temperatures, there is a sharp increase of two orders of magnitude in density (n$_{e}\sim10^{13}$ cm$^{-3}$ at $\log T=$6.9). X-ray, EUV, FUV and NUV spectra constrain the DEM from the upper chromosphere through the corona. The coronal pressures are inconsistent with the assumption of hydrostatic equilibrium, either through EM modeling or application of scaling laws, and imply large conductive loss rates and a large energy input at the highest temperatures. The timescales for radiative and conductive losses in EV Lac's upper atmosphere imply that significant continued heating must occur for the corona to maintain its quiescent properties. The high frequency radio detection requires the high temperature X-ray-emitting coronal plasma to be spatially distinct from the radio emission source. Length scales in the low-temperature corona are markedly larger than those in the high-temperature corona, further suggestions of an inhomogeneous mixture of thermal and nonthermal coronal plasma.

We report on radiative hydrodynamic simulations of M dwarf stellar flares and compare the model predictions to observations of several flares. The flares were simulated by calculating the hydrodynamic response of a model M dwarf atmosphere to a beam of non-thermal electrons. Radiative backwarming through numerous soft X-ray, extreme ultraviolet, and ultraviolet transitions are also included. The equations of radiative transfer and statistical equilibrium are treated in non-LTE for many transitions of hydrogen, helium and the Ca II ion allowing the calculation of detailed line profiles and continuum radiation. Two simulations were carried out, with electron beam fluxes corresponding to moderate and strong beam heating. In both cases we find the dynamics can be naturally divided into two phases: an initial gentle phase in which hydrogen and helium radiate away much of the beam energy, and an explosive phase characterized by large hydrodynamic waves. During the initial phase, lower chromospheric material is evaporated into higher regions of the atmosphere causing many lines and continua to brighten dramatically. The He II Lya line is especially enhanced, becoming the brightest line in the flaring spectrum. The hydrogen Balmer lines also become much brighter and show very broad line widths, in agreement with observations. We compare our predicted Balmer decrements to decrements calculated for several flare observations and find the predictions to be in general agreement with the observations. During the explosive phase both condensation and evaporation waves are produced. The moderate flare simulation predicts a peak evaporation wave of 130 km/s and a condensation wave of 30 km/s. The velocity of the condensation wave matches velocities observed in several transition region lines.

We present optical spectra of a flare on Barnard's star. Several photospheric as well as chromospheric species were enhanced by the flare heating. An analysis of the Balmer lines shows that their shapes are best explained by Stark broadening rather than chromospheric mass motions. We estimate the temperature of the flaring region in the lower atmosphere to be >8000 K and the electron density to be ~10^14 cm^{-3}, similar to values observed in other dM flares. Because Barnard's star is considered to be one of our oldest neighbors, a flare of this magnitude is probably quite rare.

We report on radiative hydrodynamic simulations of moderate and strong solar flares. The flares were simulated by calculating the atmospheric response to a beam of non-thermal electrons injected at the apex of a one-dimensional closed coronal loop, and include heating from thermal soft X-ray, extreme ultraviolet and ultraviolet (XEUV) emission. The equations of radiative transfer and statistical equilibrium were treated in non-LTE and solved for numerous transitions of hydrogen, helium, and Ca II allowing the calculation of detailed line profiles and continuum emission. This work improves upon previous simulations by incorporating more realistic non-thermal electron beam models and includes a more rigorous model of thermal XEUV heating. We find XEUV backwarming contributes less than 10% of the heating, even in strong flares. The simulations show elevated coronal and transition region densities resulting in dramatic increases in line and continuum emission in both the UV and optical regions. The optical continuum reaches a peak increase of several percent which is consistent with enhancements observed in solar white light flares. For a moderate flare (~M-class), the dynamics are characterized by a long gentle phase of near balance between flare heating and radiative cooling, followed by an explosive phase with beam heating dominating over cooling and characterized by strong hydrodynamic waves. For a strong flare (~X-class), the gentle phase is much shorter, and we speculate that for even stronger flares the gentle phase may be essentially non-existent. During the explosive phase, synthetic profiles for lines formed in the upper chromosphere and transition region show blue shifts corresponding to a plasma velocity of ~120 km/s, and lines formed in the lower chromosphere show red shifts of ~40 km/s.

We present the results of a campaign to observe flares on the M dwarf flare star EV Lacertae over the course of two days in 2001 September, utilizing a combination of radio continuum, optical photometric and spectroscopic, ultraviolet spectroscopic, and X-ray spectroscopic observations, to characterize the multi-wavelength nature of flares from this active, single late-type star. We find flares in every wavelength region in which we observed. In the multi-wavelength context, the start of the intense radio flare is coincident with an impulsive optical U-band flare, to within one minute, and yet there is no signature of an X-ray response. There are other intervals of time where optical flaring and UV flaring is occurring, but these cannot be related to the contemporaneous X-ray flaring: the time-integrated luminosities do not match the instantaneous X-ray flare luminosity, as one would expect for the Neupert effect. We investigate the probability of chance occurrences of flares from disparate wavelength regions producing temporal coincidences, but find that not all the flare associations can be explained by a superposition of flares due to a high flaring rate. We caution against making causal associations of multi-wavelength flares based solely on temporal correlations for high flaring rate stars like EV Lac.