There has been tremendous progress in the degree of realism of three-dimensional radiation magneto-hydrodynamic simulations of the solar atmosphere in the past decades. Four of the most frequently used numerical codes are Bifrost, CO5BOLD, MANCHA3D, and MURaM. Here we test and compare the wave propagation characteristics in model runs from these four codes by measuring the dispersion relation of acoustic-gravity waves at various heights. We find considerable differences between the various models. The height dependence of wave power, in particular of high-frequency waves, varies by up to two orders of magnitude between the models, and the phase difference spectra of several models show unexpected features, including $\pm180^\circ$ phase jumps.

NASA's Interface Region Imaging Spectrograph (IRIS) provides high resolution observations of the solar atmosphere through UV spectroscopy and imaging. Since the launch of IRIS in June 2013, we have conducted systematic observation campaigns in coordination with the Swedish 1-m Solar Telescope (SST) on La Palma. The SST provides complementary high-resolution observations of the photosphere and chromosphere. The SST observations include spectro-polarimetric imaging in photospheric Fe I lines and spectrally-resolved imaging in the chromospheric Ca II 8542 A, H-alpha, and Ca II K lines. We present a database of co-aligned IRIS and SST datasets that is open for analysis to the scientific community. The database covers a variety of targets including active regions, sunspots, plage, quiet Sun, and coronal holes.

We analyze a grid of radiative hydrodynamic simulations of solar flares to study the energy balance and response of the atmosphere to nonthermal electron beam heating. The appearance of chromospheric bubbles is one of the most notable features that we find in the simulations. These pockets of chromospheric plasma get trapped between the transition region and the lower atmosphere as it is superheated by the particle beam. The chromospheric bubbles are seen in the synthetic spectra, appearing as an additional component to Balmer line profiles with high Doppler velocities as high as 200 km/s. Their signatures are also visible in the wings of Ca II 8542A line profiles. These bubbles of chromospheric plasma are driven upward by a wave front that is induced by the shock of energy deposition, and require a specific heating rate and atmospheric location to manifest.

Ambipolar diffusion is a physical mechanism related to the drift between charged and neutral particles in a partially ionized plasma that is key in many different astrophysical systems. However, understanding its effects is challenging due to basic uncertainties concerning relevant microphysical aspects and the strong constraints it imposes on the numerical modeling. Our aim is to introduce a numerical tool that allows us to address complex problems involving ambipolar diffusion in which, additionally, departures from ionization equilibrium are important or high resolution is needed. The primary application of this tool is for solar atmosphere calculations, but the methods and results presented here may also have a potential impact on other astrophysical systems. We have developed a new module for the stellar atmosphere Bifrost code that improves its computational capabilities of the ambipolar diffusion term in the Generalized Ohm's Law. This module includes, among other things, collision terms adequate to processes in the coolest regions in the solar chromosphere. As a key feature of the module, we have implemented the Super Time-Stepping (STS) technique, that allows an important acceleration of the calculations. We have also introduced hyperdiffusion terms to guarantee the stability of the code. We show that to have an accurate value for the ambipolar diffusion coefficient in the solar atmosphere it is necessary to include as atomic elements in the equation of state not only hydrogen and helium but also the main electron donors like sodium, silicon and potassium. In addition, we establish a range of criteria to set up an automatic selection of the free parameters of the STS method that guarantees the best performance, optimizing the stability and speed for the ambipolar diffusion calculations. We validate the STS implementation by comparison with a self-similar analytical solution.

Context. The pursuit of more realistic spectroscopic modelling and consistent abundances has led us to begin a new series of papers designed to improve current solar and stellar abundances of various atomic species. To achieve this, we have began updating the three-dimensional (3D) non-local thermodynamic equilibrium (non-LTE) radiative transfer code, Multi3D, and the equivalent one-dimensional (1D) non-LTE radiative transfer code, MULTI. Aims. We examine our improvements to these codes by redetermining the solar barium abundance. Barium was chosen for this test as it is an important diagnostic element of the s-process in the context of galactic chemical evolution. New Ba II + H collisional data for excitation and charge exchange reactions computed from first principles had recently become available and were included in the model atom. The atom also includes the effects of isotopic line shifts and hyperfine splitting. Method. A grid of 1D LTE barium lines were constructed with MULTI and fit to the four Ba II lines available to us in the optical region of the solar spectrum. Abundance corrections were then determined in 1D non-LTE, 3D LTE, and 3D non-LTE. A new 3D non-LTE solar barium abundance was computed from these corrections. Results. We present for the first time the full 3D non-LTE barium abundance of $A({\rm Ba})=2.27\pm0.02\pm0.01$, which was derived from four individual fully consistent barium lines. Errors here represent the systematic and random errors, respectively.

We present an initial study of one of the first ALMA Band 3 observations of the Sun with the aim to characterise the diagnostic potential of brightness temperatures measured with ALMA on the Sun. The observation covers 48min at a cadence of 2s targeting a Quiet Sun region at disk-centre. Corresponding time series of brightness temperature maps are constructed with the first version of the Solar ALMA Pipeline (SoAP) and compared to simultaneous SDO observations. The angular resolution of the observations is set by the synthesized beam (1.4x2.1as). The ALMA maps exhibit network patches, internetwork regions and also elongated thin features that are connected to large-scale magnetic loops as confirmed by a comparison with SDO maps. The ALMA Band 3 maps correlate best with the SDO/AIA 171, 131 and 304 channels in that they exhibit network features and, although very weak in the ALMA maps, imprints of large-scale loops. A group of compact magnetic loops is very clearly visible in ALMA Band 3. The brightness temperatures in the loop tops reach values of about 8000-9000K and in extreme moments up to 10 000K. ALMA Band 3 interferometric observations from early observing cycles already reveal temperature differences in the solar chromosphere. The weak imprint of magnetic loops and the correlation with the 171, 131, and 304 SDO channels suggests though that the radiation mapped in ALMA Band 3 might have contributions from a larger range of atmospheric heights than previously assumed but the exact formation height of Band 3 needs to be investigated in more detail. The absolute brightness temperature scale as set by Total Power measurements remains less certain and must be improved in the future. Despite these complications and the limited angular resolution, ALMA Band 3 observations have large potential for quantitative studies of the small-scale structure and dynamics of the solar chromosphere.

We examine the capabilities of a fast and simple method to infer line-of-sight (LOS) velocities from observations of the photospheric Si I 10827 A line. This spectral line is routinely observed together with the chromospheric He I 10830 A triplet as it helps to constrain the atmospheric parameters. We study the accuracy of bisector analysis and a line core fit of Si I 10827 A. We employ synthetic profiles starting from the Bifrost enhanced network simulation. The profiles are computed solving the radiative transfer equation, including non-local thermodynamic equilibrium effects on the determination of the atomic level populations of Si I. We found a good correlation between the inferred velocities from bisectors taken at different line profile intensities and the original simulation velocity at given optical depths. This good correlation means that we can associate bisectors taken at different line-profile percentages with atmospheric layers that linearly increase as we scan lower spectral line intensities. We also determined that a fit to the line-core intensity is robust and reliable, providing information about atmospheric layers that are above those accessible through bisectors. Therefore, by combining both methods on the Si I 10827 A line, we can seamlessly trace the quiet-Sun LOS velocity stratification from the deep photosphere to higher layers until around $\log \tau = -3.5$ in a fast and straightforward way. This method is ideal for generating quick-look reference images for future missions like the Daniel K. Inoue Solar Telescope and the European Solar Telescope, for example.

The European Solar Telescope (EST) is a research infrastructure for solar physics. It is planned to be an on-axis solar telescope with an aperture of 4 m and equipped with an innovative suite of spectro-polarimetric and imaging post-focus instrumentation. The EST project was initiated and is driven by EAST, the European Association for Solar Telescopes. EAST was founded in 2006 as an association of 14 European countries. Today, as of December 2019, EAST consists of 26 European research institutes from 18 European countries. The Preliminary Design Phase of EST was accomplished between 2008 and 2011. During this phase, in 2010, the first version of the EST Science Requirement Document (SRD) was published. After EST became a project on the ESFRI roadmap 2016, the preparatory phase started. The goal of the preparatory phase is to accomplish a final design for the telescope and the legal governance structure of EST. A major milestone on this path is to revisit and update the Science Requirement Document (SRD). The EST Science Advisory Group (SAG) has been constituted by EAST and the Board of the PRE-EST EU project in November 2017 and has been charged with the task of providing with a final statement on the science requirements for EST. Based on the conceptual design, the SRD update takes into account recent technical and scientific developments, to ensure that EST provides significant advancement beyond the current state-of-the-art. The present update of the EST SRD has been developed and discussed during a series of EST SAG meetings. The SRD develops the top-level science objectives of EST into individual science cases. Identifying critical science requirements is one of its main goals. Those requirements will define the capabilities of EST and the post-focus instrument suite. The technical requirements for the final design of EST will be derived from the SRD.

The thermal structure of the chromosphere is regulated through a complex interaction of various heating processes, radiative cooling, and the ionization degree of the plasma. Here we study the impact on the thermal properties of the chromosphere when including the combined action of non-equilibrium ionization (NEI) of hydrogen and helium and ion-neutral interaction effects. We have performed a 2.5D radiative magnetohydrodynamic simulation including ion-neutral interaction effects by solving the generalized Ohm's law (GOL) as well as NEI for hydrogen and helium using the Bifrost code. The GOL equation includes ambipolar diffusion and the Hall term. We compare this simulation with another simulation that computes the ionization in local thermodynamic equilibrium (LTE) including ion-neutral interaction effects. Our numerical models reveal substantial thermal differences in magneto-acoustic shocks, the wake behind the shocks, spicules, low-lying magnetic loops, and the transition region. In particular, we find that heating through ambipolar diffusion in shock wakes is substantially less efficient, while in the shock fronts themselves it is more efficient, under NEI conditions than when assuming LTE.

Magnetic flux emergence has been shown to be a key mechanism for unleashing a wide variety of solar phenomena. However, there are still open questions concerning the rise of the magnetized plasma through the atmosphere, mainly in the chromosphere, where the plasma departs from local thermodynamic equilibrium (LTE) and is partially ionized. We aim to investigate the impact of the nonequilibrium (NEQ) ionization and recombination and molecule formation of hydrogen, as well as ambipolar diffusion, on the dynamics and thermodynamics of the flux emergence process. Using the Bifrost code, we performed 2.5D numerical experiments of magnetic flux emergence from the convection zone up to the corona. The experiments include the NEQ ionization and recombination of atomic hydrogen, the NEQ formation and dissociation of H2 molecules, and the ambipolar diffusion term of the Generalized Ohm's Law. Our experiments show that the LTE assumption substantially underestimates the ionization fraction in most of the emerged region, leading to an artificial increase in the ambipolar diffusion and, therefore, in the heating and temperatures as compared to those found when taking the NEQ effects on the hydrogen ion population into account. We see that LTE also overestimates the number density of H2 molecules within the emerged region, thus mistakenly magnifying the exothermic contribution of the H2 molecule formation to the thermal energy during the flux emergence process. We find that the ambipolar diffusion does not significantly affect the amount of total unsigned emerged magnetic flux, but it is important in the shocks that cross the emerged region, heating the plasma on characteristic times ranging from 0.1 to 100 s. We also briefly discuss the importance of including elements heavier than hydrogen in the equation of state so as not to overestimate the role of ambipolar diffusion in the atmosphere.

Context. Velocity or intensity swirls have now been shown to be widely present throughout the photosphere and chromosphere. It has been suggested that these events could contribute to the heating of the upper solar atmosphere, via exciting Alfv\'en pulses, which could carry significant amounts of energy. However, the conjectured necessary physical conditions for their excitation, that the magnetic field rotates co-spatially and co-temporally with the velocity field, has not been verified. Aims. We aim to understand whether photospheric velocity swirls exist co-spatially and co-temporally with photospheric magnetic swirls, in order to demonstrate the link between swirls and pulses. Methods. The automated swirl detection algorithm (ASDA) is applied to the photospheric horizontal velocity and vertical magnetic fields obtained from a series of realistic numerical simulations using the RMHD code Bifrost. The spatial relationship between the detected velocity and magnetic swirls is further investigated via a well-defined correlation index (CI) study. Results. On average, there are ~63 short-lived photospheric velocity swirls (with lifetimes mostly less than 20 s, and average radius of ~37 km and rotating speeds of ~2.5 km/s) detected in a field of view (FOV) of 6\times6 Mm^{-2}, implying a total population of velocity swirls of ~1.06\times10^7 in the solar photosphere. More than 80% of the detected velocity swirls are found to be accompanied by local magnetic concentrations in intergranular lanes. On average, ~71% of the detected velocity swirls have been found to co-exist with photospheric magnetic swirls with the same rotating direction. Conclusions. The co-temporal and co-spatial rotation in the photospheric velocity and magnetic fields provide evidence that the conjectured condition for the excitation of Alfv\'en pulses by photospheric swirls is fulfilled.

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.

High-energy particles enter the solar atmosphere from Galactic or solar coronal sources, producing an "albedo'' source from the quiet Sun, now observable across a wide range of photon energies. The interaction of high-energy particles in a stellar atmosphere depends essentially upon the joint variation of the magnetic field and the gas, which heretofore has been characterized parametrically as P ~ B^alpha with P the gas pressure and B the magnitude of the magnetic field. We re-examine that parametrization by using a self-consistent 3D MHD model (Bifrost) and show that this relationship tends to P ~ B^{2.9+-0.1} based on the visible portions of the sample of open-field flux tubes in such a model, but with large variations from point to point. This scatter corresponds to the strong meandering of the open-field flux tubes in the lower atmosphere, which will have a strong effect on the prediction of the emission anisotropy (limb brightening). The simulations show that much of the open flux in coronal holes originates in weak-field regions within the granular pattern of the convective motions seen in the simulations.

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

The Spectral Imaging of the Coronal Environment (SPICE) instrument is a high-resolution imaging spectrometer operating at extreme ultraviolet (EUV) wavelengths. In this paper, we present the concept, design, and pre-launch performance of this facility instrument on the ESA/NASA Solar Orbiter mission. The goal of this paper is to give prospective users a better understanding of the possible types of observations, the data acquisition, and the sources that contribute to the instrument's signal. The paper discusses the science objectives, with a focus on the SPICE-specific aspects, before presenting the instrument's design, including optical, mechanical, thermal, and electronics aspects. This is followed by a characterisation and calibration of the instrument's performance. The paper concludes with descriptions of the operations concept and data processing. The performance measurements of the various instrument parameters meet the requirements derived from the mission's science objectives. The SPICE instrument is ready to perform measurements that will provide vital contributions to the scientific success of the Solar Orbiter mission.

Manganese (Mn) is a key Fe-group elements, commonly employed in stellar population and nucleosynthesis studies to explore the role of SN Ia. We have developed a new non-local thermodynamic equilibrium (NLTE) model of Mn, including new photo-ionisation cross-sections and new transition rates caused by collisions with H and H- atoms. We applied the model in combination with 1-dimensional (1D) LTE model atmospheres and 3D hydrodynamical simulations of stellar convection to quantify the impact of NLTE and convection on the line formation. We show that the effects of NLTE are present in Mn I and, to a lesser degree, in Mn II lines, and these increase with metallicity and with effective temperature of a model. Employing 3D NLTE radiative transfer, we derive new abundance of Mn in the Sun, A(Mn)=5.52 +/- 0.03 dex, consistent with the element abundance in C I meteorites. We also apply our methods to the analysis of three metal-poor benchmark stars. We find that 3D NLTE abundances are significantly higher than 1D LTE. For dwarfs, the differences between 1D NLTE and 3D NLTE abundances are typically within 0.15 dex, however, the effects are much larger in the atmospheres of giants owing to their more vigorous convection. We show that 3D NLTE successfully solves the ionisation and excitation balance for the RGB star HD 122563 that cannot be achieved by 1D LTE or 1D NLTE modelling. For HD 84937 and HD 140283, the ionisation balance is satisfied, however, the resonance Mn I triplet lines still show somewhat lower abundances compared to the high-excitation lines. Our results for the benchmark stars confirm that 1D LTE modelling leads to significant systematic biases in Mn abundances across the full wavelength range from the blue to the IR. We also produce a list of Mn lines that are not significantly biased by 3D and can be reliably, within the 0.1 dex uncertainty, modelled in 1D NLTE.

The solar Ly$\alpha$ line is the strongest line in the ultraviolet waveband, and is greatly enhanced during solar flares. Here we present radiative hydrodynamic simulations of solar flares under different heating models, and calculate the response of this line taking into account non-equilibrium ionization of hydrogen and partial frequency redistribution. We find that in non-thermal heating models, the Ly$\alpha$ line can show a red or blue asymmetry corresponding to the chromospheric evaporation or condensation, respectively. The asymmetry may change from red to blue if the electron beam flux is large enough to produce a significant chromospheric condensation region. In the Ly$\alpha$ intensity lightcurve, there appears a dip when the change of asymmetry occurs. In thermal models, the Ly$\alpha$ line intensity peaks quickly and then falls, and the profile has an overall red asymmetry, which is similar to the profiles from heating by a soft electron beam. The Ly$\alpha$ profile shows a single red peak at the end of thermal heating, and the whole line is formed in a very small height range.

Ellerman bombs and UV bursts are transient brightenings that are ubiquitously observed in the lower atmospheres of active and emerging flux regions. Here we present inversion results of SST/CRISP and CHROMIS, as well as IRIS data of such transient events. Combining information from the Mg II h & k, Si IV and Ca II 8542A and Ca II H & K lines, we aim to characterise their temperature and velocity stratification, as well as their magnetic field configuration. We find average temperature enhancements of a few thousand kelvin close to the classical temperature minimum, but localised peak temperatures of up to 10,000-15,000 K from Ca II inversions. Including Mg II generally dampens these temperature enhancements to below 8000 K, while Si IV requires temperatures in excess of 10,000 K at low heights, but may also be reproduced with secondary temperature enhancements of 35,000-60,000 K higher up. However, reproducing Si IV comes at the expense of overestimating the Mg II emission. The line-of-sight velocity maps show clear bi-directional jet signatures and strong correlation with substructure in the intensity images, with slightly larger velocities towards the observer than away. The magnetic field parameters show an enhancement of the horizontal field co-located with the brightenings at similar heights as the temperature increase. We are thus able to largely reproduce the observational properties of Ellerman bombs with UV burst signature with temperature stratifications peaking close to the classical temperature minimum. Correctly modelling the Si IV emission in agreement with all other diagnostics is, however, an outstanding issue. Accounting for resolution differences, fitting localised temperature enhancements and/or performing spatially-coupled inversions is likely necessary to obtain better agreement between all considered diagnostics.

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.

