Search

Between 2017 and 2024, the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory has observed several white-light solar flares. Notably, for the X9.3 flare of September 6, 2017, HMI spectro-polarimetric observations reveal one or more locations within the umbra of the associated active region where the Fe I 6173 {\AA} line goes into full emission, indicating significant heating of the photosphere and lower chromosphere. For these flares, we performed a spectro-polarimetric analysis at the aforementioned locations using HMI 90s cadence Stokes data. At the Fe I emission locations, line-core emission is observed to last for a single 90 s frame and is either concurrent with or followed by increases in the line continuum intensity lasting 90 to 180 seconds. This is followed by a smooth decay to pre-flare conditions over the next three to twenty minutes. For most locations, permanent changes to the Stokes Q, U, and/or V profiles were observed, indicating long-lasting non-transient changes to the photospheric magnetic field. These emissions coincided with local maxima in hard X-ray emission observed by the Konus instrument onboard the Wind spacecraft, as well as local maxima in the time derivative of soft X-ray emission observed by GOES satellites. Comparison of the Fe I 6173 {\AA} line profile synthesis for the ad-hoc heating of the initial empirical VAL-S umbra model and quiescent Sun (VAL-C-like) model indicates that the Fe I 6173 {\AA} line emission in the white-light flare kernels could be explained by the strong heating of initially cool photospheric regions.

Solar flare ribbon fronts appear ahead of the bright structures that normally characterise solar flares, and can persist for an extended period of time in spatially localised patches before transitioning to `regular' bright ribbons. They likely represent the initial onset of flare energy deposition into the chromosphere. Chromospheric spectra (e.g. He I 10830A and the Mg II near-UV lines) from ribbon fronts exhibit properties rather different to typical flare behaviour. In prior numerical modelling efforts we were unable to reproduce the long lifetime of ribbon fronts. Here we present a series of numerical experiments that are rather simple but which have important implications. We inject a very low flux of nonthermal electrons ($F = 5\times10^{8}$ erg s$^{-1}$ cm$^{-2}$) into the chromosphere for 100 s before ramping up to standard flare energy fluxes $(F = 10^{10-11}$ erg s$^{-1}$ cm$^{-2}$). Synthetic spectra not only sustained their ribbon front-like properties for significantly longer, in the case of harder nonthermal electron spectra the ribbon front behaviour persisted for the entirety of this weak-heating phase. Lengthening or shortening the duration of the weak-heating phase commensurately lengthened or shortened the ribbon front lifetimes. Ribbon fronts transitioned to regular bright ribbons when the upper chromosphere became sufficiently hot and dense, which happened faster for softer nonthermal electron spectra. Thus, the lifetime of flare ribbon fronts are a direct measure of the duration over which a relatively low flux of high energy electrons precipitates to the chromosphere prior to the bombardment of a much larger energy flux., Comment: Accepted in the Astrophysical Journal

Optical flares have been observed from magnetically active stars for many decades; unsurprisingly, the spectra and temporal evolution are complicated. For example, the shortcomings of optically thin, static slab models have long been recognized when confronted with the observations. A less incorrect -- but equally simple -- phenomenological $T \approx 9000$ K blackbody model has instead been widely adopted in the absence of realistic (i.e., observationally-tested) time-dependent, atmospheric models that are readily available. We use the RADYN code to calculate a grid of 1D radiative-hydrodynamic stellar flare models that are driven by short pulses of electron-beam heating. The flare heating rates in the low atmosphere vary over many orders of magnitude in the grid, and we show that the models with high-energy electron beams compare well to the global trends in flux ratios from impulsive-phase stellar flare, optical spectra. The models also match detailed spectral line shape properties. We find that the pressure broadening and optical depths account for the broad components of the hydrogen Balmer $\gamma$ lines in a powerful flare with echelle spectra. The self-consistent formation of the wings and nearby continuum level provide insight into how high-energy electron beam heating evolves from the impulsive to the gradual decay phase in white-light stellar flares. The grid is publicly available, and we discuss possible applications., Comment: 25 pages, 9 figures, accepted for publication in the Astrophysical Journal (v2 matches published version)

The bulk of solar flare energy is deposited in the chromosphere. Flare ribbons and footpoints in the chromosphere therefore offer great diagnostic potential of flare energy release and transport processes. High quality observations from the IRIS spacecraft have transformed our view of the Sun's atmospheric response to flares. Since most of the chromospheric lines observed by IRIS are optically thick, forward modelling is required to fully appreciate and extract the information they carry. Reproducing certain aspects of the Mg II lines remain frustratingly out of reach in state-of-the-art flare models, which are unable to satisfactorily reproduce the very broad line profiles. A commonly proposed resolution to this is to assert that very large values of `microturbulence' is present. We asses the validity of that approach by analysing optically thin lines in the flare chromosphere from the X-class flare SOL2014-10-25T17:08:00, using the derived value of nonthermal width as a constraint to our numerical models. A nonthermal width of the order 10~km~s$^{-1}$ was found within the short-lived red wing components of three spectral lines, with relatively narrow stationary components. Simulations of this flare were produced, and in the post-processing spectral synthesis we include within the downflows a microturbulence of 10~km~s$^{-1}$. While we can reproduce the O I 1355.598~\AA\ line rather well, and we can capture the general shape and properties of the Mg II line widths, the synthetic lines are still too narrow., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Solar flares and coronal mass ejections are interrelated phenomena that together are known as solar eruptive events. These are the main drivers of space weather and understanding their origins is a primary goal of Heliophysics. In this white paper, we advocate for the allocation of sufficient resources to bring together experts in observations and modeling to construct and test next generation data-driven models of solar eruptive events. We identify the key components necessary for constructing comprehensive end-to-end models including global scale 3D MHD resolving magnetic field evolution and reconnection, small scale simulations of particle acceleration in reconnection exhausts, kinetic scale transport of flare-accelerated particles into the lower solar atmosphere, and the radiative and hydrodynamics responses of the solar atmosphere to flare heating. Using this modeling framework, long-standing questions regarding how solar eruptive events release energy, accelerate particles, and heat plasma can be explored. To address open questions in solar flare physics, we recommend that NASA and NSF provide sufficient research and analysis funds to bring together a large body of researchers and numerical tools to tackle the end-to-end modeling framework that we outline. Current dedicated theory and modeling funding programs are relatively small scale and infrequent; funding agencies must recognize that modern space physics demands the use of both observations and modeling to make rapid progress., Comment: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 9 pages, 4 figures

