Modelling the solar Transition Region with the use of an Adaptive Conduction (TRAC) method permits fast and accurate numerical solutions of the field-aligned hydrodynamic equations, capturing the enthalpy exchange between the corona and transition region, when the corona undergoes impulsive heating. The TRAC method eliminates the need for highly resolved numerical grids in the transition region and the commensurate very short time steps that are required for numerical stability. When employed with coarse spatial resolutions, typically achieved in multi-dimensional magnetohydrodynamic codes, the errors at peak density are less than 5% and the computation time is three orders of magnitude faster than fully resolved field-aligned models. This paper presents further examples that demonstrate the versatility and robustness of the method over a range of heating events, including impulsive and quasi-steady footpoint heating. A detailed analytical assessment of the TRAC method is also presented, showing that the approach works through all phases of an impulsive heating event because (i) the total radiative losses and (ii) the total heating when integrated over the transition region are both preserved at all temperatures under the broadening modifications of the method. The results from the numerical simulations complement this conclusion.

Coronal rain is ubiquitous in flare loops, forming shortly after the onset of the solar flare. Rain is thought to be caused by a thermal instability, a localized runaway cooling of material in the corona. The models that demonstrate this require extremely long duration heating on the order of the radiative cooling time, localized near the footpoints of the loops. In flares, electron beams are thought to be the primary energy transport mechanism, driving strong footpoint heating during the impulsive phase that causes evaporation, filling and heating flare loops. Electron beams, however, do not act for a long period of time, and even supposing that they did, their heating would not remain localized at the footpoints. With a series of numerical experiments, we show directly that these two issues mean that electron beams are incapable of causing the formation of rain in flare loops. This result suggests that either there is another mechanism acting in flare loops responsible for rain, or that the modeling of the cooling of flare loops is somehow deficient. To adequately describe flares, the standard model must address this issue to account for the presence of coronal rain.

The relative amount of high temperature plasma has been found to be a useful diagnostic to determine the frequency of coronal heating on sub-resolution structures. When the loops are infrequently heated, a broad emission measure (EM) over a wider range of temperatures is expected. A narrower EM is expected for high frequency heating where the loops are closer to equilibrium. The soft X-ray spectrum contains many spectral lines that provide high temperature diagnostics, including lines from Fe XVII-XIX. This region of the solar spectrum will be observed by the Marshall Grazing Incidence Spectrometer (MaGIXS) in 2020. In this paper, we derive the expected spectral lines intensity in MaGIXS to varying amounts of high temperature plasma to demonstrate that a simple line ratio of these provides a powerful diagnostic to determine the heating frequency. Similarly, we examine ratios of AIA channel intensities, filter ratios from a XRT, and energy bands from the FOXSI sounding rocket to determine their sensitivity to this parameter. We find that both FOXSI and MaGIXS provide good diagnostic capability for high-temperature plasma. We then compare the predicted line ratios to the output of a numerical model and confirm the MaGIXS ratios provide an excellent diagnostic for heating frequency.

To adequately constrain the frequency of energy deposition in active region cores in the solar corona, systematic comparisons between detailed models and observational data are needed. In this paper, we describe a pipeline for forward modeling active region emission using magnetic field extrapolations and field-aligned hydrodynamic models. We use this pipeline to predict time-dependent emission from active region NOAA 1158 as observed by SDO/AIA for low-, intermediate-, and high-frequency nanoflares. In each pixel of our predicted multi-wavelength, time-dependent images, we compute two commonly-used diagnostics: the emission measure slope and the time lag. We find that signatures of the heating frequency persist in both of these diagnostics. In particular, our results show that the distribution of emission measure slopes narrows and the mean decreases with decreasing heating frequency and that the range of emission measure slopes is consistent with past observational and modeling work. Furthermore, we find that the time lag becomes increasingly spatially coherent with decreasing heating frequency while the distribution of time lags across the whole active region becomes more broad with increasing heating frequency. In a follow up paper, we train a random forest classifier on these predicted diagnostics and use this model to classify real AIA observations of NOAA 1158 in terms of the underlying heating frequency.