Ellerman bombs (EBs) and UV bursts are both brightenings related to flux emergence regions and specifically to magnetic flux of opposite polarity that meet in the photosphere. These two reconnection-related phenomena, nominally formed far apart, occasionally occur in the same location and at the same time, thus challenging our understanding of reconnection and heating of the lower solar atmosphere. We consider the formation of an active region, including long fibrils and hot and dense coronal plasma. The emergence of a untwisted magnetic flux sheet, injected $2.5$~Mm below the photosphere, is studied as it pierces the photosphere and interacts with the preexisting ambient field. Specifically, we aim to study whether EBs and UV bursts are generated as a result of such flux emergence and examine their physical relationship. The Bifrost radiative magnetohydrodynamics code was used to model flux emerging into a model atmosphere that contained a fairly strong ambient field, constraining the emerging field to a limited volume wherein multiple reconnection events occur as the field breaks through the photosphere and expands into the outer atmosphere. Synthetic spectra of the different reconnection events were computed using the $1.5$D RH code and the fully 3D MULTI3D code. The formation of UV bursts and EBs at intensities and with line profiles that are highly reminiscent of observed spectra are understood to be a result of the reconnection of emerging flux with itself in a long-lasting current sheet that extends over several scale heights through the chromosphere. Synthetic diagnostics suggest that there are no compelling reasons to assume that UV bursts occur in the photosphere. Instead, EBs and UV bursts are occasionally formed at opposite ends of a long current sheet that resides in an extended bubble of cool gas.

We investigate the diagnostic potential of the spectral lines at 850 nm for understanding the magnetism of the lower atmosphere. For that purpose, we use a newly developed 3D simulation of a chromospheric jet to check the sensitivity of the spectral lines to this phenomenon as well as our ability to infer the atmospheric information through spectropolarimetric inversions of noisy synthetic data. We start comparing the benefits of inverting the entire spectrum at 850 nm versus only the Ca II 8542 A spectral line. We found a better match of the input atmosphere for the former case, mainly at lower heights. However, the results at higher layers were not accurate. After several tests, we determined that we need to weight more the chromospheric lines than the photospheric ones in the computation of the goodness of the fit. The new inversion configuration allows us to obtain better fits and consequently more accurate physical parameters. Therefore, to extract the most from multi-line inversions, a proper set of weights needs to be estimated. Besides that, we conclude again that the lines at 850 nm, or a similar arrangement with Ca II 8542 A plus Zeeman sensitive photospheric lines, poses the best observing configuration for examining the thermal and magnetic properties of the lower solar atmosphere.

The presence of magnetic field is crucial in the transport of energy through the solar atmosphere. Recent ground-based and space-borne observations of the quiet Sun have revealed that magnetic field accumulates at photospheric heights, via a local dynamo or from small-scale flux emergence events. However, most of this small-scale magnetic field may not expand into the chromosphere due to the entropy drop with height at the photosphere. Here we present a study that uses a high resolution 3D radiative MHD simulation of the solar atmosphere with non-grey and non-LTE radiative transfer and thermal conduction along the magnetic field to reveal that: 1) the net magnetic flux from the simulated quiet photosphere is not sufficient to maintain a chromospheric magnetic field (on average), 2) processes in the lower chromosphere, in the region dominated by magneto-acoustic shocks, are able to convert kinetic energy into magnetic energy, 3) the magnetic energy in the chromosphere increases linearly in time until the r.m.s. of the magnetic field strength saturates at roughly 4 to 30 G (horizontal average) due to conversion from kinetic energy, 4) and that the magnetic features formed in the chromosphere are localized to this region.

Solar observations with the Atacama Large Millimeter/submillimeter Array (ALMA) provide us with direct measurements of the brightness temperature in the solar chromosphere. We study the temperature distributions obtained with ALMA Band 6 (in four sub-bands at 1.21, 1.22, 1.29, and 1.3 mm) for various areas at, and in the vicinity of, a sunspot, comprising quasi-quiet and active regions with different amounts of underlying magnetic fields. We compare these temperatures with those obtained at near- and far-ultraviolet (UV) wavelengths (and with the line-core intensities of the optically-thin far-UV spectra), co-observed with the Interface Region Imaging Spectrograph (IRIS) explorer. These include the emission peaks and cores of the Mg II k 279.6 nm and Mg II h 280.4 nm lines as well as the line cores of C II 133.4 nm, O I 135.6 nm, and Si IV 139.4 nm, sampling the mid-to-high chromosphere and the low transition region. Splitting the ALMA sub-bands resulted in an slight increase of spatial resolution in individual temperature maps, thus, resolving smaller-scale structures compared to those produced with the standard averaging routines. We find that the radiation temperatures have different, though somewhat overlapping, distributions in different wavelengths and in the various magnetic regions. Comparison of the ALMA temperatures with those of the UV diagnostics should, however, be interpreted with great caution, the former is formed under the local thermodynamic equilibrium (LTE) conditions, the latter under non-LTE. The mean radiation temperature of the ALMA Band 6 is similar to that extracted from the IRIS C II line in all areas with exception of the sunspot and pores where the C II poses higher radiation temperatures. In all magnetic regions, the Mg II lines associate with the lowest mean radiation temperatures in our sample. These will provide constraints for future numerical models.

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.

We study the polarization signals and physical parameters of penumbral microjets (PMJs) by using high spatial resolution data taken in the Fe I 630 nm pair, Ca II 854.2 nm and Ca II K lines with the CRISP and CHROMIS instruments at the Swedish 1-m Solar Telescope. We infer their physical parameters, such as physical observables in the photosphere and chromospheric velocity diagnostics, by different methods, including inversions of the observed Stokes profiles with the STiC code. PMJs harbor overall brighter Ca II K line profiles and conspicuous polarization signals in Ca II 854.2 nm, specifically in circular polarization that often shows multiple lobes mainly due to the shape of Stokes I. They usually overlap photospheric regions with sheared magnetic field configuration, suggesting that magnetic reconnections could play an important role in the origin of PMJs. The discrepancy between their low LOS velocities and the high apparent speeds reported on earlier, as well as the existence of different vertical velocity gradients in the chromosphere, indicate that PMJs might not be entirely related to mass motions. Instead, PMJs could be due to perturbation fronts induced by magnetic reconnections occurring in the deep photosphere that propagate through the chromosphere. This reconnection may be associated with current heating that produces temperature enhancements from the temperature minimum region. Furthermore, enhanced collisions with electrons could also increase the coupling to the local conditions at higher layers during the PMJ phase, giving a possible explanation for the enhanced emission in the overall Ca II K profiles emerging from these transients.

The next generation of solar observatories aim to understand the magnetism of the solar chromosphere. Therefore, it is crucial to understand the polarimetric signatures of chromospheric spectral lines. For this purpose, we here examine the suitability of the three Fraunhofer Mg I b1, b2, and b4 lines at 5183.6, 5172.7, and 5167.3 A, respectively. We start by describing a simplified atomic model of only 6 levels and 3 line transitions for computing the atomic populations of the 3p-4s (multiplet number 2) levels involved in the Mg I b line transitions assuming non-local thermodynamic conditions and considering only the Zeeman effect using the field-free approximation. We test this simplified atom against more complex ones finding that, although there are differences in the computed profiles, they are small compared with the advantages provided by the simple atom in terms of speed and robustness. After comparing the three Mg I lines, we conclude that the most capable one is the b2 line as b1 forms at similar heights and always show weaker polarization signals while b4 is severely blended with photospheric lines. We also compare Mg I b2 with the K I D1 and Ca II 8542 A lines finding that the former is sensitive to the atmospheric parameters at heights that are in between those covered by the latter two lines. This makes Mg I b2 an excellent candidate for future multi-line observations that aim to seamlessly infer the thermal and magnetic properties of different features in the lower solar atmosphere.

The Chromospheric Lyman-Alpha Spectro-Polarimeter (CLASP) is a suborbital rocket experiment that on 3rd September 2015 measured the linear polarization produced by scattering processes in the hydrogen Ly-$\alpha$ line of the solar disk radiation, whose line-center photons stem from the chromosphere-corona transition region (TR). These unprecedented spectropolarimetric observations revealed an interesting surprise, namely that there is practically no center-to-limb variation (CLV) in the $Q/I$ line-center signals. Using an analytical model, we first show that the geometrical complexity of the corrugated surface that delineates the TR has a crucial impact on the CLV of the $Q/I$ and $U/I$ line-center signals. Secondly, we introduce a statistical description of the solar atmosphere based on a three-dimensional (3D) model derived from a state-of-the-art radiation magneto-hydrodynamic simulation. Each realization of the statistical ensemble is a 3D model characterized by a given degree of magnetization and corrugation of the TR, and for each such realization we solve the full 3D radiative transfer problem taking into account the impact of the CLASP instrument degradation on the calculated polarization signals. Finally, we apply the statistical inference method presented in a previous paper to show that the TR of the 3D model that produces the best agreement with the CLASP observations has a relatively weak magnetic field and a relatively high degree of corrugation. We emphasize that a suitable way to validate or refute numerical models of the upper solar chromosphere is by confronting calculations and observations of the scattering polarization in ultraviolet lines sensitive to the Hanle effect.

We study the M1.9 class solar flare SOL2015-09-27T10:40 UT using high-resolution full-Stokes imaging spectropolarimetry of the Ca ii 8542 {\AA} line obtained with the CRISP imaging spectropolarimeter at the Swedish 1-m Solar Telescope. Spectropolarimetric inversions using the non-LTE code NICOLE are used to construct semi-empirical models of the flaring atmosphere to investigate the structure and evolution of the flare temperature and magnetic field. A comparison of the temperature stratification in flaring and non-flaring areas reveals strong heating of the flare ribbon during the flare peak. The polarization signals of the ribbon in the chromosphere during the flare maximum become stronger when compared to its surroundings and to pre- and post- flare profiles. Furthermore, a comparison of the response functions to perturbations in the line-of-sight magnetic field and temperature in flaring and non-flaring atmospheres shows that during the flare the Ca ii 8542 {\AA} line is more sensitive to the lower atmosphere where the magnetic field is expected to be stronger. The chromospheric magnetic field was also determined with the weak-field approximation which led to results similar to those obtained with the NICOLE inversions.

On 3rd September 2015, the Chromospheric Lyman-Alpha SpectroPolarimeter (CLASP) successfully measured the linear polarization produced by scattering processes in the hydrogen Lyman-$\alpha$ line of the solar disk radiation, revealing conspicuous spatial variations in the $Q/I$ and $U/I$ signals. Via the Hanle effect the line-center $Q/I$ and $U/I$ amplitudes encode information on the magnetic field of the chromosphere-corona transition region (TR), but they are also sensitive to the three-dimensional structure of this corrugated interface region. With the help of a simple line formation model, here we propose a statistical inference method for interpreting the Lyman-$\alpha$ line-center polarization observed by CLASP.

The role of type II spicules in the corona has been a much debated topic in recent years. This paper aims to shed light on the impact of type II spicules in the corona using novel 2.5D radiative MHD simulations including ion-neutral interaction effects with the Bifrost code. We find that the formation of simulated type II spicules, driven by the release of magnetic tension, impacts the corona in various manners. Associated with the formation of spicules, the corona exhibits 1) magneto-acoustic shocks and flows which supply mass to coronal loops, and 2) transversal magnetic waves and electric currents that propagate at Alfv\'en speeds. The transversal waves and electric currents, generated by the spicule's driver and lasting for many minutes, are dissipated and heat the associated loop. These complex interactions in the corona can be connected with blue shifted secondary components in coronal spectral lines (Red-Blue asymmetries) observed with Hinode/EIS and SOHO/SUMER, as well as the EUV counterpart of type II spicules and propagating coronal disturbances (PCDs) observed with the 171~\AA\ and 193~\AA\ SDO/AIA channels.

The measured average velocities in solar and stellar spectral lines formed at transition region temperatures have been difficult to interpret. However, realistic three-dimensional radiation magnetohydrodynamics (3D rMHD) models of the solar atmosphere are able to reproduce the observed dominant line shifts and may thus hold the key to resolve these issues. Our new 3D rMHD simulations aim to shed light on how mass flows between the chromosphere and corona and on how the coronal mass is maintained. Passive tracer particles, so-called corks, allow the tracking of parcels of plasma over time and thus the study of changes in plasma temperature and velocity not only locally, but also in a co-moving frame. By following the trajectories of the corks, we can investigate mass and energy flows and understand the composition of the observed velocities. Our findings show that most of the transition region mass is cooling. The preponderance of transition region redshifts in the model can be explained by the higher percentage of downflowing mass in the lower and middle transition region. The average upflows in the upper transition region can be explained by a combination of both stronger upflows than downflows and a higher percentage of upflowing mass. The most common combination at lower and middle transition region temperatures are corks that are cooling and traveling downward. For these corks, a strong correlation between the pressure gradient along the magnetic field line and the velocity along the magnetic field line has been observed, indicating a formation mechanism that is related to downward propagating pressure disturbances. Corks at upper transition region temperatures are subject to a rather slow and highly variable but continuous heating process.

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

The Fe I 6173 {\AA} line is widely used in the measurements of vector magnetic fields by instruments including the Helioseismic and Magnetic Imager (HMI). We perform non-local thermodynamic equilibrium calculations of this line based on radiative hydrodynamic simulations in a flaring atmosphere. We employ both a quiet-Sun atmosphere and a penumbral atmosphere as the initial one in our simulations. We find that, in the quiet-Sun atmosphere, the line center is obviously enhanced during an intermediate flare. The enhanced emission is contributed from both radiative backwarming in the photosphere and particle beam heating in the lower chromosphere. A blue asymmetry of the line profile also appears due to an upward mass motion in the lower chromosphere. If we take a penumbral atmosphere as the initial atmosphere, the line has a more significant response to the flare heating, showing a central emission and an obvious asymmetry. The low spectral resolution of HMI would indicate some loss of information but the enhancement and line asymmetry are still kept. By calculating polarized line profiles, we find that the Stokes I and V profiles can be altered as a result of flare heating. Thus the distortion of this line has a crucial influence on the magnetic field measured from this line, and one should be cautious in interpreting the magnetic transients observed frequently in solar flares.

The heating of the solar chromosphere remains one of the most important questions in solar physics. Our current understanding is that small-scale internetwork (IN) magnetic fields play an important role as a heating agent. Indeed, cancellations of IN magnetic elements in the photosphere can produce transient brightenings in the chromosphere and transition region. These bright structures might be the signature of energy release and plasma heating, probably driven by magnetic reconnection of IN field lines. Although single events are not expected to release large amounts of energy, their global contribution to the chromosphere may be significant due to their ubiquitous presence in quiet Sun regions. In this paper we study cancellations of IN elements and analyze their impact on the energetics and dynamics of the quiet Sun atmosphere. We use high resolution, multiwavelength, coordinated observations obtained with the Interface Region Imaging Spectrograph (IRIS) and the Swedish 1-m Solar Telescope (SST) to identify cancellations of IN magnetic flux patches and follow their evolution. We find that, on average, these events live for ~3 minutes in the photosphere and ~12 minutes in the chromosphere and/or transition region. Employing multi-line inversions of the Mg II h & k lines we show that cancellations produce clear signatures of heating in the upper atmospheric layers. However, at the resolution and sensitivity accessible to the SST, their number density still seems to be one order of magnitude too low to explain the global chromospheric heating.

CHROMIS, a new imaging spectrometer at the Swedish 1-m Solar Telescope (SST), can observe the chromosphere in the H and K lines of Ca II at high spatial and spectral resolution. Accurate modeling as well as an understanding of the formation of these lines are needed to interpret the SST/CHROMIS observations. Such modeling is computationally challenging because these lines are influenced by strong departures from local thermodynamic equilibrium, three-dimensional radiative transfer, and partially coherent resonance scattering of photons. We aim to model the CaII H&K lines in 3D model atmospheres to understand their formation and to investigate their diagnostic potential for probing the chromosphere. We model the synthetic spectrum of Ca II using the radiative transfer code Multi3D in three different radiation-magnetohydrodynamic model atmospheres computed with the Bifrost code. We classify synthetic intensity profiles according to their shapes and study how their features are related to the physical properties in the model atmospheres. We investigate whether the synthetic data reproduce the observed spatially-averaged line shapes, center-to-limb variation and compare with SST/CHROMIS images. The spatially-averaged synthetic line profiles show too low central emission peaks, and too small separation between the peaks. The trends of the observed center-to-limb variation of the profiles properties are reproduced by the models. The Ca II H&K line profiles provide a temperature diagnostic of the temperature minimum and the temperature at the formation height of the emission peaks. The Doppler shift of the central depression is an excellent probe of the velocity in the upper chromosphere.

Context. The radiative losses in the solar chromosphere vary from 4~kW~m$^{-2}$ in the quiet Sun, to 20~kW~m$^{-2}$ in active regions. The mechanisms that transport non-thermal energy to and deposit it in the chromosphere are still not understood. Aims. We aim to investigate the atmospheric structure and heating of the solar chromosphere in an emerging flux region. Methods. We use observations taken with the CHROMIS and CRISP instruments on the Swedish 1-m Solar Telescope in the Ca II K, Ca II 854.2 nm, H$\alpha$, and Fe I 630.1 nm and 630.2 nm lines. We analyse the various line profiles and in addition perform multi-line, multi-species, non-Local Thermodynamic Equilibrium (non-LTE) inversions to estimate the spatial and temporal variation of the chromospheric structure. Results. We investigate which spectral features of Ca II K contribute to the frequency-integrated Ca II K brightness, which we use as a tracer of chromospheric radiative losses. The majority of the radiative losses are not associated with localized high-Ca II K-brightness events, but instead with a more gentle, spatially extended, and persistent heating. The frequency-integrated Ca II K brightness correlates strongly with the total linear polarization in the Ca II 854.2 nm line, while the Ca II K profile shapes indicate that the bulk of the radiative losses occur in the lower chromosphere. Non-LTE inversions indicate a transition from heating concentrated around photospheric magnetic elements below $\log{\tau_{500}} =-3$ to a more space-filling and time-persistent heating above $\log{\tau_{500}} =-4$. The inferred gas temperature at $\log{\tau_{500}} =-3.8$ correlates strongly with the total linear polarization in the Ca II 854.2 nm line, suggesting that that the heating rate correlates with the strength of the horizontal magnetic field in the low chromosphere.

Magnetic reconnection is thought to drive a wide variety of dynamic phenomena in the solar atmosphere. Yet the detailed physical mechanisms driving reconnection are difficult to discern in the remote sensing observations that are used to study the solar atmosphere. In this paper we exploit the high-resolution instruments Interface Region Imaging Spectrograph (IRIS) and the new CHROMIS Fabry-Perot instrument at the Swedish 1-m Solar Telescope (SST) to identify the intermittency of magnetic reconnection and its association with the formation of plasmoids in so-called UV bursts in the low solar atmosphere. The Si IV 1403A UV burst spectra from the transition region show evidence of highly broadened line profiles with often non-Gaussian and triangular shapes, in addition to signatures of bidirectional flows. Such profiles had previously been linked, in idealized numerical simulations, to magnetic reconnection driven by the plasmoid instability. Simultaneous CHROMIS images in the chromospheric Ca II K 3934A line now provide compelling evidence for the presence of plasmoids, by revealing highly dynamic and rapidly moving brightenings that are smaller than 0.2 arcsec and that evolve on timescales of order seconds. Our interpretation of the observations is supported by detailed comparisons with synthetic observables from advanced numerical simulations of magnetic reconnection and associated plasmoids in the chromosphere. Our results highlight how subarcsecond imaging spectroscopy sensitive to a wide range of temperatures combined with advanced numerical simulations that are realistic enough to compare with observations can directly reveal the small-scale physical processes that drive the wide range of phenomena in the solar atmosphere.