FOXSI is a direct-imaging, hard X-ray (HXR) telescope optimized for solar flare observations. It detects hot plasma and energetic electrons in and near energy release sites in the solar corona via bremsstrahlung emission, measuring both spatial structure and particle energy distributions. It provides two orders of magnitude faster imaging spectroscopy than previously available, probing physically relevant timescales (<1s) never before accessible to address fundamental questions of energy release and efficient particle acceleration that have importance far beyond their solar application (e.g., planetary magnetospheres, flaring stars, accretion disks). FOXSI measures not only the bright chromospheric X-ray emission where electrons lose most of their energy, but also simultaneous emission from electrons as they are accelerated in the corona and propagate along magnetic field lines. FOXSI detects emission from high in the tenuous corona, where previous instruments have been blinded by nearby bright features and will fully characterizes the accelerated electrons and hottest plasmas as they evolve in energy, space, and time to solve the mystery of how impulsive energy release leads to solar eruptions, the primary drivers of space weather at Earth, and how those eruptions are energized and evolve., Comment: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 14 pages, 4 figures, 1 table

SDO/HMI observations reveal a class of solar flares with substantial energy and momentum impacts in the photosphere, concurrent with white-light emission and helioseismic responses, known as sunquakes. Previous radiative hydrodynamic modeling has demonstrated the challenges of explaining sunquakes in the framework of the standard flare model of `electron beam' heating. One of the possibilities to explain the sunquakes and other signatures of the photospheric impact is to consider additional heating mechanisms involved in solar flares, for example, via flare-accelerated protons. In this work, we analyze a set of single-loop Fokker-Planck and radiative hydrodynamics RADYN+FP simulations where the atmosphere is heated by non-thermal power-law-distributed proton beams which can penetrate deeper than the electron beams into the low atmospheric layers. Using the output of the RADYN models, we calculate synthetic Fe I 6173 A line Stokes profiles and from those the line-of-sight (LOS) observables of the SDO/HMI instrument, as well as the 3D helioseismic response and compare them with the corresponding observational characteristics. These initial results show that the models with proton beam heating can produce the enhancement of the HMI continuum observable and explain qualitatively generation of sunquakes. The continuum observable enhancement is evident in all models but is more prominent in ones with $E_{c}\geq$500 keV. In contrast, the models with $E_{c}\leq$100 keV provide a stronger sunquake-like helioseismic impact according to the 3D acoustic modeling, suggesting that low-energy (deka- and hecto-keV) protons have an important role in the generation of sunquakes., Comment: 22 pages, 11 figures, 2 tables

It is essential that there be coordinated and co-optimized observations in X-rays, gamma-rays, and EUV during the peak of solar cycle 26 (~2036) to significantly advance our understanding of impulsive energy release in the corona. The open questions include: What are the physical origins of space-weather events? How are particles accelerated at the Sun? How is impulsively released energy transported throughout the solar atmosphere? How is the solar corona heated? Many of the processes involved in triggering, driving, and sustaining solar eruptive events -- including magnetic reconnection, particle acceleration, plasma heating, and energy transport in magnetized plasmas -- also play important roles in phenomena throughout the Universe. This set of observations can be achieved through a single flagship mission or, with foreplanning, through a combination of major missions (e.g., the previously proposed FIERCE mission concept)., Comment: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 5 pages, 1 figure

Understanding the nature of energetic particles in the solar atmosphere is one of the most important outstanding problems in heliophysics. Flare-accelerated particles compose a huge fraction of the flare energy budget; they have large influences on how events develop; they are an important source of high-energy particles found in the heliosphere; and they are the single most important corollary to other areas of high-energy astrophysics. Despite the importance of this area of study, this topic has in the past decade received only a small fraction of the resources necessary for a full investigation. For example, NASA has selected no new Explorer-class instrument in the past two decades that is capable of examining this topic. The advances that are currently being made in understanding flare-accelerated electrons are largely undertaken with data from EOVSA (NSF), STIX (ESA), and NuSTAR (NASA Astrophysics). This is despite the inclusion in the previous Heliophysics decadal survey of the FOXSI concept as part of the SEE2020 mission, and also despite NASA's having invested heavily in readying the technology for such an instrument via four flights of the FOXSI sounding rocket experiment. Due to that investment, the instrumentation stands ready to implement a hard X-ray mission to investigate flare-accelerated electrons. This white paper describes the scientific motivation for why this venture should be undertaken soon., Comment: White paper submitted to the Decadal Survey for Solar and Space Physics (Heliophysics) 2024-2033; 15 pages, 5 figures

Context: Solar flares are the result of the sudden release of magnetic energy in the corona. Much of this energy goes into accelerating charged particles to high velocity. These particles travel along the magnetic field and the energy is dissipated when the density gets high enough, primarily in the solar chromosphere. Modelling this region is difficult because the radiation energy balance is dominated by strong, optically thick spectral lines. Aims: Our aim is to provide the community with realistic simulations of a flaring loop with an emphasis on the detailed treatment of the chromospheric energy balance. This will enable a detailed comparison of existing and upcoming observations with synthetic observables from the simulations, thereby elucidating the complex interactions in a flaring chromosphere. Methods: We used the 1D radiation hydrodynamics code RADYN to perform simulations of the effect of a beam of electrons injected at the apex of a solar coronal loop. A grid of models was produced, varying the total energy input, the steepness, and low-energy cutoff of the beam energy spectrum. Results: The full simulation results for a grid of models are made available online. Some general properties of the simulations are discussed., Comment: 11 pages, 9 figures. Accepted for publication in Astronomy & Astrophysics

Solar flares are efficient particle accelerators, with a substantial fraction of the energy released manifesting as non-thermal particles. While the role that non-thermal electrons play in transporting flare energy is well studied, the properties and importance of non-thermal protons is rather less well understood. This is in large part due to the paucity of diagnostics, particularly at the lower-energy (deka-keV) range of non-thermal proton distributions in flares. One means to identify the presence of deka-keV protons is by an effect originally described by \cite{1976ApJ...208..618O}. In the Orrall-Zirker effect, non-thermal protons interact with ambient neutral hydrogen, and via charge exchange produce a population of energetic neutral atoms (ENAs) in the chromosphere. These ENAs subsequently produce an extremely redshifted photon in the red wings of hydrogen spectral lines. We revisit predictions of the strength of this effect using modern interaction cross-sections, and numerical models capable of self-consistently simulating the flaring non-equilibrium ionization stratification, and the non-thermal proton distribution (and, crucially, their feedback on each other). We synthesize both the thermal and non-thermal emission from \lya\ and \lyb, the most promising lines that may exhibit a detectable signal. These new predictions are are weaker and more transient than prior estimates, but the effects should be detectable in fortuitous circumstances. We degrade the \lyb\ emission to the resolution of the Spectral Imaging of the Coronal Environment (SPICE) instrument on board Solar Orbiter, demonstrating that though likely difficult, it should be possible to detect the presence of non-thermal protons in flares observed by SPICE., Comment: Accepted for publication in The Astrophysical Journal

During solar flares plasma is typically heated to very high temperatures, and the resulting redistribution of energy via thermal conduction is a primary mechanism transporting energy throughout the flaring solar atmosphere. The thermal flux is usually modeled using Spitzer's theory, which is based on local Coulomb collisions between the electrons carrying the thermal flux and those in the background. However, often during flares, temperature gradients become sufficiently steep that the collisional mean free path exceeds the temperature gradient scale size, so that thermal conduction becomes inherently non-local. Further, turbulent angular scattering, which is detectable in nonthermal widths of atomic emission lines, can also act to increase the collision frequency and so suppress the heat flux. Recent work by Emslie & Bian (2018) extended Spitzer's theory of thermal conduction to account for both non-locality and turbulent suppression. We have implemented their theoretical expression for the heat flux (which is a convolution of the Spitzer flux with a kernel function) into the RADYN flare-modeling code and performed a parameter study to understand how the resulting changes in thermal conduction affect flare dynamics and hence the radiation produced. We find that models with reduced heat fluxes predict slower bulk flows, less intense line emission, and longer cooling times. By comparing features of atomic emission lines predicted by the models with Doppler velocities and nonthermal line widths deduced from a particular flare observation, we find that models with suppression factors between 0.3 to 0.5 relative to the Spitzer value best reproduce observed Doppler velocities across emission lines forming over a wide range of temperatures. Interestingly, the model that best matches observed nonthermal line widths has a kappa-type velocity distribution function., Comment: Accepted by the Astrophysical Journal. 12 pages, 6 figures

Red-shifted components of chromospheric emission lines in the hard X-ray impulsive phase of solar flares have recently been studied through their 30 s evolution with the high resolution of IRIS. Radiative-hydrodynamic flare models show that these redshifts are generally reproduced by electron-beam generated chromospheric condensations. The models produce large ambient electron densities, and the pressure broadening of hydrogen Balmer series should be readily detected in observations. To accurately interpret upcoming spectral data of flares with the DKIST, we incorporate non-ideal, non-adiabatic line broadening profiles of hydrogen into the RADYN code. These improvements allow time-dependent predictions for the extreme Balmer line wing enhancements in solar flares. We study two chromospheric condensation models, which cover a range of electron beam fluxes ($1-5 \times 10^{11}$ erg s$^{-1}$ cm$^{-2}$) and ambient electron densities ($1 - 60 \times 10^{13}$ cm$^{-3}$) in the flare chromosphere. Both models produce broadening and redshift variations within 10 s of the onset of beam heating. In the chromospheric condensations, there is enhanced spectral broadening due to large optical depths at H$\alpha$, H$\beta$, and H$\gamma$, while the much lower optical depth of the Balmer series H12$-$H16 provides a translucent window into the smaller electron densities in the beam-heated layers below the condensation. The wavelength ranges of typical DKIST/ViSP spectra of solar flares will be sufficient to test the predictions of extreme hydrogen wing broadening and accurately constrain large densities in chromospheric condensations., Comment: 33 pages, 10 figures, accepted for publication in the Astrophysical Journal

Current state-of-the-art spectrographs cannot resolve the fundamental spatial (sub-arcseconds) and temporal scales (less than a few tens of seconds) of the coronal dynamics of solar flares and eruptive phenomena. The highest resolution coronal data to date are based on imaging, which is blind to many of the processes that drive coronal energetics and dynamics. As shown by IRIS for the low solar atmosphere, we need high-resolution spectroscopic measurements with simultaneous imaging to understand the dominant processes. In this paper: (1) we introduce the Multi-slit Solar Explorer (MUSE), a spaceborne observatory to fill this observational gap by providing high-cadence (<20 s), sub-arcsecond resolution spectroscopic rasters over an active region size of the solar transition region and corona; (2) using advanced numerical models, we demonstrate the unique diagnostic capabilities of MUSE for exploring solar coronal dynamics, and for constraining and discriminating models of solar flares and eruptions; (3) we discuss the key contributions MUSE would make in addressing the science objectives of the Next Generation Solar Physics Mission (NGSPM), and how MUSE, the high-throughput EUV Solar Telescope (EUVST) and the Daniel K Inouye Solar Telescope (and other ground-based observatories) can operate as a distributed implementation of the NGSPM. This is a companion paper to De Pontieu et al. (2021; arXiv:2106.15584), which focuses on investigating coronal heating with MUSE., Comment: 42 pages, 22 figures, submitted to ApJ

While solar flares are predominantly characterised by an intense broadband enhancement to the solar radiative output, certain spectral lines and continua will, in theory, exhibit flare-induced dimmings. Observations of transitions of orthohelium He I $\lambda\lambda 10830${\AA} and the He I D3 lines have shown evidence of such dimming, usually followed by enhanced emission. It has been suggested that non-thermal collisional ionisation of helium by an electron beam, followed by recombinations to orthohelium, is responsible for overpopulating the those levels, leading to stronger absorption. However it has not been possible observationally to preclude the possibility of overpopulating orthohelium via enhanced photoionisation of He I by EUV irradiance from the flaring corona followed by recombinations. Here we present radiation hydrodynamics simulations of non-thermal electron beam-driven flares where (1) both non-thermal collisional ionisation of Helium and coronal irradiance are included, and (2) only coronal irradiance is included. A grid of simulations covering a range of total energies deposited by the electron beam, and a range of non-thermal electron beam low-energy cutoff values, were simulated. In order to obtain flare-induced dimming of the He I 10830{\AA} line it was necessary for non-thermal collisional ionisations to be present. The effect was more prominent in flares with larger low-energy cutoff values and longer lived in weaker flares and flares with a more gradual energy deposition timescale. These results demonstrate the usefulness of orthohelium line emission as a diagnostic of flare energy transport., Comment: 18 pages, 17 figures, accepted for publication in The Astrophysical Journal

During solar flares, a large flux of energetic electrons propagate from the tops of reconnecting magnetic flux tubes toward the lower atmosphere. Over the course of the electrons' transport, a co-spatial counter-streaming return current is induced, thereby balancing the current density. In response to the return current electric field, a fraction of the ambient electrons will be accelerated into the runaway regime. However, models describing the accelerated electron beam/return-current system have generally failed to take these suprathermal runaway electrons into account self-consistently. We develop a model in which an accelerated electron beam drives a steady-state, sub-Dreicer co-spatial return-current electric field, which locally balances the direct beam current and freely accelerates a fraction of background (return-current) electrons. The model is self-consistent, i.e., the electric field induced by the co-evolution of the direct beam and the runaway current is considered. We find that (1) the return current electric field can return a significant number of suprathermal electrons to the acceleration region, where they can be further accelerated to higher energies, runaway electrons can be a few tens of percent of the return current flux returning to the nonthermal beam's acceleration region, (2) the energy gain of the suprathermal electrons can be up to $10-35$\,keV, (3) the heating rate in the corona can be reduced by an order of magnitude in comparison to models which neglect the runaway component. The results depend on the injected beam flux density, the temperature and density of the background plasma., Comment: 19 pages, 10 figures, submitted to ApJ

Solar flares are a fundamental component of solar eruptive events (SEEs; along with solar energetic particles, SEPs, and coronal mass ejections, CMEs). Flares are the first component of the SEE to impact our atmosphere, which can set the stage for the arrival of the associated SEPs and CME. Magnetic reconnection drives SEEs by restructuring the solar coronal magnetic field, liberating a tremendous amount of energy which is partitioned into various physical manifestations: particle acceleration, mass and magnetic-field eruption, atmospheric heating, and the subsequent emission of radiation as solar flares. To explain and ultimately predict these geoeffective events, the heliophysics community requires a comprehensive understanding of the processes that transform and distribute stored magnetic energy into other forms, including the broadband radiative enhancement that characterises flares. This white paper, submitted to the Heliophysics 2050 Workshop, discusses the flare gradual phase part of SEEs, setting out the questions that need addressing via a combination of theoretical, modelling, and observational research. In short, the flare gradual phase persists much longer than predicted so, by 2050, we must identify the characteristics of the significant energy deposition sustaining the gradual phase, and address the fundamental processes of turbulence and non-local heat flux., Comment: White paper submitted to the Heliophysics 2050 Workshop

Solar flares are a fundamental component of solar eruptive events (SEEs; along with solar energetic particles, SEPs, and coronal mass ejections, CMEs). Flares are the first component of the SEE to impact our atmosphere, which can set the stage for the arrival of the associated SEPs and CME. Magnetic reconnection drives SEEs by restructuring the solar coronal magnetic field, liberating a tremendous amount of energy which is partitioned into various physical manifestations: particle acceleration, mass and magnetic-field eruption, atmospheric heating, and the subsequent emission of radiation as solar flares. To explain and ultimately predict these geoeffective events, the heliophysics community requires a comprehensive understanding of the processes that transform and distribute stored magnetic energy into other forms, including the broadband radiative enhancement that characterises flares. This white paper, submitted to the Heliophysics 2050 Workshop, discusses the flare impulsive phase part of SEEs, setting out the questions that need addressing via a combination of theoretical, modelling, and observational research. In short, by 2050 we must determine the mechanisms of particle acceleration and propagation, and must push beyond the paradigm of energy transport via nonthermal electron beams, to also account for accelerated protons & ions and downward directed Alfven waves., Comment: White paper submitted to the Heliophysics 2050 Workshop

We describe a new approach for modeling the transport of high energy particles accelerated during flares from the acceleration region in the solar corona until their eventual thermalization in the flare footpoint. Our technique numerically solves the Fokker-Planck equation and includes forces corresponding to Coulomb collisions in a flux loop with nonuniform ionization, synchrotron emission reaction, magnetic mirroring and a return current electric field. Our solution to the Fokker-Planck equation includes second-order pitch angle and momentum diffusion. It is applicable to particles of arbitrary mass and charge. By tracking the collisions, we predict the bremsstrahlung produced as these particles interact with the ambient stellar atmosphere. This can be compared directly with observations and used to constrain the accelerated particle energy distribution. We have named our numerical code FP and have distributed it for general use. We demonstrate its effectiveness in several test cases., Comment: Accepted for publication in the Astrophysical Journal, 19 pages, 7 figures

Solar flares are 3D phenomenon but modelling a flare in 3D, including many of the important processes in the chromosphere, is a computational challenge. Accurately modelling the chromosphere is important, even if the transition region and corona are the areas of interest, due to the flow of energy, mass, and radiation through the interconnected layers. We present a solar flare arcade model, that aims to bridge the gap between 1D and 3D modelling. Our approach is limited to the synthesis of optically thin emission. Using observed active region loop structures in a 3D domain we graft simulated 1D flare atmospheres onto each loop, synthesise the emission and then project that emission onto to the 2D observational plane. Emission from SDO/AIA, GOES/XRS, and IRIS/SG Fe XXI 1354.1A was forward modelled. We analyse the temperatures, durations, mass flows, and line widths associated with the flare, finding qualitative agreement but certain quantitative differences. Compared to observations, the Doppler shifts are of similar magnitude but decay too quickly. They are not as ordered, containing a larger amount of scatter compared to observations. The duration of gradual phase emission from GOES and AIA emission is also too short. Fe XXI lines are broadened, but not sufficiently. These findings suggest that additional physics is required in our model. The arcade model that we show here as a proof-of-concept can be extended to investigate other lines and global aspects of solar flares, providing a means to better test the coronal response to models of flare energy injection., Comment: Accepted for publication in the Astrophysical Journal, 22 pages, 18 figures

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., Comment: 18 pages, 17 figures, Accepted for Publication in PASJ

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., Comment: Accepted for publication by the Astrophysical Journal - 9th April 2020

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., Comment: 12 pages, 13 figures, accepted in the Astrophysical Journal

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., Comment: Accepted for publication in the Astrophysical Journal

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., Comment: 19 pg, 10 fig, 2 appendices. Accepted for publication in ApJ

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., Comment: Science white paper submitted to the Astro2020 Decadal Survey. 8 pages, 2 figures

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., Comment: Science white paper submitted to the Astro2020 Decadal Survey. 7 pages, 2 figures

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., Comment: Accepted in the Astrophysical Journal, 23 pages (3 appendices, 15 figures)

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., Comment: This is a preliminary placeholder draft. It will certainly be replaced with an updated version

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., Comment: 19 pages, 11 figures

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., Comment: 29 pages, 9 figures, accepted for publication in the Astrophysical Journal

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., Comment: 31 pages, 13 figures, 2 tables, accepted for publication in the Astrophysical Journal

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., Comment: 13 pages

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., Comment: Next Generation Solar Physics Mission, 2 figures

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., Comment: Next Generation Solar Physics Mission white paper, 2 figures

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., Comment: 42 pages, 14 figures, 5 tables, accepted for publication in the Astrophysical Journal

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., Comment: 17 pages, 9 figures, submitted for peer-review

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., Comment: 14 pages, 18 figures, accepted in ApJ

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., Comment: Accepted for publication by the Astrophysical Journal

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., Comment: 12 pages, 12 figures

(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., Comment: 50 pages, 2 tables, 13 figures. Accepted for publication in the Solar Physics Topical Issue, "Solar and Stellar Flares". Version 2 (June 22, 2015): updated to include comments by Guest Editor. The final publication is available at Springer via http://dx.doi.org/10.1007/s11207-015-0708-x

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., Comment: Published in Science on October 17: http://www.sciencemag.org/content/346/6207/1255724 . 26 pages, 10 figures. Movies are available at: http://www.lmsal.com/~ptesta/iris_science_mov/

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., Comment: 14 pages, 18 figures. Accepted for publication in Astrophysics Journal

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., Comment: 29 pages, 12 figures; accepted to PASP

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., Comment: Accepted for publication in the ApJ