The solar atmosphere is dominated by loops of magnetic flux which connect the multi-million-degree corona to the much cooler chromosphere. The temperature and density structure of quasi-static loops is determined by the continuous flow of energy from the hot corona to the lower solar atmosphere. Loop scaling laws provide relationships between global properties of the loop (such as peak temperature, pressure, and length); they follow from the physical variable dependencies of various terms in the energy equation, and hence the form of the loop scaling law provides insight into the key physics that controls the loop structure. Traditionally, scaling laws have been derived under the assumption of collision-dominated thermal conduction. Here we examine the impact of different regimes of thermal conduction -- collision-dominated, turbulence-dominated, and free-streaming -- on the form of the scaling laws relating the loop temperature and heating rate to its pressure and half-length. We show that the scaling laws for turbulence-dominated conduction are fundamentally different than those for collision-dominated and free-streaming conduction, inasmuch as the form of the scaling laws now depend primarily on conditions at the low-temperature, rather than high-temperature, part of the loop. We also establish regimes in temperature and density space in which each of the applicable scaling laws prevail.

Thermal non-equilibrium (TNE) is believed to be a potentially important process in understanding some properties of the magnetically closed solar corona. Through one-dimensional hydrodynamic models, this paper addresses the importance of the numerical spatial resolution, footpoint heating timescales and background heating on TNE. Inadequate transition region (TR) resolution can lead to significant discrepancies in TNE cycle behaviour, with TNE being suppressed in under-resolved loops. A convergence on the periodicity and plasma properties associated with TNE required spatial resolutions of less than 2 km for a loop of length 180 Mm. These numerical problems can be resolved using an approximate method that models the TR as a discontinuity using a jump condition, as proposed by Johnston et al. (2017a,b). The resolution requirements (and so computational cost) are greatly reduced while retaining good agreement with fully resolved results. Using this approximate method we (i) identify different regimes for the response of coronal loops to time-dependent footpoint heating including one where TNE does not arise and (ii) demonstrate that TNE in a loop with footpoint heating is suppressed unless the background heating is sufficiently small. The implications for the generality of TNE are discussed.

The brightness of the emission from coronal loops in the solar atmosphere is strongly dependent on the temperature and density of the confined plasma. After a release of energy, these loops undergo a heating and upflow phase, followed by a cooling and downflow cycle. Throughout, there are significant variations in the properties of the coronal plasma. In particular, the increased coronal temperature leads to an excess downward heat flux that the transition region is unable to radiate. This generates an enthalpy flux from the transition region to the corona, increasing the coronal density. The enthalpy exchange is highly sensitive to the transition region resolution in numerical simulations. With a numerically under-resolved transition region, major errors occur in simulating the coronal density evolution and, thus, the predicted loop emission. This Letter presents a new method that addresses the difficulty of obtaining the correct interaction between the corona and corona/chromosphere interface. In the transition region, an adaptive thermal conduction approach is used that broadens any unresolved parts of the atmosphere. We show that this approach, referred to as TRAC, successfully removes the influence of numerical resolution on the coronal density response to heating while maintaining high levels of agreement with fully resolved models. When employed with coarse spatial resolutions, typically achieved in multi-dimensional MHD codes, the peak density errors are less than 3% and the computation time is three orders of magnitude faster than fully resolved field-aligned models. The advantages of using TRAC in field-aligned hydrodynamic and multi-dimensional magnetohydrodynamic simulations are discussed.

In the work described here we investigate atomic processes leading to the formation of emission lines within the IRIS wavelength range at temperatures near $10^5$~K. We focus on (1) non-equilibrium and (2) density-dependent effects influencing the formation and radiative properties of S IV and O IV. These two effects have significant impacts on spectroscopic diagnostic measurements of quantities associated with the plasma that emission lines from S IV and O IV provide. We demonstrate this by examining nanoflare-based coronal heating to determine what the detectable signatures are in transition region emission. A detailed comparison between predictions from numerical experiments and several sets of observational data is presented to show how one can ascertain when non-equilibrium ionization and/or density-dependent atomic processes are important for diagnosing nanoflare properties, the magnitude of their contribution, and what information can be reliably extracted from the spectral data. Our key findings are the following. (1) The S/O intensity ratio is a powerful diagnostic of non-equilibrium ionization. (2) Non-equilibrium ionization has a strong effect on the observed line intensities even in the case of relatively weak nanoflare heating. (3) The density-dependence of atomic rate coefficients is only important when the ion population is out of equilibrium. (4) In the sample of active regions we examined, weak nanoflares coupled with non-equilibrium ionization and density-dependent atomic rates were required to explain the observed properties (e.g. the S/O intensity ratios). (5) Enhanced S/O intensity ratios cannot be due solely to the heating strength and must depend on other processes (e.g. heating frequency, non-Maxwellian distributions).