In the lower solar atmosphere, the chromosphere is permeated by jets known as spicules, in which plasma is propelled at speeds of 50 to 150 kilometers per second into the corona. The origin of the spicules is poorly understood, although they are expected to play a role in heating the million-degree corona and are associated with Alfv\'enic waves that help drive the solar wind. We compare magnetohydrodynamic simulations of spicules with observations from the Interface Region Imaging Spectrograph and the Swedish 1-m Solar Telescope. Spicules are shown to occur when magnetic tension is amplified and transported upward through interactions between ions and neutrals or ambipolar diffusion. The tension is impulsively released to drive flows, heat plasma (through ambipolar diffusion), and generate Alfv\'enic waves.

Observations of the upper chromosphere shows an enormous amount of intricate fine structure. Much of this comes in the form of linear features which are most often assumed to be well aligned with the direction of the magnetic field in the low plasma beta regime thought to dominate the upper chromosphere. We use advanced radiative MHD simulations including the effects of ion-neutral interactions (using the generalized Ohm's law) in the partially ionized chromosphere to show that the magnetic field is often not well aligned with chromospheric features. This occurs where the ambipolar diffusion is large, i.e., ions and neutral populations decouple as the ion-neutral collision frequency drops allowing the field to slip through the neutral population, currents perpendicular to the field are strong, and thermodynamic timescales are longer than or similar to the those of ambipolar diffusion. We find this often happens in dynamic spicule or fibril-like features at the top of the chromosphere. This has important consequences for field extrapolation methods which increasingly use such upper chromospheric features to help constrain the chromospheric magnetic field: our results invalidate the underlying assumption that these features are aligned with the field. In addition, our results cast doubt on results from 1D hydrodynamic models, which assume that plasma remains on the same field lines. Finally, our simulations show that ambipolar diffusion significantly alters the amount of free energy available in the coronal part of our simulated volume, which is likely to have consequences for studies of flare initiation.

The Hydrogen Balmer series is a basic radiative loss channel from the flaring solar chromosphere. We report here on the analysis of an extremely rare set of simultaneous observations of a solar flare in the H${\alpha}$ and H${\beta}$ lines at high spatial and temporal resolution, which were acquired at the Dunn Solar Telescope. Images of the C3.3 flare (SOL2014-04-22T15:22) made at various wavelengths along the H${\alpha}$ line profile by the Interferometric Bidimensional Spectrometer (IBIS) and in the H${\beta}$ with the Rapid Oscillations in the Solar Atmosphere (ROSA) broadband imager are analyzed to obtain the intensity evolution. The H${\alpha}$ and H${\beta}$ intensity excesses in three identified flare footpoints are well correlated in time. We examine the ratio of H${\alpha}$ to H${\beta}$ flare excess, which was proposed by previous authors as a possible diagnostic of the level of electron beam energy input. In the stronger footpoints, the typical value of the the H${\alpha}$/H${\beta}$ intensity ratio observed is $\sim 0.4-0.5$, in broad agreement with values obtained from a RADYN non-LTE simulation driven by an electron beam with parameters constrained (as far as possible) by observation. The weaker footpoint has a larger H${\alpha}$/H${\beta}$ ratio, again consistent with a RADYN simulation but with a smaller energy flux. The H${\alpha}$ line profiles observed have a less prominent central reversal than is predicted by the RADYN results, but can be brought into agreement if the H${\alpha}$-emitting material has a filling factor of around 0.2--0.3.

The C I 135.58 line is located in the wavelength range of NASA's Interface Region Imagin Spectrograph (IRIS) small explorer mission. We here study the formation and diagnostic potential of this line by means of non local-thermodynamic-equilibrium modeling, employing both 1D and 3D radiation-magnetohydrodynamic models. The C I/C II ionization balance is strongly influenced by photoionization by Ly-alpha emission. The emission in the C I 135.58 line is dominated by a recombination cascade and the line forming region is optically thick. The Doppler shift of the line correlates strongly with the vertical velocity in its line forming region, which is typically located at 1.5 Mm height. With IRIS the C I 135.58 line is usually observed together with the O I 135.56 line, and from the Doppler shift of both lines, we obtain the velocity difference between the line forming regions of the two lines. From the ratio of the C I/O I line core intensity, we can determine the distance between the C I and the O I forming layers. Combined with the velocity difference, the velocity gradient at mid-chromospheric heights can be derived. The C I/O I total intensity line ratio is correlated with the inverse of the electron density in the mid-chromosphere. We conclude that the C I 135.58 line is an excellent probe of the middle chromosphere by itself, and together with the O I 135.56 line the two lines provide even more information, which complements other powerful chromospheric diagnostics of IRIS such as the Mg II h and k lines and the C II lines around 133.5 nm.

We investigate the effects of interactions between ions and neutrals on the chromosphere and overlying corona using 2.5D radiative MHD simulations with the Bifrost code. We have extended the code capabilities implementing ion-neutral interaction effects using the Generalized Ohm's Law, i.e., we include the Hall term and the ambipolar diffusion (Pedersen dissipation) in the induction equation. Our models span from the upper convection zone to the corona, with the photosphere, chromosphere and transition region partially ionized. Our simulations reveal that the interactions between ionized particles and neutral particles have important consequences for the magneto-thermodynamics of these modeled layers: 1) ambipolar diffusion increases the temperature in the chromosphere; 2) sporadically the horizontal magnetic field in the photosphere is diffused into the chromosphere due to the large ambipolar diffusion; 3) ambipolar diffusion concentrates electrical currents leading to more violent jets and reconnection processes, resulting in 3a) the formation of longer and faster spicules, 3b) heating of plasma during the spicule evolution, and 3c) decoupling of the plasma and magnetic field in spicules. Our results indicate that ambipolar diffusion is a critical ingredient for understanding the magneto-thermo-dynamic properties in the chromosphere and transition region. The numerical simulations have been made publicly available, similar to previous Bifrost simulations. This will allow the community to study realistic numerical simulations with a wider range of magnetic field configurations and physics modules than previously possible.

In this publication we continue the work started in Quintero Noda et al. (2017) examining this time a numerical simulation of a magnetic flux tube concentration. Our goal is to study if the physical phenomena that take place in it, in particular, the magnetic pumping, leaves a specific imprint on the examined spectral lines. We find that the profiles from the interior of the flux tube are periodically dopplershifted following an oscillation pattern that is also reflected in the amplitude of the circular polarization signals. In addition, we analyse the properties of the Stokes profiles at the edges of the flux tube discovering the presence of linear polarization signals for the Ca II lines, although they are weak with an amplitude around 0.5% of the continuum intensity. Finally, we compute the response functions to perturbations in the longitudinal field and we estimate the field strength using the weak field approximation. Our results indicate that the height of formation of the spectral lines changes during the magnetic pumping process which makes the interpretation of the inferred magnetic field strength and its evolution more difficult. These results complement those from previous works demonstrating the capabilities and limitations of the 850 nm spectrum for chromospheric Zeeman polarimetry in a very dynamic and complex atmosphere.

Ellerman bombs (EBs) are brightenings in the H$\alpha$ line wings that are believed to be caused by magnetic reconnection in the lower atmosphere. To study the response and evolution of the chromospheric line profiles, we perform radiative hydrodynamic simulations of EBs using both non-thermal and thermal models. Overall, these models can generate line profiles that are similar to observations. However, in non-thermal models we find dimming in the H$\alpha$ line wings and continuum when the heating begins, while for the thermal models dimming occurs only in the H$\alpha$ line core, and with a longer lifetime. This difference in line profiles can be used to determine whether an EB is dominated by non-thermal heating or thermal heating. In our simulations, if a higher heating rate is applied, the H$\alpha$ line will be unrealistically strong, while there are still no clear UV burst signatures.

We characterize the K I D1 & D2 lines in order to determine whether they could complement the 850 nm window, containing the Ca II infrared triplet lines and several Zeeman sensitive photospheric lines, that was studied previously. We investigate the effect of partial redistribution on the intensity profiles, their sensitivity to changes in different atmospheric parameters, and the spatial distribution of Zeeman polarization signals employing a realistic magnetohydrodynamic simulation. The results show that these lines form in the upper photosphere at around 500 km and that they are sensitive to the line of sight velocity and magnetic field strength at heights where neither the photospheric lines nor the Ca II infrared lines are. However, at the same time, we found that their sensitivity to the temperature essentially comes from the photosphere. Then, we conclude that the K I lines provide a complement to the lines in the 850 nm window for the determination of atmospheric parameters in the upper photosphere, especially for the line of sight velocity and the magnetic field.

There is a thin transition region (TR) in the solar atmosphere where the temperature rises from 10,000 K in the chromosphere to millions of degrees in the corona. Little is known about the mechanisms that dominate this enigmatic region other than the magnetic field plays a key role. The magnetism of the TR can only be detected by polarimetric measurements of a few ultraviolet (UV) spectral lines, the Lyman-$\alpha$ line of neutral hydrogen at 121.6 nm (the strongest line of the solar UV spectrum) being of particular interest given its sensitivity to the Hanle effect (the magnetic-field-induced modification of the scattering line polarization). We report the discovery of linear polarization produced by scattering processes in the Lyman-$\alpha$ line, obtained with the Chromospheric Lyman-Alpha Spectro-Polarimeter (CLASP) rocket experiment. The Stokes profiles observed by CLASP in quiet regions of the solar disk show that the $Q/I$ and $U/I$ linear polarization signals are of the order of 0.1 % in the line core and up to a few percent in the nearby wings, and that both have conspicuous spatial variations with scales of $\sim 10$ arcsec. These observations help constrain theoretical models of the chromosphere-corona TR and extrapolations of the magnetic field from photospheric magnetograms. In fact, the observed spatial variation from disk to limb of polarization at the line core and wings already challenge the predictions from three-dimensional magnetohydrodynamical models of the upper solar chromosphere.

A spectacular manifestation of solar activity is the appearance of transient brightenings in the far wings of the H$\alpha$ line, known as Ellerman bombs (EBs). Recent observations obtained by the Interface Region Imaging Spectrograph (IRIS) have revealed another type of plasma "bombs" (UV bursts) with high temperatures of perhaps up to $8\times 10^4$ K within the cooler lower solar atmosphere. Realistic numerical modeling showing such events is needed to explain their nature. Here, we report on 3D radiative magneto-hydrodynamic simulations of magnetic flux emergence in the solar atmosphere. We find that ubiquitous reconnection between emerging bipolar magnetic fields can trigger EBs in the photosphere, UV bursts in the mid/low chromosphere and small (nano-/micro-) flares ($10^6$ K) in the upper chromosphere. These results provide new insights on the emergence and build up of the coronal magnetic field and the dynamics and heating of the solar surface and lower atmosphere.

We use state-of-the-art, three-dimensional non-local thermodynamic equilibrium (non-LTE) radiative magnetohydrodynamic simulations of the quiet solar atmosphere to carry out detailed tests of chromospheric magnetic field diagnostics from free-free radiation at millimeter and submillimeter wavelengths (mm/submm). The vertical component of the magnetic field was deduced from the mm/submm brightness spectra and the degree of circular polarization synthesized at millimeter frequencies. We used the frequency bands observed by the Atacama Large Millimeter/Submillimeter Array (ALMA) as a convenient reference. The magnetic field maps obtained describe the longitudinal magnetic field at the effective formation heights of the relevant wavelengths in the solar chromosphere. The comparison of the deduced and model chromospheric magnetic fields at the spatial resolution of both the model and current observations demonstrates a good correlation, but has a tendency to underestimate the model field. The systematic discrepancy of about 10 percent is probably due to averaging of the restored field over the heights contributing to the radiation, weighted by the strength of the contribution. On the whole, the method of probing the longitudinal component of the magnetic field with free-free emission at mm/submm wavelengths is found to be applicable to measurements of the weak quiet-Sun magnetic fields. However, successful exploitation of this technique requires very accurate measurements of the polarization properties (primary beam and receiver polarization response) of the antennas, which will be the principal factor that determines the level to which chromospheric magnetic fields can be measured. Consequently, high-resolution and high-precision observations of circularly polarized radiation at millimeter wavelengths can be a powerful tool for producing chromospheric longitudinal magnetograms.

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.

Future solar missions and ground-based telescopes aim to understand the magnetism of the solar chromosphere. We performed a supporting study in Quintero Noda et al. (2016) focused on the infrared Ca II 8542 A line and we concluded that is one of the best candidates because it is sensitive to a large range of atmospheric heights, from the photosphere to the middle chromosphere. However, we believe that it is worth to try improving the results produced by this line observing additional spectral lines. In that regard, we examined the neighbour solar spectrum looking for spectral lines that could increase the sensitivity to the atmospheric parameters. Interestingly, we discovered several photospheric lines that greatly improve the photospheric sensitivity to the magnetic field vector. Moreover, they are located close to a second chromospheric line that also belongs to the Ca II infrared triplet, i.e. the Ca II 8498 A line, and enhances the sensitivity to the atmospheric parameters at chromospheric layers. We conclude that the lines in the vicinity of the Ca II 8542 A line not only increase its sensitivity to the atmospheric parameters at all layers, but also they constitute an excellent spectral window for chromospheric polarimetry.

Context: While classical models successfully reproduce intensities of many transition region lines, they predict helium EUV line intensities roughly an order of magnitude lower than the observed value. Aims: To determine the relevant formation mechanism(s) of the helium EUV resonance lines, capable of explaining the high intensities under quiet sun conditions. Methods: We synthesise and study the emergent spectra from a 3D radiation-magnetohydrodynamics simulation model. The effects of coronal illumination and non-equilibrium ionisation of hydrogen and helium are included self-consistently in the numerical simulation. Results: Radiative transfer calculations result in helium EUV line intensities that are an order of magnitude larger than the intensities calculated under the classical assumptions. The enhanced intensity of He I 584 is primarily caused by He II recombination cascades. The enhanced intensity of He II 304 and He II 256 is caused primarily by non-equilibrium helium ionisation. Conclusion: The analysis shows that the long standing problem of the high helium EUV line intensities disappears when taking into account optically thick radiative transfer and non-equilibrium ionisation effects.

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.

Context. The He i 1083 nm is a powerful diagnostic for inferring properties of the upper solar chromosphere, in particular for the magnetic field. The basic formation of the line in one-dimensional models is well understood, but the influence of the complex 3D structure of the chromosphere and corona has however never been investigated. This structure must play an essential role because images taken in He i 1083 nm show structures with widths down to 100 km. Aims. To understand the effect of the three-dimensional temperature and density structure in the solar atmosphere on the formation of the He i 1083 nm line. Methods. We solve the non-LTE radiative transfer problem assuming statistical equilibrium for a simple 9-level helium atom that nevertheless captures all essential physics. As a model atmosphere we use a snapshot from a 3D radiation-MHD simulation computed with the Bifrost code. Ionising radiation from the corona is self-consistently taken into account. Results. The emergent intensity in the He i 1083 nm is set by the source function and the opacity in the upper chromosphere. The former is dominated by scattering of photospheric radiation and does not vary much with spatial location. The latter is determined by the photonionisation rate in the He i ground state continuum, as well as the electron density in the chromosphere. The spatial variation of the flux of ionising radiation is caused by the spatially-structured emissivity of the ionising photons from material at T = 100 kK in the transition region. The hotter coronal material produces more ionising photons, but the resulting radiation field is smooth and does not lead to small-scale variation of the UV flux. The corrugation of the transition region further increases the spatial variation of the amount of UV radiation in the chromosphere.

Using radiation magnetohydrodynamic simulations of the solar atmospheric layers from the upper convection zone to the lower corona, we investigate the self-consistent excitation of slow magneto-acoustic body waves (slow modes) in a magnetic flux concentration. We find that the convective downdrafts in the close surroundings of a two-dimensional flux slab "pump" the plasma inside it in the downward direction. This action produces a downflow inside the flux slab, which encompasses ever higher layers, causing an upwardly propagating rarefaction wave. The slow mode, excited by the adiabatic compression of the downflow near the optical surface, travels along the magnetic field in the upward direction at the tube speed. It develops into a shock wave at chromospheric heights, where it dissipates, lifts the transition region, and produces an offspring in the form of a compressive wave that propagates further into the corona. In the wake of downflows and propagating shock waves, the atmosphere inside the flux slab in the chromosphere and higher tends to oscillate with a period of $\nu\approx 4$~mHz. We conclude that this process of "magnetic pumping" is a most plausible mechanism for the direct generation of longitudinal chromospheric and coronal compressive waves within magnetic flux concentrations, and it may provide an important heat source in the chromosphere. It may also be responsible for certain types of dynamic fibrils.

Solar spicules are chromospheric fibrils that appear everywhere on the Sun, yet their origin is not understood. Using high resolution observations of spicules obtained with the Swedish 1-m Solar Telescope, we aim to understand how spicules appear in filtergrams and Dopplergrams, how they compare in Ca II H and H-alpha, and what can make them appear and disappear. We find thatspicules display a rich and detailed spatial structure, and show a distribution of transverse velocities that when aligned with the line of sight can make them appear at different H-alpha wing positions. They become more abundant at positions closer to the line core, reflecting a distribution of Doppler shifts and widths. In H-alpha width maps they stand out as bright features both on disk and off limb, reflecting their large Doppler motions and possibly higher temperatures than in the typical H-alpha formation region. Spicule lifetimes measured from narrowband images at only a few positions will be an underestimate because Doppler shifts can make them disappear prematurely from such images; for such cases, width maps are a more robust tool. In H-alpha and Ca II H filtergrams, off limb spicules essentially have the same properties, appearance, and evolution. We find that the sudden appearance of spicules can be explained by Doppler shifts from their transverse motions, and does not require other convoluted explanations.

This document was created by the Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) in preparation of the first regular observations of the Sun with the Atacama Large Millimeter/submillimeter Array (ALMA), which are anticipated to start in ALMA Cycle 4 in October 2016. The science cases presented here demonstrate that a large number of scientifically highly interesting observations could be made already with the still limited solar observing modes foreseen for Cycle 4 and that ALMA has the potential to make important contributions to answering long-standing scientific questions in solar physics. With the proposal deadline for ALMA Cycle 4 in April 2016 and the Commissioning and Science Verification campaign in December 2015 in sight, several of the SSALMON Expert Teams composed strategic documents in which they outlined potential solar observations that could be feasible given the anticipated technical capabilities in Cycle 4. These documents have been combined and supplemented with an analysis, resulting in recommendations for solar observing with ALMA in Cycle 4. In addition, the detailed science cases also demonstrate the scientific priorities of the solar physics community and which capabilities are wanted for the next observing cycles. The work on this White Paper effort was coordinated in close cooperation with the two international solar ALMA development studies led by T. Bastian (NRAO, USA) and R. Brajsa, (ESO). This document will be further updated until the beginning of Cycle 4 in October 2016. In particular, we plan to adjust the technical capabilities of the solar observing modes once finally decided and to further demonstrate the feasibility and scientific potential of the included science cases by means of numerical simulations of the solar atmosphere and corresponding simulated ALMA observations.

We study the ascent of granular-sized magnetic bubbles from the solar photosphere through the chromosphere into the transition region and above, for the first time. Such events occurred in a flux emerging region in NOAA 11850 on September 25, 2013. During that time, the first co-observing campaign between the Swedish 1-m Solar Telescope and the IRIS spacecraft was carried out. Simultaneous observations of the chromospheric H$\alpha$ 656.28 nm and \ion{Ca}{2} 854.2 nm lines, plus the photospheric \ion{Fe}{1} 630.25 nm line, were made with the CRISP spectropolarimeter at the SST reaching a spatial resolution of 0."14. At the same time, IRIS was performing a four-step dense raster of the said emerging flux region, taking slit-jaw images at 133 (C~{\sc ii}, transition region), 140 (\ion{Si}{4}, transition region), 279.6 (\ion{Mg}{2} k, core, upper chromosphere), and 283.2 nm (\ion{Mg}{2} k, wing, photosphere). Spectroscopy of several lines was performed by the IRIS spectrograph in the far and near ultraviolet, of which we have used the \ion{Si}{4} 140.3 and the \ion{Mg}{2} k 279.6 nm lines. Coronal images from the Atmospheric Imaging Assembly of the Solar Dynamics Observatory were used to investigate the possible coronal signatures of the flux emergence events. The photospheric and chromospheric properties of small-scale emerging magnetic bubbles have been described in detail in Ortiz et al. (2014). Here we are able to follow such structures up to the transition region. We describe the properties, including temporal delays, of the observed flux emergence in all layers. We believe this may be an important mechanism of transporting energy and magnetic flux from subsurface layers to the transition region and corona.

Active regions (AR) appearing on the surface of the Sun are classified into $\alpha$, $\beta$, $\gamma$, and $\delta$ by the rules of the Mount Wilson Observatory, California on the basis of their topological complexity. Amongst these, the $\delta$-sunspots are known to be super-active and produce the most X-ray flares. Here, we present results from a simulation of the Sun by mimicking the upper layers and the corona, but starting at a more primitive stage than any earlier treatment. We find that this initial state consisting of only a thin sub-photospheric magnetic sheet breaks into multiple flux-tubes which evolve into a colliding-merging system of spots of opposite polarity upon surface emergence, similar to those often seen on the Sun. The simulation goes on to produce many exotic $\delta$-sunspot associated phenomena: repeated flaring in the range of typical solar flare energy release and ejective helical flux ropes with embedded cool-dense plasma filaments resembling solar coronal mass ejections.

The ionization state of the gas in the dynamic solar chromosphere can depart strongly from the instantaneous statistical equilibrium commonly assumed in numerical modeling. We improve on earlier simulations of the solar atmosphere that only included non-equilbrium hydrogen ionization by performing a 2D radiation-magneto-hydrodynamics simulation featuring non-equilibrium ionization of both hydrogen and helium. The simulation includes the effect of hydrogen Lyman-$\alpha$ and the EUV radiation from the corona on the ionization and heating of the atmosphere. Details on code implementation are given. We obtain helium ion fractions that are far from their equilibrium values. Comparison with models with LTE ionization shows that non-equilibrium helium ionization leads to higher temperatures in wave fronts and lower temperatures in the gas between shocks. Assuming LTE ionization results in a thermostat-like behaviour with matter accumulating around the temperatures where the LTE ionization fractions change rapidly. Comparison of DEM curves computed from our models shows that non-equilibrium ionization leads to more radiating material in the temperature range 11-18 kK compared to models with LTE helium ionization. We conclude that non-equilbrium helium ionization is important for the dynamics and thermal structure of the upper chromosphere and transition region. It might also help resolve the problem that intensities of chromospheric lines computed from current models are smaller than those observed.

The C II 133.5 nm multiplet has been observed by NASA's Interface Region Imaging Spectrograph (IRIS) in unprecedented spatial resolution. The aims of this work are to characterize these new observations of the C II lines, place them in context with previous work, and to identify any additional value the C II lines bring when compared with other spectral lines. We make use of wide, long exposure IRIS rasters covering the quiet Sun and an active region. Line properties such as velocity shift and width are extracted from individual spectra and analyzed. The lines have a variety of shapes (mostly single-peak or double-peak), are strongest in active regions and weaker in the quiet Sun. The ratio between the 133.4 nm and 133.5 nm components is always less than 1.8, indicating that their radiation is optically thick in all locations. Maps of the C II line widths are a powerful new diagnostic of chromospheric structures, and their line shifts are a robust velocity diagnostic. Compared with earlier quiet Sun observations, we find similar absolute intensities and mean line widths, but smaller red shifts; this difference can perhaps be attributed to differences in spectral resolution and spatial coverage. The C II intensity maps are somewhat similar to those of transition region lines, but also share some features with chromospheric maps such as those from the Mg II k line, indicating that they are formed between the upper chromosphere and transition region. C II intensity, width, and velocity maps can therefore be used to gather additional information about the upper chromosphere.

Context. The solar chromosphere is the interface between the solar surface and the solar corona. Modelling of this region is difficult because it represents the transition from optically thick to thin radiation escape, from gas-pressure domination to magnetic-pressure domination, from a neutral to an ionised state, from MHD to plasma physics, and from near-equilibrium (LTE) to non-equilibrium conditions. Aims. Our aim is to provide the community with realistic simulations of the magnetic solar outer atmosphere. This will enable detailed comparison of existing and upcoming observations with synthetic observables from the simulations, thereby elucidating the complex interactions of magnetic fields and plasma that are crucial for our understanding of the dynamic outer atmosphere. Methods. We used the radiation magnetohydrodynamics code Bifrost to perform simulations of a computational volume with a magnetic field topology similar to an enhanced network area on the Sun. Results. The full simulation cubes are made available online. The general properties of the simulation are discussed, and limitations are discussed.

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

The O I 135.56 nm line is covered by NASA's Interface Region Imaging Spectrograph (IRIS) small explorer mission which studies how the solar atmosphere is energized. We here study the formation and diagnostic potential of this line by means of non-LTE modelling employing both 1D semi-empirical and 3D radiation-Magneto Hydrodynamic (RMHD) models. We study the basic formation mechanisms and derive a quintessential model atom that incorporates the essential atomic physics for the formation of the O I 135.56 nm line. This atomic model has 16 levels and describes recombination cascades through highly excited levels by effective recombination rates. The ionization balance O I/O II is set by the hydrogen ionization balance through charge exchange reactions. The emission in the O I 135.56 nm line is dominated by a recombination cascade and the line is optically thin. The Doppler shift of the maximum emission correlates strongly with the vertical velocity in its line forming region, which is typically located at 1.0 - 1.5 Mm height. The total intensity of the line emission is correlated with the square of the electron density. Since the O I 135.56 nm line is optically thin, the width of the emission line is a very good diagnostic of non-thermal velocities. We conclude that the O I 135.56 nm line is an excellent probe of the middle chromosphere, and compliments other powerful chromospheric diagnostics of IRIS such as the Mg II h & k lines and the C II lines around 133.5 nm.

Penumbral microjets are short-lived, fine-structured and bright jets that are generally observed in chromospheric imaging of the penumbra of sunspots. Here we investigate their potential transition region signature, by combining observations with the Swedish 1-m Solar Telescope (SST) in the Ca II H and Ca II 8542{\AA} lines with ultraviolet imaging and spectroscopy obtained with the Interface Region Imaging Spectrograph (IRIS), which includes the C II 1334/1335{\AA}, Si IV 1394/1403{\AA} and Mg II h & k 2803/2796{\AA} lines. We find a clear corresponding signal in the IRIS Mg II k, C II and Si IV slit-jaw images, typically offset spatially from the Ca II signature in the direction along the jets: from base to top, the penumbral microjets are predominantly visible in Ca II, Mg II k and C II/Si IV, suggesting progressive heating to transition region temperatures along the jet extent. Hence, these results support the suggestion from earlier studies that penumbral microjets may heat to transition region temperatures.

Advanced 3D radiative MHD simulations now reproduce many properties of the outer solar atmosphere. When including a domain from the convection zone into the corona, a hot chromosphere and corona are self-consistently maintained. Here we study two realistic models, with different simulated area, magnetic field strength and topology, and numerical resolution. These are compared in order to characterize the heating in the 3D-MHD simulations which self-consistently maintains the structure of the atmosphere. We analyze the heating at both large and small scales and find that heating is episodic and highly structured in space, but occurs along loop shaped structures, and moves along with the magnetic field. On large scales we find that the heating per particle is maximal near the transition region and that widely distributed opposite-polarity field in the photosphere leads to a greater heating scale height in the corona. On smaller scales, heating is concentrated in current sheets, the thicknesses of which are set by the numerical resolution. Some current sheets fragment in time, this process occurring more readily in the higher-resolution model leading to spatially highly intermittent heating. The large scale heating structures are found to fade in less than about five minutes, while the smaller, local, heating shows time scales of the order of 2 minutes in one model and 1 minutes in the other, higher-resolution, model.

Solar flares involve complex processes that are coupled and span a wide range of temporal, spatial, and energy scales. Modeling such processes self-consistently has been a challenge in the past. Here we present results from simulations that couple particle kinetics with hydrodynamics of the atmospheric plasma. We combine the Stanford unified Fokker-Planck code that models particle acceleration and transport with the RADYN hydrodynamic code that models the atmospheric response to collisional heating by accelerated electrons through detailed radiative transfer calculations. We perform simulations using two different electron spectra, one an {\it ad hoc} power law and the other predicted by the model of stochastic acceleration by turbulence or plasma waves. Surprisingly, the later model, even with energy flux $\ll 10^{10}$ erg s$^{-1}$ cm$^{-2}$, can cause "explosive" chromospheric evaporation and drive stronger up- and downflows (and hydrodynamic shocks). This is partly because our acceleration model, like many others, produces a spectrum consisting of a quasi-thermal component plus a power-law tail. We synthesize emission line profiles covering different heights in the lower atmosphere, including H$\alpha$ 6563 \AA, HeII 304 \AA, CaII K 3934 \AA\ and SiIV 1393 \AA. One interesting result is the unusual high temperature (up to a few $10^5$ K) of the formation site of HeII 304 \AA, which is expected due to photoionization-recombination under flare conditions, compared to those in the quiet Sun dominated by collisional excitation. When compared with observations, our results can constrain the properties of non-thermal electrons and thus the poorly understood particle acceleration mechanism.

ALMA will open a new chapter in the study of the Sun by providing a leap in spatial resolution and sensitivity compared to currently available mm wave- length observations. In preparation of ALMA, we have carried out a large number of observational tests and state-of-the-art radiation MHD simulations. Here we review the best available observations of the Sun at millimeter wavelengths.Using state of the art radiation MHD simulations of the solar atmosphere we demonstrate the huge potential of ALMA observations for uncovering the nature of the solar chromosphere. We show that ALMA will not only provide a reliable probe of the thermal structure and dynamics of the chromosphere, it will also open up a powerful new diagnostic of magnetic field at chromospheric heights, a fundamentally important, but so far poorly known parameter.

We analyze observations from the Interface Region Imaging Spectrograph of the Mg II k line, the Mg II UV subordinate lines, and the O I 135.6 nm line to better understand the solar plage chromosphere. We also make comparisons with observations from the Swedish 1 m Solar Telescope of the H{\alpha} line, the Ca II 8542 line, and Solar Dynamics Observatory/Atmospheric Imaging Assembly observations of the coronal 19.3 nm line. To understand the observed Mg II profiles, we compare these observations to the results of numerical experiments. The single-peaked or flat-topped Mg II k profiles found in plage imply a transition region at a high column mass and a hot and dense chromosphere of about 6500 K. This scenario is supported by the observed large-scale correlation between moss brightness and filled-in profiles with very little or absent self-reversal. The large wing width found in plage also implies a hot and dense chromosphere with a steep chromospheric temperature rise. The absence of emission in the Mg II subordinate lines constrain the chromospheric temperature and the height of the temperature rise while the width of the O I 135.6 nm line sets a limit to the non-thermal velocities to around 7 km/s.

The Mg II h&k doublet are two of the primary spectral lines observed by the Sun-pointing Interface Region Imaging Spectrograph (IRIS). These lines are tracers of the magnetic and thermal environment that spans from the photosphere to the upper chromosphere. We use a double gaussian model to fit the Mg II h profile for a full-Sun mosaic dataset taken 24-Aug-2014. We use the ensemble of high-quality profile fits to conduct a statistical study on the variability of the line profile as it relates the magnetic structure, dynamics, and center-to-limb viewing angle. The average internetwork profile contains a deeply reversed core and is weakly asymmetric at h2. In the internetwork, we find a strong correlation between h3 wavelength and profile asymmetry as well h1 width and h2 width. The average reversal depth of the h3 core is inversely related to the magnetic field. Plage and sunspots exhibit many profiles which do not contain a reversal. These profiles also occur infrequently in the internetwork. We see indications of magnetically aligned structures in plage and network in statistics associated with the line core, but these structures are not clear or extended in the internetwork. The center-to-limb variations are compared with predictions of semi-empirical model atmospheres. We measure a pronounced limb darkening in the line core which is not predicted by the model. The aim of this work is to provide a comprehensive measurement baseline and preliminary analysis on the observed structure and formation of the Mg II profiles observed by IRIS.

We use 3D radiation magnetohydrodynamic models to investigate how the thermodynamic quantities in the simulation are encoded in observable quantities, thus exploring the diagnostic potential of the 133.5 nm lines. We find that the line core intensity is correlated with the temperature at the formation height but the correlation is rather weak, especially when the lines are strong. The line core Doppler shift is a good measure of the line-of-sight velocity at the formation height. The line width is both dependent on the width of the absorption profile (thermal and non-thermal width) and an opacity broadening factor of 1.2-4 due to the optically thick line formation with a larger broadening for double peak profiles. The 133.5 nm lines can be formed both higher and lower than the core of the Mg II k line depending on the amount of plasma in the 14-50 kK temperature range. More plasma in this temperature range gives a higher 133.5 nm formation height relative to the Mg II k line core. The synthetic line profiles have been compared with IRIS observations. The derived parameters from the simulated line profiles cover the parameter range seen in observations but on average the synthetic profiles are too narrow. We interpret this discrepancy as a combination of a lack of plasma at chromospheric temperatures in the simulation box and too small non-thermal velocities. The large differences in the distribution of properties between the synthetic profiles and the observed ones show that the 133.5 nm lines are powerful diagnostics of the upper chromosphere and lower transition region.

The 133.5 nm lines are important observables for the NASA/SMEX mission Interface Region Imaging Spectrograph (IRIS). To make 3D non-LTE radiative transfer computationally feasible it is crucial to have a model atom with as few levels as possible while retaining the main physical processes. We here develop such a model atom and we study the general formation properties of the C II lines. We find that a nine-level model atom of C I-C III with the transitions treated assuming complete frequency redistribution (CRD) suffices to describe the 133.5 nm lines. 3D scattering effects are important for the intensity in the core of the line. The lines are formed in the optically thick regime. The core intensity is formed in layers where the temperature is about 10kK at the base of the transition region. The lines are 1.2-4 times wider than the atomic absorption profile due to the formation in the optically thick regime. The smaller opacity broadening happens for single peak intensity profiles where the chromospheric temperature is low with a steep source function increase into the transition region, the larger broadening happens when there is a temperature increase from the photosphere to the low chromosphere leading to a local source function maximum and a double peak intensity profile with a central reversal. Assuming optically thin formation with the standard coronal approximation leads to several errors: Neglecting photoionization severly underestimates the amount of C II at temperatures below 16kK, erroneously shifts the formation from 10kK to 25kK and leads to too low intensities.

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.

Ellerman bombs are transient brightenings of the extended wings of the solar Balmer lines in emerging active regions. We describe their properties in the ultraviolet lines sampled by the Interface Region Imaging Spectrograph (IRIS), using simultaneous imaging spectroscopy in H$\alpha$ with the Swedish 1-m Solar Telescope (SST) and ultraviolet images from the Solar Dynamics Observatory for Ellerman bomb detection and identification. We select multiple co-observed Ellerman bombs for detailed analysis. The IRIS spectra strengthen the view that Ellerman bombs mark reconnection between bipolar kilogauss fluxtubes with the reconnection and the resulting bi-directional jet located within the solar photosphere and shielded by overlying chromospheric fibrils in the cores of strong lines. The spectra suggest that the reconnecting photospheric gas underneath is heated sufficiently to momentarily reach stages of ionization normally assigned to the transition region and the corona. We also analyze similar outburst phenomena that we classify as small flaring arch filaments and ascribe to higher-located reconnection. They have different morphology and produce hot arches in million-Kelvin diagnostics.

(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.

A triplet of subordinate lines of Mg II exists in the region around the h&k lines. In solar spectra these lines are seen mostly in absorption, but in some cases can become emission lines. The aim of this work is to study the formation of this triplet, and investigate any diagnostic value they can bring. Using 3D radiative magnetohydrodynamic simulations of quiet Sun and flaring flux emergence, we synthesize spectra and investigate how spectral features respond to the underlying atmosphere. We find that emission in the lines is rare and is typically caused by a steep temperature increase in the lower chromosphere (above 1500 K, with electron densities above 10$^{18}$ m$^{-3}$). In both simulations the lines are sensitive to temperature increases taking place at column masses >= 5e-4 g cm$^{-2}$. Additional information can also be inferred from the peak-to-wing ratio and shape of the line profiles. Using observations from NASA's Interface Region Imaging Spectrograph we find both absorption and emission line profiles with similar shapes to the synthetic spectra, which suggests that these lines represent a useful diagnostic that complements the MgII h&k lines.

The energy for the coronal heating must be provided from the convection zone. The amount and the method by which this energy is transferred into the corona depends on the properties of the lower atmosphere and the corona itself. We review: 1) how the energy could be built in the lower solar atmosphere; 2) how this energy is transferred through the solar atmosphere; and 3) how the energy is finally dissipated in the chromosphere and/or corona. Any mechanism of energy transport has to deal with the various physical processes in the lower atmosphere. We will focus on a physical process that seems to be highly important in the chromosphere and not deeply studied until recently: the ion-neutral interaction effects (INIE) in the chromosphere. We review the relevance and the role of the partial ionization in the chromosphere and show that this process actually impacts considerably the outer solar atmosphere. We include analysis of our 2.5D radiative MHD simulations with the Bifrost code (Gudiksen et al. 2011) including the partial ionization effects on the chromosphere and corona and thermal conduction along magnetic field lines. The photosphere, chromosphere and transition region are partially ionized and the interaction between ionized particles and neutral particles has important consequences on the magneto-thermodynamics of these layers. The INIE are treated using generalized Ohm's law, i.e., we consider the Hall term and the ambipolar diffusion in the induction equation. The interaction between the different species affects the modeled atmosphere as follows: 1) the ambipolar diffusion dissipates magnetic energy and increases the minimum temperature in the chromosphere; 2) the upper chromosphere may get heated and expanded over a greater range of heights. These processes reveal appreciable differences between the modeled atmospheres of simulations with and without INIE.

Spectroscopic observations of prominence eruptions associated with coronal mass ejections (CMEs), although relatively rare, can provide valuable plasma and 3D geometry diagnostics. We report the first observations by the Interface Region Imaging Spectrograph (IRIS) mission of a spectacular fast CME/prominence eruption associated with an equivalent X1.6 flare on 2014 May 9. The maximum plane-of-sky and Doppler velocities of the eruption are 1200 and 460 km/s, respectively. There are two eruption components separated by ~200 km/s in Doppler velocity: a primary, bright component and a secondary, faint component, suggesting a hollow, rather than solid, cone-shaped distribution of material. The eruption involves a left-handed helical structure undergoing counter-clockwise (viewed top-down) unwinding motion. There is a temporal evolution from upward eruption to downward fallback with less-than-free-fall speeds and decreasing nonthermal line widths. We find a wide range of Mg II k/h line intensity ratios (less than ~2 expected for optically-thin thermal emission): the lowest ever-reported median value of 1.17 found in the fallback material and a comparably high value of 1.63 in nearby coronal rain and intermediate values of 1.53 and 1.41 in the two eruption components. The fallback material exhibits a strong ($> 5 \sigma$) linear correlation between the k/h ratio and the Doppler velocity as well as the line intensity. We demonstrate that Doppler dimming of scattered chromospheric emission by the erupted material can potentially explain such characteristics.