Understanding the dynamics of the solar chromosphere is crucial to understanding energy transport across the solar atmosphere. The chromosphere is optically thick at many wavelengths and described by non-local thermodynamic equilibrium (NLTE), making it difficult to interpret observations. Furthermore, there is considerable evidence that the atmosphere is filamented, and that current instruments do not sufficiently resolve small scale features. In flares, it is likely that multithreaded models are required to describe the heating. The combination of NLTE effects and multithreaded modeling requires computationally demanding calculations, which has motivated the development of a model that can efficiently treat both. We describe the implementation of a solver in a hydrodynamic code for the hydrogen level populations that approximates the NLTE solutions. We derive an accurate electron density across the atmosphere, that includes the effects of non-equilibrium ionization for helium and metals. We show the effects on hydrodynamic simulations, which are used to synthesize light curves using a post-processing radiative transfer code. We demonstrate the utility of this model on IRIS observations of a small flare. We show that the Doppler shifts in Mg II, C II, and O I can be explained with a multithreaded model of loops subjected to electron beam heating, so long as NLTE effects are treated. The intensities, however, do not match observed values very well, which is due to assumptions about the initial atmosphere. We briefly show how altering the initial atmosphere can drastically alter line profiles, and therefore derived quantities, and suggest that it should be tuned to pre-flare observations.

Much evidence suggests that the solar corona is heated impulsively, meaning that nanoflares may be ubiquitous in quiet and active regions (ARs). Hard X-ray (HXR) observations with unprecedented sensitivity $>$3~keV are now enabled by focusing instruments. We analyzed data from the \textit{Focusing Optics X-ray Solar Imager (FOXSI)} rocket and the \textit{Nuclear Spectroscopic Telescope Array (NuSTAR)} spacecraft to constrain properties of AR nanoflares simulated by the EBTEL field-line-averaged hydrodynamics code. We generated model X-ray spectra by computing differential emission measures for homogeneous nanoflare sequences with heating amplitudes $H_0$, durations $\tau$, delay times between events $t_N$, and filling factors $f$. The single quiescent AR observed by \textit{FOXSI-2} on 2014 December 11 is well fit by nanoflare sequences with heating amplitudes 0.02 erg cm$^{-3}$ s$^{-1}$ $<$ $H_0$ $<$ 13 erg cm$^{-3}$ s$^{-1}$ and a wide range of delay times and durations. We exclude delays between events shorter than $\sim$900 s at the 90\% confidence level for this region. Three of five regions observed by {\nustar} on 2014 November 1 are well fit by homogeneous nanoflare models, while two regions with higher fluxes are not. Generally, the {\nustar} count spectra are well fit by nanoflare sequences with smaller heating amplitudes, shorter delays, and shorter durations than the allowed \textit{FOXSI-2} models. These apparent discrepancies are likely due to differences in spectral coverage between the two instruments and intrinsic differences among the regions. Steady heating ($t_N$ = $\tau$) was ruled out with $>$99\% confidence for all regions observed by either instrument.