The Interface Region Imaging Spectrograph (IRIS) reveals small-scale rapid brightenings in the form of bright grains all over coronal holes and the quiet sun. These bright grains are seen with the IRIS 1330 \AA, 1400 \AA\ and 2796 \AA\ slit-jaw filters. We combine coordinated observations with IRIS and from the ground with the Swedish 1-m Solar Telescope (SST) which allows us to have chromospheric (Ca II 8542 \AA, Ca II H 3968 \AA, H\alpha, and Mg II k 2796 \AA), and transition region (C II 1334 \AA, Si IV 1402) spectral imaging, and single-wavelength Stokes maps in Fe I 6302 \AA at high spatial (0.33"), temporal and spectral resolution. We conclude that the IRIS slit-jaw grains are the counterpart of so-called acoustic grains, i.e., resulting from chromospheric acoustic waves in a non-magnetic environment. We compare slit-jaw images with spectra from the IRIS spectrograph. We conclude that the grain intensity in the 2796 \AA\ slit-jaw filter comes from both the Mg II k core and wings. The signal in the C II and Si IV lines is too weak to explain the presence of grains in the 1300 and 1400 \AA\ slit-jaw images and we conclude that the grain signal in these passbands comes mostly from the continuum. Even though weak, the characteristic shock signatures of acoustic grains can often be detected in IRIS C II spectra. For some grains, spectral signature can be found in IRIS Si IV. This suggests that upward propagating acoustic waves sometimes reach all the way up to the transition region.

Observations of the solar chromosphere in the line-core of the \Halpha\ line show dark elongated structures called fibrils that show swaying motion. We performed a 3D radiation-MHD simulation of a network region, and computed synthetic \Halpha\ images from this simulation to investigate the relation between fibrils and the magnetic field lines in the chromosphere. The periods, amplitudes and phase-speeds of the synthetic fibrils are consistent with those observed. We analyse the relation between the synthetic fibrils and the field lines threading through them, and find that some fibrils trace out the same field line along the fibril's length, but there are also fibrils that sample different field lines at different locations along their length. Fibrils sample the same field lines on a time scale of $\sim200$~s. This is shorter than their own lifetime. We analysed the evolution of the atmosphere along a number of field lines that thread through fibrils and find that they carry slow-mode waves that load mass into the field line, as well as transverse waves that propagate with the Alfv\'en speed. Transverse waves propagating in opposite directions cause an interference pattern with complex apparent phase speeds. The relationship between fibrils and field lines is complex. It is governed by constant migration and swaying of the field lines, their mass loading by slow modes and subsequent draining, and their actual visibility in \Halpha. Field lines are visible where they lie close to the optical depth unity surface. The location of the latter is governed by the height at which the column mass in the chromosphere reaches a certain value. We conclude that using the swaying motion of fibrils as a tracer of chromospheric transverse oscillations must be done with caution.

Probing the magnetism of the upper solar chromosphere requires measuring and modeling the scattering polarization produced by anisotropic radiation pumping in UV spectral lines. Here we apply PORTA (a novel radiative transfer code) to investigate the hydrogen Ly$\alpha$ line in a 3D model of the solar atmosphere resulting from a state of the art MHD simulation. At full spatial resolution the linear polarization signals are very significant all over the solar disk, with a large fraction of the field of view showing line-center amplitudes well above the 1% level. Via the Hanle effect the line-center polarization signals are sensitive to the magnetic field of the model's transition region, even when its mean field strength is only 15 G. The breaking of the axial symmetry of the radiation field produces significant forward-scattering polarization in Ly$\alpha$, without the need of an inclined magnetic field. Interestingly, the Hanle effect tends to decrease such forward-scattering polarization signals in most of the points of the field of view. When the spatial resolution is degraded, the line-center polarization of Ly$\alpha$ drops below the 1% level, reaching values similar to those previously found in 1D semi-empirical models (i.e., up to about 0.5%). The center to limb variation of the spatially-averaged polarization signals is qualitatively similar to that found in 1D models, with the largest line-center amplitudes at $\mu=\cos\theta\approx 0.4$ ($\theta$ being the heliocentric angle). These results are important, both for designing the needed space-based instrumentation and for a reliable interpretation of future observations of the Ly$\alpha$ polarization.

Aims. We use advanced 3D NLTE radiative magnetohydrodynamic simulations of the solar atmosphere to carry out detailed tests of chromospheric diagnostics at millimeter and submillimeter wavelengths. Methods. We focused on the diagnostics of the thermal structure of the chromosphere in the wavelength bands from 0.4 mm up to 9.6 mm that can be accessed with the Atacama Large Millimeter/Submillimeter Array (ALMA) and investigated how these diagnostics are affected by the instrumental resolution. Results. We find that the formation height range of the millimeter radiation depends on the location in the simulation domain and is related to the underlying magnetic structure. Nonetheless, the brightness temperature is a reasonable measure of the gas temperature at the effective formation height at a given location on the solar surface. There is considerable scatter in this relationship, but this is significantly reduced when very weak magnetic fields are avoided. Our results indicate that although instrumental smearing reduces the correlation between brightness and temperature, millimeter brightness can still be used to reliably diagnose electron temperature up to a resolution of 1". If the resolution is more degraded, then the value of the diagnostic diminishes rapidly. Conclusions. We conclude that millimeter brightness can image the chromospheric thermal structure at the height at which the radiation is formed. Thus multiwavelength observations with ALMA with a narrow step in wavelength should provide sufficient information for a tomographic imaging of the chromosphere.

We report on observations of recurrent jets by instruments onboard the Interface Region Imaging Spectrograph (IRIS), Solar Dynamics Observatory (SDO) and Hinode spacecrafts. Over a 4-hour period on July 21st 2013, recurrent coronal jets were observed to emanate from NOAA Active Region 11793. FUV spectra probing plasma at transition region temperatures show evidence of oppositely directed flows with components reaching Doppler velocities of +/- 100 km/s. Raster Doppler maps using a Si IV transition region line show all four jets to have helical motion of the same sense. Simultaneous observations of the region by SDO and Hinode show that the jets emanate from a source region comprising a pore embedded in the interior of a supergranule. The parasitic pore has opposite polarity flux compared to the surrounding network field. This leads to a spine-fan magnetic topology in the coronal field that is amenable to jet formation. Time-dependent data-driven simulations are used to investigate the underlying drivers for the jets. These numerical experiments show that the emergence of current-carrying magnetic field in the vicinity of the pore supplies the magnetic twist needed for recurrent helical jet formation.

As the interface between the Sun's photosphere and corona, the chromosphere and transition region play a key role in the formation and acceleration of the solar wind. Observations from the Interface Region Imaging Spectrograph reveal the prevalence of intermittent small-scale jets with speeds of 80-250 km/s from the narrow bright network lanes of this interface region. These jets have lifetimes of 20-80 seconds and widths of 300 km or less. They originate from small-scale bright regions, often preceded by footpoint brightenings and accompanied by transverse waves with ~20 km/s amplitudes. Many jets reach temperatures of at least ~100000 K and constitute an important element of the transition region structures. They are likely an intermittent but persistent source of mass and energy for the solar wind.

We use coordinated observations with the Interface Region Imaging Spectrograph (IRIS) and the Swedish 1-m Solar Telescope (SST) to identify the disk counterpart of type II spicules in upper-chromospheric and transition region (TR) diagnostics. These disk counterparts were earlier identified through short-lived asymmetries in chromospheric spectral lines: rapid blue- or red-shifted excursions (RBEs or RREs). We find clear signatures of RBEs and RREs in Mg II h & k, often with excursions of the central h3 and k3 absorption features in concert with asymmetries in co-temporal and co-spatial H-alpha spectral profiles. We find spectral signatures for RBEs and RREs in C II 1335 and 1336 A and Si IV 1394 and 1403 A spectral lines and interpret this as a sign that type II spicules are heated to at least TR temperatures, supporting other recent work. These C II and Si IV spectral signals are weaker for a smaller network region than for more extended network regions in our data. A number of bright features around extended network regions observed in IRIS slit-jaw imagery SJI 1330 and 1400, recently identified as network jets, can be clearly connected to H-alpha RBEs and/or RREs in our coordinated data. We speculate that at least part of the diffuse halo around network regions in the IRIS SJI 1330 and 1400 images can be attributed to type II spicules with insufficient opacity in the C II and Si IV lines to stand out as single features in these passbands.

The heating of the outer solar atmospheric layers, i.e., the transition region and corona, to high temperatures is a long standing problem in solar (and stellar) physics. Solutions have been hampered by an incomplete understanding of the magnetically controlled structure of these regions. The high spatial and temporal resolution observations with the Interface Region Imaging Spectrograph (IRIS) at the solar limb reveal a plethora of short, low lying loops or loop segments at transition-region temperatures that vary rapidly, on the timescales of minutes. We argue that the existence of these loops solves a long standing observational mystery. At the same time, based on comparison with numerical models, this detection sheds light on a critical piece of the coronal heating puzzle.

The solar chromosphere and transition region (TR) form an interface between the Sun's surface and its hot outer atmosphere. Here most of the non-thermal energy that powers the solar atmosphere is transformed into heat, although the detailed mechanism remains elusive. High-resolution (0.33-arcsec) observations with NASA's Interface Region Imaging Spectrograph (IRIS) reveal a chromosphere and TR that are replete with twist or torsional motions on sub-arcsecond scales, occurring in active regions, quiet Sun regions, and coronal holes alike. We coordinated observations with the Swedish 1-m Solar Telescope (SST) to quantify these twisting motions and their association with rapid heating to at least TR temperatures. This view of the interface region provides insight into what heats the low solar atmosphere.

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

The solar atmosphere was traditionally represented with a simple one-dimensional model. Over the past few decades, this paradigm shifted for the chromosphere and corona that constitute the outer atmosphere, which is now considered a dynamic structured envelope. Recent observations by IRIS (Interface Region Imaging Spectrograph) reveal that it is difficult to determine what is up and down even in the cool 6000-K photosphere just above the solar surface: this region hosts pockets of hot plasma transiently heated to almost 100,000 K. The energy to heat and accelerate the plasma requires a considerable fraction of the energy from flares, the largest solar disruptions. These IRIS observations not only confirm that the photosphere is more complex than conventionally thought, but also provide insight into the energy conversion in the process of magnetic reconnection.

Solar spicules have eluded modelers and observers for decades. Since the discovery of the more energetic type II, spicules have become a heated topic but their contribution to the energy balance of the low solar atmosphere remains unknown. Here we give a first glimpse of what quiet Sun spicules look like when observed with NASA's recently launched Interface Region Imaging Spectrograph (IRIS). Using IRIS spectra and filtergrams that sample the chromosphere and transition region we compare the properties and evolution of spicules as observed in a coordinated campaign with Hinode and the Atmospheric Imaging Assembly. Our IRIS observations allow us to follow the thermal evolution of type II spicules and finally confirm that the fading of Ca II H spicules appears to be caused by rapid heating to higher temperatures. The IRIS spicules do not fade but continue evolving, reaching higher and falling back down after 500-800 s. Ca II H type II spicules are thus the initial stages of violent and hotter events that mostly remain invisible in Ca II H filtergrams. These events have very different properties from type I spicules, which show lower velocities and no fading from chromospheric passbands. The IRIS spectra of spicules show the same signature as their proposed disk counterparts, reinforcing earlier work. Spectroheliograms from spectral rasters also confirm that quiet Sun spicules originate in bushes from the magnetic network. Our results suggest that type II spicules are indeed the site of vigorous heating (to at least transition region temperatures) along extensive parts of the upward moving spicular plasma.

IRIS data allow us to study the solar transition region (TR) with an unprecedented spatial resolution of 0.33 arcsec. On 2013 August 30, we observed bursts of high Doppler shifts suggesting strong supersonic downflows of up to 200 km/s and weaker, slightly slower upflows in the spectral lines Mg II h and k, C II 1336 \AA, Si IV 1394 \AA, and 1403 \AA, that are correlated with brightenings in the slitjaw images (SJIs). The bursty behavior lasts throughout the 2 hr observation, with average burst durations of about 20 s. The locations of these short-lived events appear to be the umbral and penumbral footpoints of EUV loops. Fast apparent downflows are observed along these loops in the SJIs and in AIA, suggesting that the loops are thermally unstable. We interpret the observations as cool material falling from coronal heights, and especially coronal rain produced along the thermally unstable loops, which leads to an increase of intensity at the loop footpoints, probably indicating an increase of density and temperature in the TR. The rain speeds are on the higher end of previously reported speeds for this phenomenon, and possibly higher than the free-fall velocity along the loops. On other observing days, similar bright dots are sometimes aligned into ribbons, resembling small flare ribbons. These observations provide a first insight into small-scale heating events in sunspots in the TR.

We present the first results of sunspot oscillations from observations by the Interface Region Imaging Spectrograph. The strongly nonlinear oscillation is identified in both the slit-jaw images and the spectra of several emission lines formed in the transition region and chromosphere. We first apply a single Gaussian fit to the profiles of the Mgii 2796.35 {\AA}, Cii 1335.71 {\AA}, and Si iv 1393.76 {\AA} lines in the sunspot. The intensity change is about 30%. The Doppler shift oscillation reveals a sawtooth pattern with an amplitude of about 10 km/s in Si iv. In the umbra the Si iv oscillation lags those of Cii and Mgii by about 3 and 12 s, respectively. The line width suddenly increases as the Doppler shift changes from redshift to blueshift. However, we demonstrate that this increase is caused by the superposition of two emission components. We then perform detailed analysis of the line profiles at a few selected locations on the slit. The temporal evolution of the line core is dominated by the following behavior: a rapid excursion to the blue side, accompanied by an intensity increase, followed by a linear decrease of the velocity to the red side. The maximum intensity slightly lags the maximum blueshift in Si iv, whereas the intensity enhancement slightly precedes the maximum blueshift in Mgii. We find a positive correlation between the maximum velocity and deceleration, a result that is consistent with numerical simulations of upward propagating magnetoacoustic shock waves.

The Ca II 854.2 nm spectral line is a common diagnostic of the solar chromosphere. The average line profile shows an asymmetric core, and its bisector shows a characteristic inverse-C shape. The line actually consists of six components with slightly different wavelengths depending on the isotope of calcium. This isotopic splitting of the line has been taken into account in studies of non-solar stars, but never for the Sun. We performed non-LTE radiative transfer computations from three models of the solar atmosphere and show that the asymmetric line-core and inverse C-shape of the bisector of the 854.2 nm line can be explained by isotopic splitting. We confirm this finding by analysing observations and showing that the line asymmetry is present irrespective of conditions in the solar atmosphere. Finally, we show that inversions based on the Ca II 854.2 nm line should take the isotopic splitting into account, otherwise the inferred atmospheres will contain erroneous velocity gradients and temperatures.

Helium ionization plays an important role in the energy balance of the upper chromosphere and transition region. Helium spectral lines are also often used as diagnostics of these regions. We carry out 1D radiation-hydrodynamics simulations of the solar atmosphere and find that the helium ionization is mostly set by photoionization and direct collisional ionization, counteracted by radiative recombination cascades. By introducing an additional recombination rate mimicking the recombination cascades, we construct a simplified 3 level helium model atom consisting of only the ground states. This model atom is suitable for modeling non-equilibrium helium ionization in 3D numerical models. We perform a brief investigation of the formation of the He I 10830 and He II 304 spectral lines. Both lines show non-equilibrium features that are not recovered with statistical equilibrium models, and caution should therefore be exercised when such models are used as a basis in the interpretation of observations.

The Interface Region Imaging Spectrograph (IRIS) small explorer spacecraft provides simultaneous spectra and images of the photosphere, chromosphere, transition region, and corona with 0.33-0.4 arcsec spatial resolution, 2 s temporal resolution and 1 km/s velocity resolution over a field-of-view of up to 175 arcsec x 175 arcsec. IRIS was launched into a Sun-synchronous orbit on 27 June 2013 using a Pegasus-XL rocket and consists of a 19-cm UV telescope that feeds a slit-based dual-bandpass imaging spectrograph. IRIS obtains spectra in passbands from 1332-1358, 1389-1407 and 2783-2834 Angstrom including bright spectral lines formed in the chromosphere (Mg II h 2803 Angstrom and Mg II k 2796 Angstrom) and transition region (C II 1334/1335 Angstrom and Si IV 1394/1403 Angstrom). Slit-jaw images in four different passbands (C II 1330, Si IV 1400, Mg II k 2796 and Mg II wing 2830 Angstrom) can be taken simultaneously with spectral rasters that sample regions up to 130 arcsec x 175 arcsec at a variety of spatial samplings (from 0.33 arcsec and up). IRIS is sensitive to emission from plasma at temperatures between 5000 K and 10 MK and will advance our understanding of the flow of mass and energy through an interface region, formed by the chromosphere and transition region, between the photosphere and corona. This highly structured and dynamic region not only acts as the conduit of all mass and energy feeding into the corona and solar wind, it also requires an order of magnitude more energy to heat than the corona and solar wind combined. The IRIS investigation includes a strong numerical modeling component based on advanced radiative-MHD codes to facilitate interpretation of observations of this complex region. Approximately eight Gbytes of data (after compression) are acquired by IRIS each day and made available for unrestricted use within a few days of the observation.

The Mg II h&k lines are the prime chromospheric diagnostics of NASA's Interface Region Imaging Spectrograph (IRIS). In the previous papers of this series we used a realistic three-dimensional radiative magnetohydrodynamics model to calculate the h&k lines in detail and investigated how their spectral features relate to the underlying atmosphere. In this work, we employ the same approach to investigate how the h&k diagnostics fare when taking into account the finite resolution of IRIS and different noise levels. In addition, we investigate the diagnostic potential of several other photospheric lines and near-continuum regions present in the near-ultraviolet (NUV) window of IRIS and study the formation of the NUV slit-jaw images. We find that the instrumental resolution of IRIS has a small effect on the quality of the h&k diagnostics; the relations between the spectral features and atmospheric properties are mostly unchanged. The peak separation is the most affected diagnostic, but mainly due to limitations of the simulation. The effects of noise start to be noticeable at a signal-to-noise ratio (S/N) of 20, but we show that with noise filtering one can obtain reliable diagnostics at least down to a S/N of 5. The many photospheric lines present in the NUV window provide velocity information for at least eight distinct photospheric heights. Using line-free regions in the h&k far wings we derive good estimates of photospheric temperature for at least three heights. Both of these diagnostics, in particular the latter, can be obtained even at S/Ns as low as 5.

NASA's Interface Region Imaging Spectrograph (IRIS) small explorer mission will study how the solar atmosphere is energized. IRIS contains an imaging spectrograph that covers the Mg II h&k lines as well as a slit-jaw imager centered at Mg II k. Understanding the observations requires forward modeling of Mg II h&k line formation from 3D radiation-MHD models. We compute the vertically emergent h&k intensity from a snapshot of a dynamic 3D radiation-MHD model of the solar atmosphere, and investigate which diagnostic information about the atmosphere is contained in the synthetic line profiles. We find that the Doppler shift of the central line depression correlates strongly with the vertical velocity at optical depth unity, which is typically located less than 200 km below the transition region (TR). By combining the Doppler shifts of the h and the k line we can retrieve the sign of the velocity gradient just below the TR. The intensity in the central line depression is anticorrelated with the formation height, especially in subfields of a few square Mm. This intensity could thus be used to measure the spatial variation of the height of the transition region. The intensity in the line-core emission peaks correlates with the temperature at its formation height, especially for strong emission peaks. The peaks can thus be exploited as a temperature diagnostic. The wavelength difference between the blue and red peaks provides a diagnostic of the velocity gradients in the upper chromosphere. The intensity ratio of the blue and red peaks correlates strongly with the average velocity in the upper chromosphere. We conclude that the Mg II h&k lines are excellent probes of the very upper chromosphere just below the transition region, a height regime that is impossible to probe with other spectral lines.

NASA's Interface Region Imaging Spectrograph (IRIS) space mission will study how the solar atmosphere is energized. IRIS contains an imaging spectrograph that covers the Mg II h&k lines as well as a slit-jaw imager centered at Mg II k. Understanding the observations will require forward modeling of Mg II h&k line formation from 3D radiation-MHD models. This paper is the first in a series where we undertake this forward modeling. We discuss the atomic physics pertinent to h&k line formation, present a quintessential model atom that can be used in radiative transfer computations and discuss the effect of partial redistribution (PRD) and 3D radiative transfer on the emergent line profiles. We conclude that Mg II h&k can be modeled accurately with a 4-level plus continuum Mg II model atom. Ideally radiative transfer computations should be done in 3D including PRD effects. In practice this is currently not possible. A reasonable compromise is to use 1D PRD computations to model the line profile up to and including the central emission peaks, and use 3D transfer assuming complete redistribution to model the central depression.

We performed a 3D radiative MHD simulation of the solar atmosphere. This simulation shows a jet-like feature that shows similarities to the type II spicules observed for the first time with Hinode. Rapid Blueshifted Events (RBEs) on the solar disk are associated with these spicules. Observational results suggest they may contribute significantly in supplying the corona with hot plasma. We perform a detailed comparison of the properties of the simulated jet with those of type II spicules (observed with Hinode) and RBEs (with ground-based instruments). We analyze variety of synthetic emission and absorption lines from the simulations including chromospheric Ca II and Ha to TR and coronal temperatures (10E4 to several 10E6K). We compare their synthetic intensities, line profiles, Doppler shifts, line widths and asymmetries with observations from Hinode/SOT and EIS, SOHO/SUMER, SST and SDO/AIA. Many properties of the synthetic observables resemble the observations, and we describe in detail the physical processes that lead to these observables. Detailed analysis of the synthetic observables provides insight into how observations should be analyzed to derive information about physical variables in such a dynamic event. For example, we find that LOS superposition in the optically thin atmosphere requires the combination of Doppler shifts and spectral line asymmetry to determine the velocity in the jet. Other properties differ from the observations, especially in the chromospheric lines. The mass density of the part of the spicule with a chromospheric temperature is too low to produce significant opacity in chromospheric lines. These and other discrepancies are described in detail, and we discuss which mechanisms and physical processes may need to be included in the MHD simulations to mimic the thermodynamic processes of the chromosphere and corona, in particular to reproduce type II spicules.

The heating of the Sun's chromosphere remains poorly understood. While progress has been made on understanding what drives the quiet Sun internetwork chromosphere, chromospheric heating in strong magnetic field regions continues to present a difficult challenge, mostly because of a lack of observational constraints. We use high-resolution spectropolarimetric data from the Swedish 1-m Solar Telescope to identify the location and spatio-temporal properties of heating in the magnetic chromosphere. In particular, we report the existence of raised-core spectral line profiles in the Ca II 8542 line. These profiles are characterized by the absence of an absorption line core, showing a quasi-flat profile between +/- 0.5 {\AA}, and are abundant close to magnetic bright-points and plage. Comparison with 3D MHD simulations indicates that such profiles occur when the line-of-sight goes through an "elevated temperature canopy" associated with the expansion with height of the magnetic field of flux concentrations. This temperature canopy in the simulations is caused by ohmic dissipation where there are strong magnetic field gradients. The raised-core profiles are thus indicators of locations of increased chromospheric heating. We characterize the location and temporal and spatial properties of such profiles in our observations, thus providing much stricter constraints on theoretical models of chromospheric heating mechanisms than before.

Spicules have been observed on the sun for more than a century, typically in chromospheric lines such as H-alpha and Ca II H. Recent work has shown that so-called 'type II' spicules may have a role in providing mass to the corona and the solar wind. In chromospheric filtergrams these spicules are not seen to fall back down, and they are shorter-lived and more dynamic than the spicules that have been classically reported in ground-based observations. Observations of type II spicules with Hinode show fundamentally different properties from what was previously measured. In earlier work we showed that these dynamic type II spicules are the most common type, a view that was not properly identified by early observations.The aim of this work is to investigate the effects of spatio-temporal resolution in the classical spicule measurements. Making use of Hinode data degraded to match the observing conditions of older ground-based studies, we measure the properties of spicules with a semi-automated algorithm. These results are then compared to measurements using the original Hinode data. We find that degrading the data has a significant effect on the measured properties of spicules. Most importantly, the results from the degraded data agree well with older studies (e.g. mean spicule duration more than 5 minutes, and upward apparent velocities of about 25 km/s). These results illustrate how the combination of spicule superposition, low spatial resolution and cadence affect the measured properties of spicules, and that previous measurements can be misleading.

In order to understand the heating of the solar corona it is crucial to obtain empirical information on the magnetic field in its lower boundary (the transition region). To this end, we need to measure and model the linear polarization produced by scattering processes in strong UV lines, such as the hydrogen Ly$\alpha$ line. The interpretation of the observed Stokes profiles will require taking into account that the outer solar atmosphere is highly structured and dynamic, and that the height of the transition region may well vary from one place in the atmosphere to another. Here we report on the Ly$\alpha$ scattering polarization signals we have calculated in a realistic model of an enhanced network region, resulting from a state-of-the-art radiation MHD simulation. This model is characterized by spatially complex variations of the physical quantities at transition region heights. The results of our investigation lead us to emphasize that scattering processes in the upper solar chromosphere should indeed produce measurable linear polarization in Ly$\alpha$. More importantly, we show that via the Hanle effect the model's magnetic field produces significant changes in the emergent $Q/I$ and $U/I$ profiles. Therefore, we argue that by measuring the polarization signals produced by scattering processes and the Hanle effect in Ly$\alpha$ and contrasting them with those computed in increasingly realistic atmospheric models, we should be able to decipher the magnetic, thermal and dynamic structure of the upper chromosphere and transition region of the Sun.

Understanding the dynamic solar chromosphere is fundamental in solar physics. Spicules are an important feature of the chromosphere, connecting the photosphere to the corona, potentially mediating the transfer of energy and mass. The aim of this work is to study the properties of spicules over different regions of the sun. Our goal is to investigate if there is more than one type of spicules, and how spicules behave in the quiet sun, coronal holes, and active regions. We make use of high-cadence and high-spatial resolution Ca II H observations taken by Hinode/SOT. Making use of a semi-automated detection algorithm, we self-consistently track and measure the properties of 519 spicules over different regions. We find clear evidence of two types of spicules. Type I spicules show a rise and fall and have typical lifetimes of 150-400 s and maximum ascending velocities of 15-40 km/s, while type II spicules have shorter lifetimes of 50-150 s, faster velocities of 30-110 km/s, and are not seen to fall down, but rather fade at around their maximum length. Type II spicules are the most common, seen in quiet sun and coronal holes. Type I spicules are seen mostly in active regions. There are regional differences between quiet sun and coronal hole spicules, likely attributable to the different field configurations. The properties of type II spicules are consistent with published results of Rapid Blueshifted Events (RBEs), supporting the hypothesis that RBEs are their disk counterparts. For type I spicules we find the relations between their properties to be consistent with a magnetoacoustic shock wave driver, and with dynamic fibrils as their disk counterpart. The driver of type II spicules remains unclear from limb observations.

Determining the temperature distribution of coronal plasmas can provide stringent constraints on coronal heating. Current observations with the Extreme ultraviolet Imaging Spectrograph onboard Hinode and the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory provide diagnostics of the emission measure distribution (EMD) of the coronal plasma. Here we test the reliability of temperature diagnostics using 3D radiative MHD simulations. We produce synthetic observables from the models, and apply the Monte Carlo Markov chain EMD diagnostic. By comparing the derived EMDs with the "true" distributions from the model we assess the limitations of the diagnostics, as a function of the plasma parameters and of the signal-to-noise of the data. We find that EMDs derived from EIS synthetic data reproduce some general characteristics of the true distributions, but usually show differences from the true EMDs that are much larger than the estimated uncertainties suggest, especially when structures with significantly different density overlap along the line-of-sight. When using AIA synthetic data the derived EMDs reproduce the true EMDs much less accurately, especially for broad EMDs. The differences between the two instruments are due to the: (1) smaller number of constraints provided by AIA data, (2) broad temperature response function of the AIA channels which provide looser constraints to the temperature distribution. Our results suggest that EMDs derived from current observatories may often show significant discrepancies from the true EMDs, rendering their interpretation fraught with uncertainty. These inherent limitations to the method should be carefully considered when using these distributions to constrain coronal heating.

Spicules are long, thin, highly dynamic features that jut out ubiquitously from the solar limb. They dominate the interface between the chromosphere and corona and may provide significant mass and energy to the corona. We use high-quality observations with the Swedish 1-m Solar Telescope to establish that so-called type II spicules are characterized by the simultaneous action of three different types of motion: 1. field-aligned flows of order 50-100 km/s, 2. swaying motions of order 15-20 km/s, and 3. torsional motions of order 25-30 km/s. The first two modes have been studied in detail before, but not the torsional motions. Our analysis of many near-limb and off-limb spectra and narrow-band images yields strong evidence that most, if not all, type-II spicules undergo large torsional modulation and that these motions, like spicule swaying, represent Alfvenic waves propagating outward at several hundred km/s. The combined action of the different motions explains the similar morphology of spicule bushes in the outer red and blue wings of chromospheric lines, and needs to be taken into account when interpreting Doppler motions to derive estimates for field-aligned flows in spicules and determining the Alfvenic wave energy in the solar atmosphere. Our results also suggest that large torsional motion is an ingredient in the production of type-II spicules and that spicules play an important role in the transport of helicity through the solar atmosphere.

The structure of the solar chromosphere is believed to be governed by magnetic fields, even in quiet-Sun regions that have a relatively weak photospheric field. During the past decade inversion methods have emerged as powerful tools for analyzing the chromosphere of active regions. The applicability of inversions to infer the stratification of the physical conditions in a dynamic 3D solar chromosphere has not yet been studied in detail. This study aims to establish the diagnostic capabilities of non-local thermodynamical equilibrium (NLTE) inversion techniques of Stokes profiles induced by the Zeeman effect in the Ca II 8542 line. We computed the Ca II atomic level populations in a snapshot from a 3D radiation-MHD simulation of the quiet solar atmosphere in non-LTE using the 3D radiative transfer code Multi3d. These populations were used to compute synthetic full-Stokes profiles in the Ca II 8542 line using 1.5D radiative transfer and the inversion code Nicole. The profiles were then spectrally degraded to account for finite filter width and Gaussian noise was added to account for finite photon flux. These profiles were inverted using Nicole and the results were compared with the original model atmosphere. Our NLTE inversions applied to quiet-Sun synthetic observations provide reasonably good estimates of the chromospheric magnetic field, line-of-sight velocities and somewhat less accurate, but still very useful, estimates of the temperature. Three dimensional scattering of photons cause cool pockets in the chromosphere to be invisible in the line profile and consequently they are also not recovered by the inversions. To successfully detect Stokes linear polarization in this quiet snapshot, a noise level below 10^{-3.5} is necessary.

The magnetic field of the quiet-Sun chromosphere remains a mystery for solar physicists. The reduced number of chromospheric lines are intrinsically hard to model and only a few of them are magnetically sensitive. In this work, we use a 3D numerical simulation of the outer layers of the solar atmosphere, to asses the reliability of non-LTE inversions, in this case applied to the Ca II 8542 \AA\ line. We show that NLTE inversions provide realistic estimates of physical quantities from synthetic observations.

Context. The radiative energy balance in the solar chromosphere is dominated by strong spectral lines that are formed out of LTE. It is computationally prohibitive to solve the full equations of radiative transfer and statistical equilibrium in 3D time dependent MHD simulations. Aims. To find simple recipes to compute the radiative energy balance in the dominant lines under solar chromospheric conditions. Methods. We use detailed calculations in time-dependent and 2D MHD snapshots to derive empirical formulae for the radiative cooling and heating. Results. The radiative cooling in neutral hydrogen lines and the Lyman continuum, the H and K and intrared triplet lines of singly ionized calcium and the h and k lines of singly ionized magnesium can be written as a product of an optically thin emission (dependent on temperature), an escape probability (dependent on column mass) and an ionization fraction (dependent on temperature). In the cool pockets of the chromosphere the same transitions contribute to the heating of the gas and similar formulae can be derived for these processes. We finally derive a simple recipe for the radiative heating of the chromosphere from incoming coronal radiation. We compare our recipes with the detailed results and comment on the accuracy and applicability of the recipes.

We use state-of-the-art radiation-MHD simulations and 3D non-LTE radiative transfer computations to investigate \Halpha\ line formation in the solar chromosphere and apply the results of this investigation to develop the potential of \Halpha\ as diagnostic of the chromosphere. We show that one can accurately model \Halpha\ line formation assuming statistical equilibrium and complete frequency redistribution provided the computation of the model atmosphere included non-equilibrium ionization of hydrogen, and the Lyman-$\alpha$ and Lyman-$\beta$ line profiles are described by Doppler profiles. We find that 3D radiative transfer is essential in modeling hydrogen lines due to the low photon destruction probability in \Halpha. The \Halpha\ opacity in the upper chromosphere is mainly sensitive to the mass density and only weakly sensitive to temperature. We find that the \Halpha\ line-core intensity is correlated with the average formation height: the larger the average formation height, the lower the intensity. The line-core width is a measure of the gas temperature in the line-forming region. The fibril-like dark structures seen in \Halpha\ line-core images computed from our model atmosphere are tracing magnetic field lines. These structures are caused by field-aligned ridges of enhanced chromospheric mass density that raise their average formation height, and therefore makes them appear dark against their deeper-formed surroundings. We compare with observations, and find that the simulated line-core widths are very similar to the observed ones, without the need for additional microturbulence.

The tropical wisdom that when it is hot and dense we can expect rain might also apply to the Sun. Indeed, observations and numerical simulations have shown that strong heating at footpoints of loops, as is the case for active regions, puts their coronae out of thermal equilibrium, which can lead to a phenomenon known as catastrophic cooling. Following local pressure loss in the corona, hot plasma locally condenses in these loops and dramatically cools down to chromospheric temperatures. These blobs become bright in H-alpha and Ca II H in time scales of minutes, and their dynamics seem to be subject more to internal pressure changes in the loop rather than to gravity. They thus become trackers of the magnetic field, which results in the spectacular coronal rain that is observed falling down coronal loops. In this work we report on high resolution observations of coronal rain with the Solar Optical Telescope (SOT) on Hinode and CRISP at the Swedish Solar Telescope (SST). A statistical study is performed in which properties such as velocities and accelerations of coronal rain are derived. We show how this phenomenon can constitute a diagnostic tool for the internal physical conditions inside loops. Furthermore, we analyze transverse oscillations of strand-like condensations composing coronal rain falling in a loop, and discuss the possible nature of the wave. This points to the important role that coronal rain can play in the fields of coronal heating and coronal seismology.

We present the results of numerical simulations of wave propagation and jet formation in solar atmosphere models with different magnetic field configurations. The presence in the chromosphere of waves with periods longer than the acoustic cutoff period has been ascribed to either strong inclined magnetic fields, or changes in the radiative relaxation time. Our simulations include a sophisticated treatment of radiative losses, as well as fields with different strengths and inclinations. Using Fourier and wavelet analysis techniques, we investigate the periodicity of the waves that travel through the chromosphere. We find that the velocity signal is dominated by waves with periods around 5 minutes in regions of strong, inclined field, including at the edges of strong flux tubes where the field expands, whereas 3-minute waves dominate in regions of weak or vertically oriented fields. Our results show that the field inclination is very important for long-period wave propagation, whereas variations in the radiative relaxation time have little effect. Furthermore, we find that atmospheric conditions can vary significantly on timescales of a few minutes, meaning that a Fourier analysis of wave propagation can be misleading. Wavelet techniques take variations with time into account and are more suitable analysis tools. Finally, we investigate the properties of jets formed by the propagating waves once they reach the transition region, and find systematic differences between the jets in inclined field regions and those in vertical field regions, in agreement with observations of dynamic fibrils.

Context: Numerical simulations of stellar convection and photospheres have been developed to the point where detailed shapes of observed spectral lines can be explained. Stellar atmospheres are very complex, and very different physical regimes are present in the convection zone, photosphere, chromosphere, transition region and corona. To understand the details of the atmosphere it is necessary to simulate the whole atmosphere since the different layers interact strongly. These physical regimes are very diverse and it takes a highly efficient massively parallel numerical code to solve the associated equations. Aims: The design, implementation and validation of the massively parallel numerical code Bifrost for simulating stellar atmospheres from the convection zone to the corona. Methods: The code is subjected to a number of validation tests, among them the Sod shock tube test, the Orzag-Tang colliding shock test, boundary condition tests and tests of how the code treats magnetic field advection, chromospheric radiation, radiative transfer in an isothermal scattering atmosphere, hydrogen ionization and thermal conduction. Results: Bifrost completes the tests with good results and shows near linear efficiency scaling to thousands of computing cores.

aims: We aim to provide an estimate of the minimum temperature of the quiet solar chromosphere. methods: We perform a 2D radiation-MHD simulation spanning the upper convection zone to the lower corona. The simulation includes non-LTE radiative transfer and {an equation-of-state that includes non-equilibrium ionization of hydrogen and non-equilibrium H_2 molecule formation}. We analyze the reliability of the various assumptions made in our model in order to assess the realism of the simulation. results: Our simulation contains pockets of cool gas with down to 1660 K from 1 Mm up to 3.2 Mm height. It overestimates the radiative heating, and contains non-physical heating below 1660 K. Therefore we conclude that cool pockets in the quiet solar chromosphere might have even lower temperatures than in the simulation, provided that there exist areas in the chromosphere without significant magnetic heating. We suggest off-limb molecular spectroscopy to look for such cool pockets and 3D simulations including a local dynamo and a magnetic carpet to investigate Joule heating in the quiet chromosphere.

Imaging spectroscopy of the solar atmosphere using the NaI D1 line yields marked asymmetry between the blue and red line wings: sampling a quiet-Sun area in the blue wing displays reversed granulation, whereas sampling in the red wing displays normal granulation. The MgI b2 line of comparable strength does not show this asymmetry, nor does the stronger CaII 8542 line. We demonstrate the phenomenon with near-simultaneous spectral images in NaI D1, MgI b2, and CaII 8542 from the Swedish 1-m Solar Telescope. We then explain it with line-formation insights from classical 1D modeling and with a 3D magnetohydrodynamical simulation combined with NLTE spectral line synthesis that permits detailed comparison with the observations in a common format. The cause of the imaging asymmetry is the combination of correlations between intensity and Dopplershift modulation in granular overshoot and the sensitivity to these of the steep profile flanks of the NaI D1 line. The MgI b2 line has similar core formation but much wider wings due to larger opacity buildup and damping in the photosphere. Both lines obtain marked core asymmetry from photospheric shocks in or near strong magnetic concentrations, less from higher-up internetwork shocks that produce similar asymmetry in the spatially averaged CaII 8542 profile.

The measurement of Doppler velocities in spectroscopic solar observations requires a reference for the local frame of rest. The rotational and radial velocities of the Earth and the rotation of the Sun introduce velocity offsets in the observations. Normally, good references for velocities are missing (e.g. telluric lines), especially in filter-based spectropolarimetric observations. We determine an absolute reference for line-of-sight velocities measured from solar observations for any heliocentric angle, calibrating the convective line shift of spatially-averaged profiles on quiet sun from a 3D hydrodynamical simulation. This method works whenever there is quiet sun in the field-of-view, and it has the advantage of being relatively insensitive to uncertainties in the atomic data. We carry out radiative transfer computations in LTE for selected C I and Fe I lines, whereas the Ca II infrared lines are synthesized in non-LTE. The degree of asymmetry and the line shifts, however, show a clear dependence on the heliocentric angle and the properties of the lines. The profiles at \mu = 1 are compared with observed profiles to prove their reliability, and they are tested against errors induced by the LTE calculations, inaccuracies in the atomic data and the 3D simulation. Theoretical quiet-sun profiles of lines commonly used by solar observers are provided to the community. Those can be used as absolute references for line-of-sight velocities. The limb effect is produced by the projection of the 3D atmosphere along the line of sight. We estimate the precision of the disk-center line shifts to be approximately 50 m s^-1, but the off-center profiles remain to be tested against observations.

Our aim is to determine the dominant processes and timescales for the ionisation equilibrium of calcium under solar chromospheric conditions. The study is based on numerical simulations with the RADYN code, which includes hydrodynamics, radiative transfer, and a detailed non-equilibrium treatment of hydrogen, calcium, and helium. The simulations are characterised by upwards propagating shock waves, which cause strong temperature fluctuations and variations of the ionisation degree of calcium. The passage of a hot shock front leads to a strong net ionisation of Ca II, rapidly followed by net recombination. The relaxation timescale of the Ca ionisation state is found to be of the order of a few seconds at the top of the photosphere and 10 to 30 s in the upper chromosphere. Generally, the timescales are significantly reduced in the wakes of hot shock fronts. The timescales can be reliably determined from a simple analysis of the eigenvalues of the transition rate matrix. The timescales are dominated by the radiative recombination from Ca III into the metastable Ca II energy levels of the 4d 2D term. These transitions depend strongly on the density of free electrons and therefore on the (non-equilibrium) ionisation degree of hydrogen, which is the main electron donor. The ionisation/recombination timescales derived here are too long for the assumption of an instantaneous ionisation equilibrium to be valid and, on the other hand, are not long enough to warrant an assumption of a constant ionisation fraction. Fortunately, the ionisation degree of Ca II remains small in the height range, where the cores of the H, K, and the infrared triplet lines are formed. We conclude that the difference due to a detailed treatment of Ca ionisation has only negligible impact on the modelling of spectral lines of Ca II and the plasma properties under the conditions in the quiet solar chromosphere.

We present the implementation of a radiative transfer solver with coherent scattering in the new BIFROST code for radiative magneto-hydrodynamical (MHD) simulations of stellar surface convection. The code is fully parallelized using MPI domain decomposition, which allows for large grid sizes and improved resolution of hydrodynamical structures. We apply the code to simulate the surface granulation in a solar-type star, ignoring magnetic fields, and investigate the importance of coherent scattering for the atmospheric structure. A scattering term is added to the radiative transfer equation, requiring an iterative computation of the radiation field. We use a short-characteristics-based Gauss-Seidel acceleration scheme to compute radiative flux divergences for the energy equation. The effects of coherent scattering are tested by comparing the temperature stratification of three 3D time-dependent hydrodynamical atmosphere models of a solar-type star: without scattering, with continuum scattering only, and with both continuum and line scattering. We show that continuum scattering does not have a significant impact on the photospheric temperature structure for a star like the Sun. Including scattering in line-blanketing, however, leads to a decrease of temperatures by about 350\,K below log tau < -4. The effect is opposite to that of 1D hydrostatic models in radiative equilibrium, where scattering reduces the cooling effect of strong LTE lines in the higher layers of the photosphere. Coherent line scattering also changes the temperature distribution in the high atmosphere, where we observe stronger fluctuations compared to a treatment of lines as true absorbers.

Broad-band images in the Ca II H line, from the BFI instrument on the Hinode spacecraft, show emission from spicules emerging from and visible right down to the observed limb. Surprisingly, little absorption of spicule light is seen along their lengths. We present formal solutions to the transfer equation for given (ad-hoc) source functions, including a stratified chromosphere from which spicules emanate. The model parameters are broadly compatible with earlier studies of spicules. The visibility of Ca II spicules down to the limb in Hinode data seems to require that spicule emission be Doppler shifted relative to the stratified atmosphere, either by supersonic turbulent or organized spicular motion. The non-spicule component of the chromosphere is almost invisible in the broad band BFI data, but we predict that it will be clearly visible in high spectral resolution data. Broad band Ca II H limb images give the false impression that the chromosphere is dominated by spicules. Our analysis serves as a reminder that the absence of a signature can be as significant as its presence.

In our studies of the dynamics and energetics of the solar atmosphere, we have detected, in high-quality observations from Hinode SOT/NFI, ubiquitous small-scale upflows which move horizontally with supersonic velocities in the quiet Sun. We present the properties of these fast moving upflows (FMUs) and discuss different interpretations.

The present study addresses the following questions: How representative of the actual velocities in the solar atmosphere are the Doppler shifts of spectral lines? How reliable is the velocity signal derived from narrowband filtergrams? How well defined is the height of the measured Doppler signal? Why do phase difference spectra always pull to 0 degrees phase lag at high frequencies? Can we actually observe high frequency waves (P<= 70s)? What is the atmospheric MTF of high frequency waves? How reliably can we determine the energy flux of high frequency waves? We address these questions by comparing observations obtained with Hinode/NFI with results from two 3D numerical simulations (Oslo Stagger and CO5BOLD). Our results suggest that the observed high frequency Doppler velocity signal is caused by rapid height variations of the velocity response function in an atmosphere with strong velocity gradients and cannot be interpreted as evidence of propagating high frequency acoustic waves. Estimates of the energy flux of high frequency waves should be treated with caution, in particular those that apply atmospheric MTF corrections.

We present evidence that transition region red-shifts are naturally produced in episodically heated models where the average volumetric heating scale height lies between that of the chromospheric pressure scale height of 200 km and the coronal scale height of 50 Mm. In order to do so we present results from 3d MHD models spanning the upper convection zone up to the corona, 15 Mm above the photosphere. Transition region and coronal heating in these models is due both the stressing of the magnetic field by photospheric and convection `zone dynamics, but also in some models by the injection of emerging magnetic flux.

We have performed 3D radiation MHD simulations extending from the convection zone to the corona covering a box 16 Mm$^3$ at 32 km spatial resolution. The simulations show very fine structure in the chromosphere with acoustic shocks interacting with the magnetic field. Magnetic flux concentrations have a temperature lower than the surroundings in the photosphere but higher in the low chromosphere. The heating is there mostly through ohmic dissipation preferentially at the edges of the flux concentrations. The magnetic field is often wound up around the flux concentrations. When acoustic waves travel up along the field this topology leads to swirling motions seen in chromospheric diagnostic lines such as the calcium infrared triplet.

The NaD1 line in the solar spectrum is sometimes attributed to the solar chromosphere. We study its formation in quiet-Sun network and internetwork. We first present high-resolution profile-resolved images taken in this line with the imaging spectrometer IBIS at the Dunn Solar Telescope and compare these to simultaneous chromospheric images taken in Ca 8542 and Halpha. We then model NaD1 formation by performing 3D NLTE profile synthesis for a snapshot from a 3D radiation-magnetohydrodynamics simulation. We find that most NaD1 brightness is not chromospheric but samples the magnetic concentrations that make up the quiet-Sun network in the photosphere, well below the height where they merge into chromospheric canopies, with aureoles from 3D resonance scattering. The line core is sensitive to magneto-acoustic shocks in and near magnetic concentrations, where shocks occur deeper than elsewhere, and may provide evidence of heating sited deep within magnetic concentrations.

Recently a second type of spicules was discovered at the solar limb with the Solar Optical Telescope (SOT) onboard the Japanese Hinode spacecraft. These previously unrecognized type II spicules are thin chromospheric jets that are shorter-lived (10-60 s) and that show much higher apparent upward velocities (of order 50-100 km/s) than the classical spicules. Since they have been implicated in providing hot plasma to coronal loops, their formation, evolution and properties are important ingredients for a better understanding of the mass and energy balance of the low solar atmosphere. Here we report on the discovery of the disk counterparts of type II spicules using spectral imaging data in the Ca II 854.2 nm and Halpha lines with the CRisp Imaging SpectroPolarimeter (CRISP) at the Swedish Solar Telescope (SST) in La Palma. We find rapid blueward excursions in the line profiles of both chromospheric lines that correspond to thin, jet-like features that show apparent velocities of order 50 km/s. These blueward excursions seem to form a separate absorbing component with Doppler shifts of order 20 and 50 km/s for the Ca II 854.2 nm and Halpha line respectively. We show that the appearance, lifetimes, longitudinal and transverse velocities and occurrence rate of these rapid blue excursions on the disk are very similar to those of the type II spicules at the limb. A detailed study of the spectral line profiles in these events suggests that plasma is accelerated along the jet, and plasma is being heated throughout the short lifetime of the event.

3D simulations of magnetic flux emergence are carried out in a computational domain spanning the upper layers of the convection zone to the lower corona. We use the Oslo Staggered Code to solve the full MHD equations with non-grey and NLTE radiative transfer and thermal conduction along the magnetic field lines. In this paper we concentrate on the later stages of the simulations and study the evolution of the structure of the rising flux in the upper chromosphere and corona, the interaction between the emerging flux and the weak coronal magnetic field initially present, and the associated dynamics. The flux tube injected at the bottom boundary rises to the photosphere where it largely remains. However, some parts of the flux tube become unstable and expand in patches into the upper chromosphere. The flux rapidly expands towards the corona, pushing the coronal and transition region material aside, lifting and maintaining the transition region at heights greater than 5 Mm above the photosphere for extensive periods of time. The pre-existing magnetic field in the corona and transition region is perturbed by the incoming flux and reoriented by a series of high Joule heating events. Low density structures form in the corona while at later times a high density filamentary structure appears in the lower part of the expanding flux. The dynamics of these and other structures is discussed. While Joule heating due to the expanding flux is episodic, it increases in relative strength as fresh magnetic field rises and becomes energetically important in the upper chromosphere and corona at later times. Chromospheric, transition region and coronal lines are computed and their response to the perturbation caused by the expanding emerging flux is discussed.

We analyze features that resemble type i spicules in two different 3D numerical simulations in which we include horizontal magnetic flux emergence in a computational domain spanning the upper layers of the convection zone to the lower corona. The two simulations differ mainly in the preexisting ambient magnetic field strength and in the properties of the inserted flux tube. We use the Oslo Staggered Code (OSC) to solve the full MHD equations with non-grey and non-LTE radiative transfer and thermal conduction along the magnetic field lines. We find a multitude of features that show a spatiotemporal evolution that is similar to that observed in type i spicules, which are characterized by parabolic height vs. time profiles, and are dominated by rapid upward motion at speeds of 10-30 km/s, followed by downward motion at similar velocities. We measured the parameters of the parabolic profile of the spicules and find similar correlations between the parameters as those found in observations. The values for height (or length) and duration of the spicules found in the simulations are more limited in range than those in the observations. The spicules found in the simulation with higher preexisting ambient field have shorter length and smaller velocities. From the simulations, it appears that these kinds of spicules can, in principle, be driven by a variety of mechanisms that include p-modes, collapsing granules, magnetic energy release in the photosphere and lower chromosphere and convective buffeting of flux concentrations.

We analyze the important formation processes for the Mg I 457.1 nm line. This line is an intercom- bination line and the source function is close to the local thermodynamic equilibrium (LTE) value. The strong coupling to the local temperature and the relatively high population of the lower level (the ground state of Mg I) makes this line an ideal candidate for temperature diagnostics in the lower chromosphere/temperature minimum region. Linking the temperature probed to an absolute phys- ical height is non trivial because of Non-LTE ionization. We analyze the Non-LTE effects and find that photo-ionization from the lower energy levels together with cascading collisional recombination dominate the ionization balance. Taking properly into account the line-blanketing in the UV is essen- tial for obtaining the right photoionization rates. The identification of the main Non-LTE effects in the line allows us to construct a quintessential model atom, ideal for computationally demanding tasks, e.g. full 3D and/or time-dependent radiative transfer. Furthermore we analyze the diagnostic potential to temperature of this line in solar-like atmospheres, by synthesizing the line from a series of parametrized atmospheric models. These models have been constructed with fixed effective temper- ature, but with a variable heat term in the energy equation to obtain a chromospheric temperature rise at different heights. We conclude that the line has a significant potential in the diagnostics of the lower chromosphere temperature structure, especially for cooler atmospheres, such as sunspots.

Interpretation of imagery of the solar chromosphere in the widely used \CaIIIR infrared line is hampered by its complex, three-dimensional and non-LTE formation. Forward modelling is required to aid understanding. We use a 3D non-LTE radiative transfer code to compute synthetic \CaIIIR images from a radiation-MHD simulation of the solar atmosphere spanning from the convection zone to the corona. We compare the simulation with observations obtained with the CRISP filter at the Swedish 1--m Solar Telescope. We find that the simulation reproduces dark patches in the blue line wing caused by Doppler shifts, brightenings in the line core caused by upward-propagating shocks and thin dark elongated structures in the line core that form the interface between upward and downward gas motion in the chromosphere. The synthetic line core is narrower than the observed one, indicating that the sun exhibits both more vigorous large-scale dynamics as well as small scale motions that are not resolved within the simulation, presumably owing to a lack of spatial resolution.

We present a detailed study of the non-Local Thermodynamic Equilibrium (non-LTE) formation of the high-excitation neutral oxygen 777 nm triplet in MARCS model atmospheres representative of late-type stars with spectral types F to K. We carried out the calculations using the statistical equilibrium code MULTI, including estimates of the impact on elemental abundance analysis. The atomic model employed includes, in particular, recent quantum-mechanical electron collision data. We confirm that the O I triplet lines form under non-LTE conditions in late-type stars, suffering negative abundance corrections with respect to LTE. At low metallicity, large line opacity stems from triplet-quintet intersystem electron collisions, a form of coupling previously not considered or seriously underestimated. The non-LTE effects become generally severe for models (both giants and dwarfs) with higher T_eff. Interestingly, in metal-poor turn-off stars, the negative non-LTE abundance corrections tend to rapidly become more severe towards lower metallicity. When neglecting H collisions, they amount to as much as ~ 0.9 dex and ~ 1.2 dex, respectively at [Fe/H]=-3 and [Fe/H]=-3.5. Even when such collisions are included, the LTE abundance remains a serious overestimate, correspondingly by ~ 0.5 dex and ~ 0.9 dex at such low metallicities. Although the poorly known inelastic hydrogen collisions thus remain an important uncertainty, the large metallicity-dependent non-LTE effects seem to point to a resulting "low" (compared to LTE) [O/Fe] in metal-poor halo stars.

Recently, De Pontieu et al. (2007b) discovered a class of spicules that evolves more rapidly than previously known spicules, with rapid apparent motions of 50--150 km s${}^{-1}$, thickness of a few 100 km, and lifetimes of order 10--60 seconds. These so-called type II spicules have been difficult to study because of limited spatio-temporal and thermal resolution. Here we use the IBIS instrument to search for the high velocities in the disk counterpart of type II spicules. We have detected rapidly evolving events, with lifetimes that are less than a minute and often equal to the cadence of the instrument (19 secs). These events are characterized by a Doppler shift that only appears in the blue wing of the Ca II IR line. Furthermore the spatial extent, lifetime, and location near network, all suggest a link to type II spicules. However, the magnitude of the measured Doppler velocity is significantly lower than the apparent motions seen at the limb. We use Monte Carlo simulations to show that this discrepancy can be explained by a forward model in which the visibility on the disk of the high-velocity flows in these events is limited by a combination of line-of-sight projection and reduced opacity in upward propelled plasma, especially in reconnection driven jets that are powered by a roughly constant energy supply.

3D numerical simulations of a horizontal magnetic flux tube emergence with different twist are carried out in a computational domain spanning the upper layers of the convection zone to the lower corona. We use the Oslo Staggered Code to solve the full MHD equations with non-grey and non-LTE radiative transfer and thermal conduction along the magnetic field lines. The emergence of the magnetic flux tube input at the bottom boundary into a weakly magnetized atmosphere is presented. The photospheric and chromospheric response is described with magnetograms, synthetic images and velocity field distributions. The emergence of a magnetic flux tube into such an atmosphere results in varied atmospheric responses. In the photosphere the granular size increases when the flux tube approaches from below. In the convective overshoot region some 200km above the photosphere adiabatic expansion produces cooling, darker regions with the structure of granulation cells. We also find collapsed granulation in the boundaries of the rising flux tube. Once the flux tube has crossed the photosphere, bright points related with concentrated magnetic field, vorticity, high vertical velocities and heating by compressed material are found at heights up to 500km above the photosphere. At greater heights in the magnetized chromosphere, the rising flux tube produces a cool, magnetized bubble that tends to expel the usual chromospheric oscillations. In addition the rising flux tube dramatically increases the chromospheric scale height, pushing the transition region and corona aside such that the chromosphere extends up to 6Mm above the photosphere. The emergence of magnetic flux tubes through the photosphere to the lower corona is a relatively slow process, taking of order 1 hour.

We use coordinated Hinode SOT/EIS observations that include high-resolution magnetograms, chromospheric and TR imaging and TR/coronal spectra in a first test to study how the dynamics of the TR are driven by the highly dynamic photospheric magnetic fields and the ubiquitous chromospheric waves. Initial analysis shows that these connections are quite subtle and require a combination of techniques including magnetic field extrapolations, frequency-filtered time-series and comparisons with synthetic chromospheric and TR images from advanced 3D numerical simulations. As a first result, we find signatures of magnetic flux emergence as well as 3 and 5 mHz wave power above regions of enhanced photospheric magnetic field in both chromospheric, transition region and coronal emission.

We use high-resolution observations of the Sun in Ca II H 3968 A from the Solar Optical Telescope on Hinode to show that there are at least two types of spicules that dominate the structure of the magnetic solar chromosphere. Both types are tied to the relentless magnetoconvective driving in the photosphere, but have very different dynamic properties. ``Type-I'' spicules are driven by shock waves that form when global oscillations and convective flows leak into the upper atmosphere along magnetic field lines on 3-7 minute timescales. ``Type-II'' spicules are much more dynamic: they form rapidly (in ~10s), are very thin (<200km wide), have lifetimes of 10-150s (at any one height) and seem to be rapidly heated to (at least) transition region temperatures, sending material through the chromosphere at speeds of order 50-150 km/s. The properties of Type II spicules suggest a formation process that is a consequence of magnetic reconnection, typically in the vicinity of magnetic flux concentrations in plage and network. Both types of spicules are observed to carry Alfven waves with significant amplitudes of order 20 km/s.

We present high spatial resolution spectroscopic measurements of dynamic fibrils (DFs) in the Ca {\small{II}} 8662 {\AA} line. These data show clear Doppler shifts in the identified DFs, which demonstrates that at least a subset of DFs are actual mass motions in the chromosphere. A statistical analysis of 26 DFs reveals a strong and statistically significant correlation between the maximal velocity and the deceleration. The range of the velocities and the decelerations are substantially lower, about a factor two, in our spectroscopic observations compared to the earlier results based on proper motion in narrow band images. There are fundamental differences in the different observational methods; when DFs are observed spectroscopically the measured Doppler shifts are a result of the atmospheric velocity, weighted with the response function to velocity over an extended height. When the proper motion of DFs is observed in narrow band images, the movement of the top of the DF is observed. This point is sharply defined because of the high contrast between the DF and the surroundings. The observational differences between the two methods are examined by several numerical experiments using both numerical simulations and a time series of narrow band H$\alpha$ images. With basis in the simulations we conclude that the lower maximal velocity is explained by the low formation height of the Ca IR line. We conclude that the present observations support the earlier result that DFs are driven by magneto-acoustic shocks exited by convective flows and p-modes.

The ionization of hydrogen in the solar chromosphere and transition region does not obey LTE or instantaneous statistical equilibrium because the timescale is long compared with important hydrodynamical timescales, especially of magneto-acoustic shocks. We implement an algorithm to compute non-equilibrium hydrogen ionization and its coupling into the MHD equations within an existing radiation MHD code, and perform a two-dimensional simulation of the solar atmosphere from the convection zone to the corona. Analysis of the simulation results and comparison to a companion simulation assuming LTE shows that: a) Non-equilibrium computation delivers much smaller variations of the chromospheric hydrogen ionization than for LTE. The ionization is smaller within shocks but subsequently remains high in the cool intershock phases. As a result, the chromospheric temperature variations are much larger than for LTE because in non-equilibrium, hydrogen ionization is a less effective internal energy buffer. The actual shock temperatures are therefore higher and the intershock temperatures lower. b) The chromospheric populations of the hydrogen n = 2 level, which governs the opacity of Halpha, are coupled to the ion populations. They are set by the high temperature in shocks and subsequently remain high in the cool intershock phases. c) The temperature structure and the hydrogen level populations differ much between the chromosphere above photospheric magnetic elements and above quiet internetwork. d) The hydrogen n = 2 population and column density are persistently high in dynamic fibrils, suggesting that these obtain their visibility from being optically thick in Halpha also at low temperature.

We use Hinode/SOT Ca II H-line and blue continuum broadband observations to study the presence and power of high frequency acoustic waves at high spatial resolution. We find that there is no dominant power at small spatial scales; the integrated power using the full resolution of Hinode (0.05'' pixels, 0.16'' resolution) is larger than the power in the data degraded to 0.5'' pixels (TRACE pixel size) by only a factor of 1.2. At 20 mHz the ratio is 1.6. Combining this result with the estimates of the acoustic flux based on TRACE data of Fossum & Carlsson (2006), we conclude that the total energy flux in acoustic waves of frequency 5-40 mHz entering the internetwork chromosphere of the quiet Sun is less than 800 W m$^{-2}$, inadequate to balance the radiative losses in a static chromosphere by a factor of five.

A major goal in solar physics has during the last five decades been to find how energy flux generated in the solar convection zone is transported and dissipated in the outer solar layers. Progress in this field has been slow and painstaking. However, advances in computer hardware and numerical methods, vastly increased observational capabilities and growing physical insight seem finally to be leading towards understanding. Here we present exploratory numerical MHD models that span the entire solar atmosphere from the upper convection zone to the lower corona. These models include non-grey, non-LTE radiative transport in the photosphere and chromosphere, optically thin radiative losses as well as magnetic field-aligned heat conduction in the transition region and corona.

Spectral diagnostic features formed in the solar chromosphere are few and difficult to interpret -- they are neither formed in the optically thin regime nor in local thermodynamic equilibrium (LTE). To probe the state of the chromosphere, both from observations and theory, it is therefore necessary with modeling. I discuss both traditional semi-empirical modeling, numerical experiments illustrating important ingredients necessary for a self-consistent theoretical modeling of the solar chromosphere and the first results of such models.

We present high spatial and high temporal resolution observations of the quiet Sun in H-alpha obtained with the Swedish 1-m Solar Telescope on La Palma. We observe that many mottles, jet-like features in the quiet Sun, display clear up- and downward motions along their main axis. In addition, many mottles show vigorous transverse displacements. Unique identification of the mottles throughout their lifetime is much harder than for their active region counterpart, dynamic fibrils. This is because many seem to lack a sharply defined edge at their top, and significant fading often occurs throughout their lifetime. For those mottles that can be reliably tracked, we find that the mottle tops often undergo parabolic paths. We find a linear correlation between the deceleration these mottles undergo and the maximum velocity they reach, similar to what was found earlier for dynamic fibrils. Combined with an analysis of oscillatory properties, we conclude that at least part of the quiet Sun mottles are driven by magnetoacoustic shocks. In addition, the mixed polarity environment and vigorous dynamics suggest that reconnection may play a significant role in the formation of some quiet Sun jets.

With the recent advent of the Swedish 1-m Solar Telescope (SST), advanced image processing techniques, as well as numerical simulations that provide a more realistic view of the chromosphere, a comprehensive understanding of chromospheric jets such as spicules, mottles and fibrils is now within reach. In this paper, we briefly summarize results from a recent analysis of dynamic fibrils, short-lived jet-like features that dominate the chromosphere (as imaged in H-alpha) above and about active region plage. Using extremely high-resolution observations obtained at the SST, and advanced numerical 2D radiative MHD simulations, we show that fibrils are most likely formed by chromospheric shock waves that occur when convective flows and global oscillations leak into the chromosphere along the field lines of magnetic flux concentrations. In addition, we present some preliminary observations of quiet Sun jets or mottles. We find that the mechanism that produces fibrils in active regions is most likely also at work in quiet Sun regions, although it is modified by the weaker magnetic field and the presence of more mixed-polarity. A comparison with numerical simulations suggests that the weaker magnetic field in quiet Sun allows for significantly stronger (than in active regions) transverse motions that are superposed on the field-aligned, shock-driven motions. This leads to a more dynamic, and much more complex environment than in active region plage. In addition, our observations of the mixed polarity environment in quiet Sun regions suggest that other mechanisms, such as reconnection, may well play a significant role in the formation of some quiet Sun jets.

We present unprecedented high resolution H-alpha observations, obtained with the Swedish 1-m Solar Telescope, that, for the first time, spatially and temporally resolve dynamic fibrils in active regions on the Sun. These jet-like features are similar to mottles or spicules in quiet Sun. We find that most of these fibrils follow almost perfect parabolic paths in their ascent and descent. We measure the properties of the parabolic paths taken by 257 fibrils, and present an overview of the deceleration, maximum velocity, maximum length and duration, as well as their widths and the thickness of a bright ring that often occurs above dynamic fibrils. We find that the observed deceleration of the projected path is typically only a fraction of solar gravity, and incompatible with a ballistic path at solar gravity. We report on significant differences of fibril properties between those occurring above a dense plage region, and those above a less dense plage region where the magnetic field seems more inclined from the vertical. We compare these findings to advanced numerical 2D radiative MHD simulations, and find that fibrils are most likely formed by chromospheric shock waves that occur when convective flows and global oscillations leak into the chromosphere along the field lines of magnetic flux concentrations. Detailed comparison of observed and simulated fibril properties shows striking similarities of the values for deceleration, maximum velocity, maximum length and duration. We compare our results with observations of mottles and find that a similar mechanism is most likely at work in the quiet Sun.

We have obtained high resolution spectrograms of small scale magnetic structures with the Swedish 1-m Solar Telescope. We present Doppler measurements at $0\farcs{2}$ spatial resolution of bright points, ribbons and flowers and their immediate surroundings, in the C {\small{I}} 5380.3 {\AA} line (formed in the deep photosphere) and the two Fe {\small{I}} lines at 5379.6 {\AA} and 5386.3 {\AA}. The velocity inside the flowers and ribbons are measured to be almost zero, while we observe downflows at the edges. These downflows are increasing with decreasing height. We also analyze realistic magneto-convective simulations to obtain a better understanding of the interpretation of the observed signal. We calculate how the Doppler signal depends on the velocity field in various structures. Both the smearing effect of the non-negligible width of this velocity response function along the line of sight and of the smearing from the telescope and atmospheric point spread function are discussed. These studies lead us to the conclusion that the velocity inside the magnetic elements are really upflow of the order 1--2 km s${}^{-1}$ while the downflows at the edges really are much stronger than observed, of the order 1.5--3.3 km s${}^{-1}$.

We investigate the non-Local Thermodynamic Equilibrium (non-LTE) line formation of neutral carbon in late-type stars in order to remove some of the potential systematic errors in stellar abundance analyses employing C I features. The statistical equilibrium code MULTI was used on a grid of plane-parallel 1D MARCS atmospheric models. Within the parameter space explored, the high-excitation C I lines are stronger in non-LTE due to the combined effect of line source function drop and increased lower level population; the relative importance of the two effects depends on the particular combination of Teff, log g, [Fe/H] and [C/Fe] and on the analysed C I line. As a consequence, the non-LTE abundance corrections are negative and can be substantially so, for example ~-0.4 dex in halo turn-off stars at [Fe/H]~-3. The magnitude of the non-LTE corrections is rather insensitive to whether inelastic H collisions are included or not. Our results have implications on studies of nucleosynthetic processes and on Galactic chemical evolution models. When applying our calculated corrections to recent observational data, the upturn in [C/O] at low metallicity might still be present (thus apparently still necessitating contributions from massive Pop. III stars for the carbon production), but at a lower level and possibly with a rather shallow trend of ~-0.2 dex/dex below [O/H]~-1

The formation of jets such as dynamic fibrils, mottles, and spicules in the solar chromosphere is one of the most important, but also most poorly understood, phenomena of the Sun's magnetized outer atmosphere. We use extremely high-resolution observations from the Swedish 1-m Solar Telescope combined with advanced numerical modeling to show that in active regions these jets are a natural consequence of upwardly propagating slow mode magnetoacoustic shocks. These shocks form when waves generated by convective flows and global p-mode oscillations in the lower lying photosphere leak upward into the magnetized chromosphere. We find excellent agreement between observed and simulated jet velocities, decelerations, lifetimes and lengths. Our findings suggest that previous observations of quiet sun spicules and mottles may also be interpreted in light of a shock driven mechanism.

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.

The Ca II infrared triplet around 8540 A is a good candidate for observing chromospheric magnetism. We show results from combining a radiation hydrodynamic simulation with a Stokes synthesis code. The simulation shows interesting time-varying behavior of the Stokes V profiles as waves propagate through the lines' formation regions. Disappearing and reappearing lobes in the Stokes V profiles as well as profile asymmetries are closely related to the atmospheric velocity gradients.

We present non-Local Thermodynamic Equilibrium (non-LTE) calculations for neutral carbon spectral line formation, carried out for a grid of model atmospheres covering the range of late-type stars. The results of our detailed calculations suggest that the carbon non-LTE corrections in these stars are higher than usually adopted, remaining substantial even at low metallicity. For the most metal-poor stars in the sample of Akerman et al. (2004), the non-LTE abundance corrections are of the order of -0.35...-0.45 dex (when neglecting H collisions). Applying our results to those observations, the apparent [C/O] upturn seen in their LTE analysis is no longer present, thus revealing no need to invoke contributions from Pop. III stars to the carbon nucleosynthesis.

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

We present the first center-to-limb G-band images synthesized from high resolution simulations of solar magneto-convection. Towards the limb the simulations show "hilly" granulation with dark bands on the far side, bright granulation walls and striated faculae, similar to observations. At disk center G-band bright points are flanked by dark lanes. The increased brightness in magnetic elements is due to their lower density compared with the surrounding intergranular medium. One thus sees deeper layers where the temperature is higher. At a given geometric height, the magnetic elements are cooler than the surrounding medium. In the G-band, the contrast is further increased by the destruction of CH in the low density magnetic elements. The optical depth unity surface is very corrugated. Bright granules have their continuum optical depth unity 80 km above the mean surface, the magnetic elements 200-300 km below. The horizontal temperature gradient is especially large next to flux concentrations. When viewed at an angle, the deep magnetic elements optical surface is hidden by the granules and the bright points are no longer visible, except where the "magnetic valleys" are aligned with the line of sight. Towards the limb, the low density in the strong magnetic elements causes unit line-of-sight optical depth to occur deeper in the granule walls behind than for rays not going through magnetic elements and variations in the field strength produce a striated appearance in the bright granule walls.

Investigating the reliability of the assumption of instantaneous chemical equilibrium (ICE) for calculating the CO number density in the solar atmosphere is of crucial importance for the resolution of the long-standing controversy over the existence of `cool clouds' in the chromosphere, and for determining whether the cool gas owes its existence to CO radiative cooling or to a hydrodynamical process. Here we report the first results of such an investigation in which we have carried out time-dependent gas-phase chemistry calculations in radiation hydrodynamical simulations of solar chromospheric dynamics. We show that while the ICE approximation turns out to be suitable for modeling the observed infrared CO lines at the solar disk center, it may substantially overestimate the `heights of formation' of strong CO lines synthesized close to the edge of the solar disk, especially concerning vigorous dynamic cases resulting from relatively strong photospheric disturbances. This happens because during the cool phases of the hydrodynamical simulations the CO number density in the outer atmospheric regions is smaller than what is stipulated by the ICE approximation, resulting in decreased CO opacity in the solar chromosphere. As a result, the cool CO-bearing gas which produces the observed molecular lines must be located at atmospheric heights not greater than 700 km, approximately. We conclude that taking into account the non-equilibrium chemistry improves the agreement with the available on-disk and off-limb observations, but that the hydrodynamical simulation model has to be even cooler than anticipated by the ICE approximation, and this has to be the case at the `new' (i.e. deeper) formation regions of the rovibrational CO lines.

The lithium abundances in metal-poor halo stars are of importance for cosmology, galaxy evolution and stellar structure. In an attempt to study possible systematic errors in the derived Li abundances, the line formation of LiI lines has been investigated by means of realistic 3D hydrodynamical model atmospheres of halo stars and 3D non-LTE radiative transfer calculations. These are the first detailed 3D non-LTE computations reported employing a multi-level atomic model showing that such problems are now computationally tractable. The detailed computations reveal that the LiI population has a strong influence from the radiation field rather than the local gas temperature, indicating that the low derived Li abundances found by Asplund et al. (1999) are an artifact of their assumption of LTE. Relative to 3D LTE, the detailed calculations show pronounced over-ionization. In terms of abundances the 3D non-LTE values are within 0.05 dex of the 1D non-LTE results for the particular cases of HD140283 and HD84937, which is a consequence of the dominance of the radiation in determining the population density of LiI. Although 3D non-LTE can be expected to give results rather close (+/-0.1 dex) to 1D non-LTE for this reason, there may be systematic trends with metallicity and effective temperature.

We investigate the ionization of hydrogen in a dynamic Solar atmosphere. We show that the time scale for ionization/recombination can be estimated from the eigenvalues of a modified rate matrix where the optically thick Lyman transitions that are in detailed balance have been excluded. We find that the time scale for ionization/recombination is dominated by the slow collisional leakage from the ground state to the first excited state. Throughout the chromosphere the time scale is long ($10^3$-$10^5$ s), except in shocks where the increased temperature and density shorten the time scale for ionization/recombination, especially in the upper chromosphere. Because the relaxation time scale is much longer than dynamic time scales, hydrogen ionization does not have time to reach its equilibrium value and its fluctuations are much smaller than the variation of its statistical equilibrium value appropriate for the instantaneous conditions. The ionization state tends to represent the higher temperature of the shocks, and the mean electron density is up to a factor of six higher than the electron density calculated in statistical equilibrium from the mean atmosphere. The simulations show that a static picture and a dynamic picture of the chromosphere are fundamentally different and that time variations are crucial for our understanding of the chromosphere itself and the spectral features formed there.

On 12 March 1996 we obtained observations of the quiet Sun with the SUMER instrument. The observa- tions were sequences of 15-20 second exposures of ultraviolet emission line profiles and of the neighboring continua. These data contain signatures of the dynamics of the solar chromosphere that are uniquely useful because of wavelength coverage, moderate signal-to-noise ratios, and image stability. The dominant observed phenomenon is an oscillatory behavior that is analogous to the 3 minute oscillations seen in Ca II lines. The oscillations appear to be coherent over 3-8". At any time they occur over approx. 50 % of the area studied, and they appear as large perturbations in the intensities of lines and continua. The oscillations are most clearly seen in intensity variations in the UV (lambda > 912 A) continua, and they are also seen in the intensities and velocities of chromospheric lines of C I, N I and O I. Intensity brightenings are accompanied by blueshifts of typically 5 km s$^{-1}$. Phase differences between continuum and line intensities also indicate the presence of upward propagating waves. Three minute intensity oscillations are occasionally seen in second spectra (C II 1335), but never in third spectra (C III and Si III). Third spectra and He I 584 show oscillations in velocity that are not simply related to the 3 minute oscillations. The continuum intensity variations are consistent with recent simulations of chromospheric dynamics (Carlsson & Stein 1994) while the line observations indicate that important ingredients are missing at higher layers in the simulations. The data show that time variations are crucial for our understanding of the chromosphere itself and for the spectral features formed there.

We study the formation of B I lines in a grid of cool stellar model atmospheres without the assumption of local thermodynamic equilibrium (LTE). The non-LTE modelling includes the effect of other lines blending with the B I resonance lines. Except for the cases where the B I lines are very strong, the departures from LTE relevant for the resonance lines can be described as an overionisation effect and an optical-pumping effect. This causes the lines to be weaker than in LTE so that an abundance analysis assuming LTE will underestimate stellar boron abundances. We present non-LTE abundance corrections useful to improve on abundances derived from the B I 250 nm and 209 nm lines under the LTE assumption. Application of the results on literature data indicates that the B/Fe ratio in metal-poor stars is constant.

Enhanced chromospheric emission which corresponds to an outwardly increasing semiempirical temperature structure can be produced by wave motion without any increase in the mean gas temperture. Hence, the sun may not have a classical chromosphere in magnetic field free internetwork regions. Other significant differences between the properties of dynamic and static atmospheres should be considered when analyzing chromospheric observations.

We discuss the nature of non-LTE effects affecting abundance analysis of cool stars. The departures from LTE of importance for the B I lines in solar-type stars are described and some new results are presented. Boron abundances derived under the LTE assumption have significant systematic errors, especially for metal-poor stars. For beryllium, current results suggest that departures from LTE will not affect abundance analysis significantly.