Despite its prediction over two decades ago, the detection of faint, high-temperature ("hot") emission due to nanoflare heating in non-flaring active region cores has proved challenging. Using an efficient two-fluid hydrodynamic model, this paper investigates the properties of the emission expected from repeating nanoflares (a nanoflare train) of varying frequency as well as the separate heating of electrons and ions. If the emission measure distribution ($\mathrm{EM}(T)$) peaks at $T = T_m$, we find that $\mathrm{EM}(T_m)$ is independent of details of the nanoflare train, and $\mathrm{EM}(T)$ above and below $T_m$ reflects different aspects of the heating. Below $T_m$ the main influence is the relationship of the waiting time between successive nanoflares to the nanoflare energy. Above $T_m$ power-law nanoflare distributions lead to an extensive plasma population not present in a monoenergetic train. Furthermore, in some cases characteristic features are present in $\mathrm{EM}(T)$. Such details may be detectable given adequate spectral resolution and a good knowledge of the relevant atomic physics. In the absence of such resolution we propose some metrics that can be used to infer the presence of "hot" plasma.

The properties expected of "hot" non-flaring plasmas due to nanoflare heating in active regions are investigated using hydrodynamic modeling tools, including a two-fluid development of the EBTEL code. Here we study a single nanoflare and show that while simple models predict an emission measure distribution extending well above 10 MK that is consistent with cooling by thermal conduction, many other effects are likely to limit the existence and detectability of such plasmas. These include: differential heating between electrons and ions, ionization non-equilibrium and, for short nanoflares, the time taken for the coronal density to increase. The most useful temperature range to look for this plasma, often called the "smoking gun" of nanoflare heating, lies between $10^{6.6}$ and $10^7$ K. Signatures of the actual heating may be detectable in some instances.

The million degree plasma of the solar corona must be supplied by the underlying layers of the atmosphere. The mechanism and location of energy release, and the precise source of coronal plasma, remain unresolved. In earlier work we pursued the idea that warm plasma is supplied to the corona via direct heating of the chromosphere by nanoflares, contrary to the prevailing belief that the corona is heated in-situ and the chromosphere is subsequently energized and ablated by thermal conduction. We found that single (low-frequency) chromospheric nanoflares could not explain the observed intensities, Doppler-shifts, and red/blue asymmetries in Fe XII and XIV emission lines. In the present work we follow up on another suggestion that the corona could be powered by chromospheric nanoflares that repeat on a timescale substantially shorter than the cooling/draining timescale. That is, a single magnetic strand is re-supplied with coronal plasma before the existing plasma has time to cool and drain. We perform a series of hydrodynamic experiments and predict the Fe XII and XIV line intensities, Doppler-shifts, and red/blue asymmetries. We find that our predicted quantities disagree dramatically with observations and fully developed loop structures cannot be created by intermediate- or high-frequency chromospheric nanoflares. We conclude that the mechanism ultimately responsible for producing coronal plasma operates above the chromosphere, but this does not preclude the possibility of a similar mechanism powering the chromosphere; extreme examples of which may be responsible for heating chromospheric plasma to transition region temperatures (e.g. type II spicules).

In this work we investigate the global activity patterns predicted from a model active region heated by distributions of nanoflares that have a range of frequencies. What differs is the average frequency of the distributions. The activity patterns are manifested in time lag maps of narrow-band instrument channel pairs. We combine hydrodynamic and forward modeling codes with a magnetic field extrapolation to create a model active region and apply the time lag method to synthetic observations. Our aim is not to reproduce a particular set of observations in detail, but to recover some typical properties and patterns observed in active regions. Our key findings are the following. 1. cooling dominates the time lag signature and the time lags between the channel pairs are generally consistent with observed values. 2. shorter coronal loops in the core cool more quickly than longer loops at the periphery. 3. all channel pairs show zero time lag when the line-of-sight passes through coronal loop foot-points. 4. there is strong evidence that plasma must be re-energized on a time scale comparable to the cooling timescale to reproduce the observed coronal activity, but it is likely that a relatively broad spectrum of heating frequencies are operating across active regions. 5. due to their highly dynamic nature, we find nanoflare trains produce zero time lags along entire flux tubes in our model active region that are seen between the same channel pairs in observed active regions.

Observations of solar flares with RHESSI have shown X-ray sources traveling along flaring loops, from the corona down to the chromosphere and back up. The 28 November 2002 C1.1 flare, first observed with RHESSI by Sui et al. 2006 and quantitatively analyzed by O'Flannagain et al. 2013, very clearly shows this behavior. By employing numerical experiments, we use these observations of X-ray source height motions as a constraint to distinguish between heating due to a non-thermal electron beam and in situ energy deposition in the corona. We find that both heating scenarios can reproduce the observed light curves, but our results favor non-thermal heating. In situ heating is inconsistent with the observed X-ray source morphology and always gives a height dispersion with photon energy opposite to what is observed.

In the standard model of solar flares, energy deposition by a beam of electrons drives strong chromospheric evaporation leading to a significantly denser corona and much brighter emission across the spectrum. Chromospheric evaporation was examined in great detail by Fisher, Canfield, & McClymont (1985a,b,c), who described a distinction between two different regimes, termed explosive and gentle evaporation. In this work, we examine the importance of electron energy and stopping depths on the two regimes and on the atmospheric response. We find that with explosive evaporation, the atmospheric response does not depend strongly on electron energy. In the case of gentle evaporation, lower energy electrons are significantly more efficient at heating the atmosphere and driving up-flows sooner than higher energy electrons. We also find that the threshold between explosive and gentle evaporation is not fixed at a given beam energy flux, but also depends strongly on the electron energy and duration of heating. Further, at low electron energies, a much weaker beam flux is required to drive explosive evaporation.

We consider the lifetimes of spots on the Sun and other stars from the standpoint of magnetic diffusion. While normal magnetic diffusivity predicts lifetimes of sunspots that are too large by at least two orders of magnitude, turbulent magnetic diffusivity accounts for both the functional form of the solar empirical spot-lifetime relation and for the observed sunspot lifetimes, provided that the relevant diffusion length is the supergranule size. Applying this relation to other stars, the value of turbulent diffusivity depends almost entirely on supergranule size, with very weak dependence on other variables such as magnetic field strength and density. Overall, the best observational data for other stars is consistent with the extension of the solar relation provided that stellar supergranule sizes for some stars are significantly larger than they are on the Sun.

It has been suggested that the hot plasma of the solar corona comes primarily from impulsive heating events, or nanoflares, that occur in the lower atmosphere, either in the upper part of the ordinary chromosphere or at the tips of type II spicules. We test this idea with a series of hydrodynamic simulations. We find that synthetic Fe XII (195) and Fe XIV (274) line profiles generated from the simulations disagree dramatically with actual observations. The integrated line intensities are much too faint; the blue shifts are much too fast; the blue-red asymmetries are much too large; and the emission is confined to low altitudes. We conclude that chromospheric nanoflares are not a primary source of hot coronal plasma. Such events may play an important role in producing the chromosphere and powering its intense radiation, but they do not, in general, raise the temperature of the plasma to coronal values. Those cases where coronal temperatures are reached must be relatively uncommon. The observed profiles of Fe XII and Fe XIV come primarily from plasma that is heated in the corona itself, either by coronal nanoflares or a quasi-steady coronal heating process. Chromospheric nanoflares might play a role in generating waves that provide this coronal heating.

The collisional thick target model, wherein a large number of electrons are accelerated down a flaring loop, can be used to explain many observed properties of solar flares. In this study, we focus on the sensitivity of GOES flare classification to the properties of the thick target model. Using a hydrodynamic model with RHESSI-derived electron beam parameters, we explore the effects of the beam energy flux (or total non-thermal energy), the cut-off energy, and the spectral index of the electron distribution on the soft X-rays (SXRs) observed by GOES. We conclude that (1) the GOES class is proportional to the non-thermal energy $E^{\alpha}$ for $\alpha \approx 1.7$ in the low energy passband (1-8 \AA) and $\approx 1.6$ in the high energy passband (0.5-4 \AA); (2) the GOES class is only weakly dependent on the spectral index in both passbands; (3) increases in the cut-off will increase the flux in the 0.5-4 \AA\ passband, but decrease the flux in the 1-8 \AA\ passband, while decreases in the cut-off will cause a decrease in the 0.5-4 \AA\ passband and a slight increase in the 1-8 \AA\ passband.

We examine the radiative cooling of coronal loops and demonstrate that the recently identified catastrophic cooling (Reale and Landi, 2012) is due to the inability of a loop to sustain radiative / enthalpy cooling below a critical temperature, which can be > 1 MK in flares, 0.5 - 1 MK in active regions and 0.1 MK in long tenuous loops. Catastrophic cooling is characterised by a rapid fall in coronal temperature while the coronal density changes by a small amount. Analytic expressions for the critical temperature are derived and show good agreement with numerical results. This effect limits very considerably the lifetime of coronal plasmas below the critical temperature.

The effect of the numerical spatial resolution in models of the solar corona and corona / chromosphere interface is examined for impulsive heating over a range of magnitudes using one dimensional hydrodynamic simulations. It is demonstrated that the principle effect of inadequate resolution is on the coronal density. An underresolved loop typically has a peak density of at least a factor of two lower than a resolved loop subject to the same heating, with larger discrepencies in the decay phase. The temperature for under-resolved loops is also lower indicating that lack of resolution does not "bottle up" the heat flux in the corona. Energy is conserved in the models to under 1% in all cases, indicating that this is not responsible for the low density. Instead, we argue that in under-resolved loops the heat flux "jumps across" the transition region to the dense chromosphere from which it is radiated rather than heating and ablating transition region plasma. This emphasises the point that the interaction between corona and chromosphere occurs only through the medium of the transition region. Implications for three dimensional magnetohydrodynamic coronal models are discussed.

Most of our knowledge of the physical processes in distant plasmas is obtained through measurement of the radiation they produce. Here we provide an overview of the main collisional and radiative processes and examples of diagnostics relevant to the microphysical processes in the plasma. Many analyses assume a time-steady plasma with ion populations in equilibrium with the local temperature and Maxwellian distributions of particle velocities, but these assumptions are easily violated in many cases. We consider these departures from equilibrium and possible diagnostics in detail.

The time-dependence of heating in solar active regions can be studied by analyzing the slope of the emission measure distribution cool-ward of the peak. In a previous study we showed that low-frequency heating can account for 0% to 77% of active region core emission measures. We now turn our attention to heating by a finite succession of impulsive events for which the timescale between events on a single magnetic strand is shorter than the cooling timescale. We refer to this scenario as a "nanoflare train" and explore a parameter space of heating and coronal loop properties with a hydrodynamic model. Our conclusions are: (1) nanoflare trains are consistent with 86% to 100% of observed active region cores when uncertainties in the atomic data are properly accounted for; (2) steeper slopes are found for larger values of the ratio of the train duration $\Delta_H$ to the post-train cooling and draining timescale $\Delta_C$, where $\Delta_H$ depends on the number of heating events, the event duration and the time interval between successive events ($\tau_C$); (3) $\tau_C$ may be diagnosed from the width of the hot component of the emission measure provided that the temperature bins are much smaller than 0.1 dex; (4) the slope of the emission measure alone is not sufficient to provide information about any timescale associated with heating - the length and density of the heated structure must be measured for $\Delta_H$ to be uniquely extracted from the ratio $\Delta_H/\Delta_C$.

Observational measurements of active region emission measures contain clues to the time-dependence of the underlying heating mechanism. A strongly non-linear scaling of the emission measure with temperature indicates a large amount of hot plasma relative to warm plasma. A weakly non-linear (or linear) scaling of the emission measure indicates a relatively large amount of warm plasma, suggesting that the hot active region plasma is allowed to cool and so the heating is impulsive with a long repeat time. This case is called {\it low-frequency} nanoflare heating and we investigate its feasibility as an active region heating scenario here. We explore a parameter space of heating and coronal loop properties with a hydrodynamic model. For each model run, we calculate the slope $\alpha$ of the emission measure distribution $EM(T) \propto T^\alpha$. Our conclusions are: (1) low-frequency nanoflare heating is consistent with about 36% of observed active region cores when uncertainties in the atomic data are not accounted for; (2) proper consideration of uncertainties yields a range in which as many as 77% of observed active regions are consistent with low-frequency nanoflare heating and as few as zero; (3) low-frequency nanoflare heating cannot explain observed slopes greater than 3; (4) the upper limit to the volumetric energy release is in the region of 50 erg cm$^{-3}$ to avoid unphysical magnetic field strengths; (5) the heating timescale may be short for loops of total length less than 40 Mm to be consistent with the observed range of slopes; (6) predicted slopes are consistently steeper for longer loops.

A complete understanding of Doppler shift in active region loops can help probe the basic physical mechanism involved into the heating of those loops. Here we present observations of upflows in coronal loops detected in a range of temperature temperatures (log T=5.8 - 6.2). The loop was not discernible above these temperatures. The speed of upflow was strongest at the footpoint and decreased with height. The upflow speed at the footpoint was about 20 km/s in Fe VIII which decreased with temperature being about 13 km/s in Fe X, about 8 km/s in Fe XII and about 4 km/s in FeXIII. To the best of our knowledge this is the first observation providing evidence of upflow of plasma in coronal loop structures at these temperatures. We interpret these observations as evidence of chromospheric evaporation in quasi-static coronal loops.

This paper develops the zero-dimensional (0D) hydrodynamic coronal loop model "Enthalpy-based Thermal Evolution of Loops" (EBTEL) proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating, especially impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. The new features in EBTEL are important for accurate tracking of, in particular, the density. The 0D results are compared to a 1D hydro code (Hydrad) with generally good agreement. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function, (b) extensive parameter surveys and (c) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.

We present the science case for a broadband X-ray imager with high-resolution spectroscopy, including simulations of X-ray spectral diagnostics of both active regions and solar flares. This is part of a trilogy of white papers discussing science, instrument (Bandler et al. 2010), and missions (Bookbinder et al. 2010) to exploit major advances recently made in transition-edge sensor (TES) detector technology that enable resolution better than 2 eV in an array that can handle high count rates. Combined with a modest X-ray mirror, this instrument would combine arcsecondscale imaging with high-resolution spectra over a field of view sufficiently large for the study of active regions and flares, enabling a wide range of studies such as the detection of microheating in active regions, ion-resolved velocity flows, and the presence of non-thermal electrons in hot plasmas. It would also enable more direct comparisons between solar and stellar soft X-ray spectra, a waveband in which (unusually) we currently have much better stellar data than we do of the Sun.

In models of fast magnetic reconnection, flux transfer occurs within a small portion of a current sheet triggering stored magnetic energy to be thermalized by shocks. When the initial current sheet separates magnetic fields which are not perfectly anti-parallel, i.e. they are skewed, magnetic energy is first converted to bulk kinetic energy and then thermalized in slow magnetosonic shocks. We show that the latter resemble parallel shocks or hydrodynamic shocks for all skew angles except those very near the anti-parallel limit. As for parallel shocks, the structures of reconnection-driven slow shocks are best studied using two-fluid equations in which ions and electrons have independent temperature. Time-dependent solutions of these equations can be used to predict and understand the shocks from reconnection of skewed magnetic fields. The results differ from those found using a single-fluid model such as magnetohydrodynamics. In the two-fluid model electrons are heated indirectly and thus carry a heat flux always well below the free-streaming limit. The viscous stress of the ions is, however, typically near the fluid-treatable limit. We find that for a wide range of skew angles and small plasma beta an electron conduction front extends ahead of the slow shock but remains within the outflow jet. In such cases conduction will play a more limited role in driving chromospheric evaporation than has been predicted based on single-fluid, anti-parallel models.

Context: The location of coronal heating in magnetic loops has been the subject of a long-lasting controversy: does it occur mostly at the loop footpoints, at the top, is it random, or is the average profile uniform? Aims: We try to address this question in model loops with MHD turbulence and a profile of density and/or magnetic field along the loop. Methods: We use the ShellAtm MHD turbulent heating model described in Buchlin & Velli (2006), with a static mass density stratification obtained by the HydRad model (Bradshaw & Mason 2003). This assumes the absence of any flow or heat conduction subsequent to the dynamic heating. Results: The average profile of heating is quasi-uniform, unless there is an expansion of the flux tube (non-uniform axial magnetic field) or the variation of the kinetic and magnetic diffusion coefficients with temperature is taken into account: in the first case the heating is enhanced at footpoints, whereas in the second case it is enhanced where the dominant diffusion coefficient is enhanced. Conclusions: These simulations shed light on the consequences on heating profiles of the complex interactions between physical effects involved in a non-uniform turbulent coronal loop.