We investigate the temperature enhancements and formation heights of impulsive heating phenomena in solar active regions such as Ellerman bombs (EBs), ultraviolet bursts (UVBs), and flaring active-region fibrils (FAFs) using interferometric observations in the millimeter (mm) continuum provided by the Atacama Large Millimeter/submillimeter Array (ALMA). We examined 3 mm signatures of heating events identified in Solar Dynamics Observatory (SDO) observations of an active region and compare the results with synthetic spectra from a 3D radiative magnetohydrodynamic simulation. We estimated the contribution from the corona to the mm brightness using differential emission measure analysis. We report the null detection of EBs in the 3 mm continuum at $\sim1.2$" spatial resolution, which is evidence that they are sub-canopy events that do not significantly contribute to heating the upper chromosphere. In contrast, we find the active region to be populated with multiple compact, bright, flickering mm bursts -- reminiscent of UVBs. The high brightness temperatures of up to $\sim14200$ K in some events may have a significant contribution (up to $\sim$7%) from the corona. We also detect FAF-like events in the 3 mm continuum that show rapid motions of $>10000$ K plasma launched with high plane-of-sky velocities ($37-340\rm\,km\,s^{-1}$) from nanoflare kernels. The mm FAFs are the brightest class of warm canopy fibrils that connect magnetic regions of opposite polarity. The simulation confirms that ALMA should be able to detect the mm counterparts of UVBs and small flares and thus provide a complementary diagnostic for impulsive heating in the solar chromosphere.

Broad-band images of the solar chromosphere in the Ca II H&K line cores around active regions are covered with fine bright elongated structures called bright fibrils. The mechanisms that form these structures and cause them to appear bright are still unknown. We aim to investigate the physical properties, such as temperature, line-of-sight velocity, and microturbulence, in the atmosphere that produces bright fibrils and to compare those to the properties of their surrounding atmosphere. We used simultaneous observations of a plage region in Fe I 6301-2 \r{A}, Ca II 8542 \r{A}, Ca II K, and H$\alpha$ acquired by the CRISP and CHROMIS instruments on the Swedish 1-m Solar Telescope. We manually selected a sample of 282 Ca II K bright fibrils. We compared the appearance of the fibrils in our sample to the Ca II 8542 \r{A} and H$\alpha$ data. We performed non-local thermodynamic equilibrium (non-LTE) inversions using the inversion code STiC on the Fe I 6301-2 \r{A}, Ca II 8542 \r{A}, Ca II K lines to infer the physical properties of the atmosphere. The line profiles in bright fibrils have a higher intensity in their K$_2$ peaks compared to profiles formed in the surrounding atmosphere. The inversion results show that the atmosphere in fibrils is on average $100-200$~K hotter at an optical depth log$(\tau) = -4.3$ compared to their surroundings. The line-of-sight velocity at chromospheric heights in the fibrils does not show any preference towards upflows or downflows. The microturbulence in the fibrils is on average 0.5 km s$^{-1}$ higher compared to their surroundings. Our results suggest that the fibrils have a limited extent in height, and they should be viewed as hot threads pervading the chromosphere.

Nearly all energy generated by fusion in the solar core is ultimately radiated away into space in the solar atmosphere, while the remaining energy is carried away in the form of neutrinos. The exchange of energy between the solar gas and the radiation field is thus an essential ingredient of atmospheric modeling. The equations describing these interactions are known, but their solution is so computationally expensive that they can only be solved in approximate form in multi-dimensional radiation-MHD modeling. In this review, I discuss the most commonly used approximations for energy exchange between gas and radiation in the photosphere, chromosphere, and corona.

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.

Numerical simulations of the solar chromosphere predict a diverse thermal structure with both hot and cool regions. Observations of plage regions in particular feature broader and brighter chromospheric lines, which suggest that they are formed in hotter and denser conditions than in the quiet-Sun, but also implies a non-thermal component whose source is unclear. We revisit the problem of the stratification of temperature and microturbulence in plage now adding millimeter continuum observations provided by ALMA to inversions of near-ultraviolet IRIS spectra as a powerful new diagnostic to disentangle the two parameters. We fit cool chromospheric holes and track the fast evolution of compact mm brightenings in the plage region. We use the STiC non-LTE inversion code to simultaneously fit real ultraviolet and millimeter spectra in order to infer the thermodynamic parameters of the plasma. We confirm the anticipated constraining potential of ALMA in non-LTE inversions of the solar chromosphere. We find significant differences between the inversion results of IRIS data alone compared to the results of a combination with the mm data: the IRIS+ALMA inversions have increased contrast and temperature range, and tend to favor lower values of microturbulence in the chromosphere of plage. The average brightness temperature of the plage region at 1.25 mm is 8500 K, but the ALMA maps also show much cooler ($\sim3000$ K) and hotter ($\sim11\,000$ K) evolving features partially seen in other diagnostics. To explain the former, the inversions require the existence of localized, low temperature regions in the chromosphere where molecules such as CO could form. The hot features could sustain such high temperatures due to non-equilibrium hydrogen ionization effects in a shocked chromosphere - a scenario that is supported by low-frequency shock wave patterns found in the MgII lines probed by IRIS.

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.

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.

Because of the complex physics that governs the formation of chromospheric lines, interpretation of solar chromospheric observations is difficult. The origin and characteristics of many chromospheric features are, because of this, unresolved. We focus here on studying two prominent features: long fibrils and flare ribbons. To model them, we use a 3D MHD simulation of an active region which self-consistently reproduces both of them. We model the H$\alpha$, Mg II k, Ca II K, and Ca II 8542 {\AA} lines using the 3D non-LTE radiative transfer code Multi3D. This simulation reproduces long fibrils that span between the opposite-polarity sunspots and go up to 4 Mm in height. They can be traced in all lines due to density corrugation. Opposite to previous studies, H$\alpha$, Mg II h&k, and Ca II H&K, are formed at similar height in this model. Magnetic field lines are aligned with the H$\alpha$ fibrils, but the latter holds to a lesser extent for the Ca II 8542 {\AA} line. The simulation shows structures in the H$\alpha$ line core that look like flare ribbons. The emission in the ribbons is caused by a dense chromosphere and a transition region at high column mass. The ribbons are visible in all chromospheric lines, but least prominent in Ca II 8542 {\AA} line. In some pixels, broad asymmetric profiles with a single emission peak are produced, similar to the profiles observed in flare ribbons. They are caused by a deep onset of the chromospheric temperature rise and large velocity gradients. The simulation produces long fibrils similar to what is seen in observations. It also produces structures similar to flare ribbons despite the lack of non-thermal electrons in the simulation. The latter suggests that thermal conduction might be a significant agent in transporting flare energy to the chromosphere in addition to non-thermal electrons.

The inference of the underlying state of the plasma in the solar chromosphere remains extremely challenging because of the nonlocal character of the observed radiation and plasma conditions in this layer. Inversion methods allow us to derive a model atmosphere that can reproduce the observed spectra by undertaking several physical assumptions. The most advanced approaches involve a depth-stratified model atmosphere described by temperature, line-of-sight velocity, turbulent velocity, the three components of the magnetic field vector, and gas and electron pressure. The parameters of the radiative transfer equation are computed from a solid ground of physical principles. To apply these techniques to spectral lines that sample the chromosphere, NLTE effects must be included in the calculations. We developed a new inversion code STiC to study spectral lines that sample the upper chromosphere. The code is based the RH synthetis code, which we modified to make the inversions faster and more stable. For the first time, STiC facilitates the processing of lines from multiple atoms in non-LTE, also including partial redistribution effects. Furthermore, we include a regularization strategy that allows for model atmospheres with a complex stratification, without introducing artifacts in the reconstructed physical parameters, which are usually manifested in the form of oscillatory behavior. This approach takes steps toward a node-less inversion, in which the value of the physical parameters at each grid point can be considered a free parameter. In this paper we discuss the implementation of the aforementioned techniques, the description of the model atmosphere, and the optimizations that we applied to the code. We carry out some numerical experiments to show the performance of the code and the regularization techniques that we implemented. We made STiC publicly available to the community.

Context. Magnetic reconnection during flares takes place in the corona but a substantial part of flare energy is deposited in the chromosphere. However, high-resolution spectro-polarimetric chromospheric observations of flares are very rare. The most used observables are Ca II 8542 {\AA} and He I 10830 {\AA}. Aims. We aim to study the chromosphere during a C3.6 class flare via spectro-polarimetric observations of the He I D3 line. Methods. We present the first SST/CRISP spectro-polarimetric observations of He I D3. We analyze the data using the inversion code HAZEL, and estimate the LOS velocity and the magnetic field vector. Results. Strong He I D3 emission at the flare footpoints, as well as strong He I D3 absorption profiles tracing the flaring loops are observed during the flare. The He I D3 traveling emission kernels at the flare footpoints exhibit strong chromospheric condensations of up to $\sim$ 60 km/s at their leading edge. Our observations suggest that such condensations result in shocking the deep chromosphere, causing broad and modestly blueshifted He I D3 profiles indicating subsequent upflows. A strong and rather vertical magnetic field of up to $\sim$ 2500 G is measured in the flare footpoints, confirming that the He I D3 line is likely formed in the deep chromosphere at those locations. We provide chromospheric line-of-sight velocity and magnetic field maps obtained via He I D3 inversions. We propose a fan-spine configuration as the flare magnetic field topology. Conclusions. The He I D3 line is an excellent diagnostic to study the chromosphere during flares. The impact of strong condensations on the deep chromosphere has been observed. Detailed maps of the flare dynamics and the magnetic field are obtained.

The European Solar Telescope (EST) is a project of a new-generation solar telescope. It has a large aperture of 4~m, which is necessary for achieving high spatial and temporal resolution. The high polarimetric sensitivity of the EST will allow to measure the magnetic field in the solar atmosphere with unprecedented precision. Here, we summarise the recent advancements in the realisation of the EST project regarding the hardware development and the refinement of the science requirements.

Unipolar magnetic regions are often associated with supergranular cells. The chromosphere above these regions is regulated by the magnetic field, but the field structure is poorly known. In unipolar regions, the fibrillar arrangement does not always coincide with magnetic field lines, and polarimetric observations are needed to establish the chromospheric magnetic topology. In an active region close to the limb, we observed a unipolar annular network of supergranular size. This supergranular structure harbours a radial distribution of the fibrils converging towards its centre. We observed this structure at different heights by taking data in the FeI 6301-6302 {\AA}, H-\alpha, CaII 8542 \AA\ and the CaII H&K spectral lines with the CRISP and CHROMIS instruments at the Swedish 1-m Solar Telescope. We performed Milne-Eddington inversions of the spectropolarimetric data of FeI and applied the weak field approximation to CaII 8542 \AA\ data to retrieve the magnetic field in the photosphere and chromosphere. We used magnetograms of CRISP, HINODE/SP and HMI to calculate the magnetic flux. We investigated the velocity using the line-of-sight velocities computed from the Milne-Eddington inversion and from Doppler shift of the K$_3$ feature in the CaII K spectral line. To describe the typical spectral profiles characterising the chromosphere above the supergranular structure, we performed a K-mean clustering of the spectra in CaIIK. The photospheric magnetic flux is not balanced. The magnetic field vector at chromospheric heights, retrieved by the weak field approximation, indicates that the field lines within the supegranular cell tend to point inwards, and might form a canopy above the unipolar region. In the centre of the supergranular cell hosting the unipolar region, we observe a persistent chromospheric brightening coinciding with a strong gradient in the line-of-sight velocity.

High-resolution observations of the solar chromosphere at millimeter wavelengths are now possible with the Atacama Large Millimeter Array (ALMA), promising to tackle many open problems in solar physics. Observations from other ground and space-based telescopes will greatly benefit from coordinated endeavors with ALMA, yet the diagnostic potential of combined optical, ultraviolet and mm observations has remained mostly unassessed. In this paper we investigate whether mm-wavelengths could aid current inversion schemes to retrieve a more accurate representation of the temperature structure of the solar atmosphere. We performed several non-LTE inversion experiments of the emergent spectra from a snapshot of 3D radiation-MHD simulation. We included common line diagnostics such as CaII H, K, 8542 \AA$\,$ and MgII h and k, taking into account partial frequency redistribution effects, along with the continuum around 1.2 mm and 3 mm. We found that including the mm-continuum in inversions allows a more accurate inference of temperature as function of optical depth. The addition of ALMA bands to other diagnostics should improve the accuracy of the inferred chromospheric temperatures between $\log \tau\sim[-6,-4.5]$ where the CaII and MgII lines are weakly coupled to the local conditions. However, we found that simultaneous multi-atom, non-LTE inversions of optical and UV lines present equally strong constraints in the lower chromosphere and thus are not greatly improved by the 1.2 mm band. Nonetheless, the 3 mm band is still needed to better constrain the mid-upper chromosphere.

Context: Radiation-MHD simulations have become a standard tool to investigate the physics of solar and stellar atmospheres. Aims: The aim of this paper to present a method that allows efficient and accurate analysis of flows in such simulations in the Lagrangian frame. Methods: This paper presents a method that allows the construction of pathlines given a seed point that can be chosen freely at any location and at any time during the simulation where the simulation state is stored. The method is based on passive tracer particles. Through injection of particles in expanding regions the occurrence of particle-free volumes is avoided, even in the case of strongly compressive flows. Results:The method was implemented in the solar and stellar atmosphere simulation code Bifrost. It is efficient and accurate. As examples I present an analysis of a gas parcel in the convection zone and a particle in the solar transition region.}

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.

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.

$\delta$-sunspots are known to be favourable locations for fast and energetic events like flares and CMEs. The photosphere of this type of sunspots has been thoroughly investigated in the past three decades. The atmospheric conditions in the chromosphere are not so well known, however. his study is focused on the chromosphere of a $\delta$-sunspot that harbours a series of fan-shaped jets in its penumbra . The aim of this study is to establish the magnetic field topology and the temperature distribution in the presence of jets in the photosphere and the chromosphere. We use data from the Swedish 1-m Solar Telescope (SST) and the Solar Dynamics Observatory. We invert the spectropolarimetric FeI 6302~\AA\ and CaII ~8542~\AA\ data from the SST using the the non-LTE inversion code NICOLE to estimate the magnetic field configuration, temperature and velocity structure in the chromosphere. A loop-like magnetic structure is observed to emerge in the penumbra of the sunspot. The jets are launched from the loop-like structure. Magnetic reconnection between this emerging field and the pre-existing vertical field is suggested by hot plasma patches on the interface between the two fields. The height at which the reconnection takes place is located between $\log \tau_{500} = -2$ and $\log \tau_{500} = -3$. The magnetic field vector and the atmospheric temperature maps show a stationary configuration during the whole observation.

The Chromospheric Lyman Alpha Spectropolarimeter (CLASP) observed the Sun in H I Lyman-{\alpha} during a suborbital rocket flight on September 3, 2015. The Interface Region Imaging Telescope (IRIS) coordinated with the CLASP observations and recorded nearly simultaneous and co-spatial observations in the Mg II h&k lines. The Mg II h and Ly-{\alpha} lines are important transitions, energetically and diagnostically, in the chromosphere. The canonical solar atmosphere model predicts that these lines form in close proximity to each other and so we expect that the line profiles will exhibit similar variability. In this analysis, we present these coordinated observations and discuss how the two profiles compare over a region of quiet sun at viewing angles that approach the limb. In addition to the observations, we synthesize both line profiles using a 3D radiation-MHD simulation. In the observations, we find that the peak width and the peak intensities are well correlated between the lines. For the simulation, we do not find the same relationship. We have attempted to mitigate the instrumental differences between IRIS and CLASP and to reproduce the instrumental factors in the synthetic profiles. The model indicates that formation heights of the lines differ in a somewhat regular fashion related to magnetic geometry. This variation explains to some degree the lack of correlation, observed and synthesized, between Mg II and Ly-{\alpha}. Our analysis will aid in the definition of future observatories that aim to link dynamics in the chromosphere and transition region.

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 present joint observations of the Sun by the Atacama Large Millimeter/submillimeter Array (ALMA) and the Interface Region Imaging Spectrograph (IRIS). The observations were made of a solar active region on 2015 December 18 as part of the ALMA science verification effort. A map of the Sun's continuum emission of size $2.4' \times 2.3'$ was obtained by ALMA at a wavelength of 1.25 mm (239 GHz) using mosaicing techniques. A contemporaneous map of size $1.9'\times 2.9'$ was obtained in the Mg II h doublet line at 2803.5\AA\ by IRIS. Both mm/submm$-\lambda$ continuum emission and ultraviolet (UV) line emission are believed to originate from the solar chromosphere and both have the potential to serve as powerful and complementary diagnostics of physical conditions in this poorly understood layer of the solar atmosphere. While a clear correlation between mm-$\lambda$ brightness temperature $T_B$ and the Mg II h line radiation temperature $T_{rad}$ is observed the slope is $<1$, perhaps as a result of the fact that these diagnostics are sensitive to different parts of the chromosphere and/or the Mg II h line source function includes a scattering component. There is a significant offset between the mean $T_B$(1.25 mm) and mean $T_{rad}$(Mg II), the former being $\approx 35\%$ greater than the latter. Partitioning the maps into "sunspot", "quiet regions", and "plage regions" we find that the slope of the scatter plots between the IRIS Mg II h line $T_{rad}$ and the ALMA brightness temperature $T_B$ is 0.4 (sunspot), 0.56 (quiet regions), and 0.66 (plage regions). We suggest that this change may be caused by the regional dependence of the formation heights of the IRIS and ALMA diagnostics, and/or the increased degree of coupling between the UV source function and the local gas temperature in the hotter, denser gas in plage regions.

The aim of this paper is to understand the formation of the OI line at 7772 \AA\ in the solar chromosphere. We used SST/CRISP observations to observe OI 7772 \AA\ in several places around the solar limb. We compared the observations with synthetic spectra calculated with the RH code in the one-dimension spherical geometry mode. New accurate hydrogen collisional rates were included for the RH calculations. The observations reveal a dark gap in the lower chromosphere, which is caused by variations in the line opacity as shown by our models. The lower level of the 7772 \AA\ transition is populated by a downward cascade from the continuum. We study the effect of Lyman-$\beta$ pumping and hydrogen collisions between the triplet and quintet system in OI. Both have a small but non-negligible influence on the line intensity.

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.

Our ability to model the shapes and strengths of iron lines in the solar spectrum is a critical test of the accuracy of the solar iron abundance, which sets the absolute zero-point of all stellar metallicities. We use an extensive 463-level Fe atom with new photoionisation cross-sections for FeI as well as quantum mechanical calculations of collisional excitation and charge transfer with neutral hydrogen; the latter effectively remove a free parameter that has hampered all previous line formation studies of Fe in non-local thermodynamic equilibrium (NLTE). For the first time, we use realistic 3D NLTE calculations of Fe for a quantitative comparison to solar observations. We confront our theoretical line profiles with observations taken at different viewing angles across the solar disk with the Swedish 1-m Solar Telescope. We find that 3D modelling well reproduces the observed centre-to-limb behaviour of spectral lines overall, but highlight aspects that may require further work, especially cross-sections for inelastic collisions with electrons. Our inferred solar iron abundance is log(eps(Fe))=7.48+-0.04.

3D non-LTE radiative transfer problems are computationally demanding, and this sets limits on the size of the problems that can be solved. So far Multilevel Accelerated Lambda Iteration (MALI) has been to the method of choice to perform high-resolution computations in multidimensional problems. The disadvantage of MALI is that its computing time scales as $\mathcal{O}(n^2)$, with $n$ the number of grid points. When the grid gets finer, the computational cost increases quadratically. We aim to develop a 3D non-LTE radiative transfer code that is more efficient than MALI. We implement a non-linear multigrid, fast approximation storage scheme, into the existing Multi3D radiative transfer code. We verify our multigrid implementation by comparing with MALI computations. We show that multigrid can be employed in realistic problems with snapshots from 3D radiative-MHD simulations as input atmospheres. With multigrid, we obtain a factor 3.3-4.5 speedup compared to MALI. With full-multigrid the speed-up increases to a factor 6. The speedup is expected to increase for input atmospheres with more grid points and finer grid spacing. Solving 3D non-LTE radiative transfer problems using non-linear multigrid methods can be applied to realistic atmospheres with a substantial speed-up.

Context. Ellerman bombs (EBs) are short-lived emission features, characterized by extended wing emission in hydrogen Balmer lines. Until now, no distinct signature of EBs has been found in the He I 10830 {\AA} line, and conclusive observations of EBs in He I D 3 have never been reported. Aims. We aim to study the signature of EBs in neutral helium triplet lines. Methods. The observations consist of 10 consecutive SST/TRIPPEL raster scans close to the limb, featuring the H$\beta$, He I D3 and He I 10830 {\AA} spectral regions. We also obtained raster scans with IRIS and make use of the SDO/AIA 1700 {\AA} channel. We use Hazel to invert the neutral helium triplet lines. Results. Three EBs in our data show distinct emission signatures in neutral helium triplet lines, most prominently visible in the He I D3 line. The helium lines have two components: a broad and blue-shifted emission component associated with the EB, and a narrower absorption component formed in the overlying chromosphere. One of the EBs in our data shows evidence of strong velocity gradients in its emission component. The emission component of the other two EBs could be fitted using a constant slab. Our analysis hints towards thermal Doppler motions having a large contribution to the broadening for helium and IRIS lines. We conclude that the EBs must have high temperatures to exhibit emission signals in neutral helium triplet lines. An order of magnitude estimate places our observed EBs in the range of $T\sim 2\cdot 10^4-10^5$ K.

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 Mg II h&k lines are powerful diagnostics for studying the solar chromosphere. They have become particularly popular with the launch of the IRIS satellite, and a number of studies that include these lines have lead to great progress in understanding chromospheric heating, in many cases thanks to the support from 3D MHD simulations. In this study we utilize another approach to analyze observations: non-LTE inversions of the Mg II h&k and UV triplet lines including the effects of partial redistribution. Our inversion code attempts to construct a model atmosphere that is compatible with the observed spectra. We have assessed the capabilities and limitations of the inversions using the FALC atmosphere and a snapshot from a 3D radiation-MHD simulation. We find that Mg II h&k allow reconstructing a model atmosphere from the middle photosphere to the transition region. We have also explored the capabilities of a multi-line/multi-atom setup, including the Mg II h&k, the Ca II 854.2 nm and the Fe I 630.25 lines to recover the full stratification of physical parameters, including the magnetic field vector, from the photosphere to the chromosphere. Finally, we present the first inversions of observed IRIS spectra from quiet-Sun, plage and sunspot, with very promising results.

Models of star formation in the early universe depend on the details of accretion, fragmentation and radiative feedback. Different simulations predict different initial mass functions of the first stars, ranging from predominantly low mass (0.1-10 Msol), to massive (10-100 Msol), or even supermassive (100-1000 Msol). The mass distribution of the first stars should lead to unique chemical imprints on the low-mass second and later generation metal-poor stars still in existence. The chemical composition of SMSS0313-6708, which has the lowest abundances of Ca and Fe of any star known, indicates it was enriched by a single massive supernova. However, even weak spectral lines may be affected by strong 3D and NLTE effects in metal-poor stars. If these effects are ignored or treated incorrectly, errors in the inferred abundances may significantly bias the inferred properties of the polluting supernovae. We redetermine the chemical composition of SMSS0313-6708 using 3D NLTE radiative transfer to obtain accurate abundances for Li, Na, Mg, Al, Ca and Fe. The model atoms employ realistic collisional rates, with no calibrated free parameters. We find significantly higher abundances in 3D NLTE than 1D LTE by 0.8 dex for Fe, and 0.5 dex for Mg, Al and Ca, while Li and Na are unaffected to within 0.03 dex. In particular, our upper limit for [Fe/H] is now a factor ten larger, at [Fe/H] < -6.53 (3 sigma), than previous estimates based on <3D> NLTE (i.e. using averaged 3D models). This higher estimate is due to a conservative upper limit estimation, updated NLTE data, and 3D-<3D> NLTE differences, all of which lead to a higher abundance determination. We find that the revised chemical composition of SMSS0313-6708 matches supernova yields for massive progenitors of 20-60 Msol exploding with low energies (1-2 x 10^51 erg), as well as progenitors of 10 Msol with very low explosion energies (< 10^51 erg).

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.

Resonance spectral lines such as H I Ly {\alpha}, Mg II h&k, and Ca II H&K that form in the solar chromosphere are influenced by the effects of 3D radiative transfer as well as partial redistribution (PRD). So far no one has modeled these lines including both effects simultaneously owing to the high computing demands of existing algorithms. Such modeling is however indispensable for accurate diagnostics of the chromosphere. We present a computationally tractable method to treat PRD scattering in 3D model atmospheres using a 3D non-LTE radiative transfer code. To make the method memory-friendly, we use the hybrid approximation of Leenaarts et al. (2012) for the redistribution integral. To make it fast, we use linear interpolation on equidistant frequency grids. We verify our algorithm against computations with the RH code and analyze it for stability, convergence, and usefulness of acceleration using model atoms of Mg II with the h&k lines and H I with the Ly {\alpha} line treated in PRD. A typical 3D PRD solution can be obtained in a model atmosphere with $252 \times 252 \times 496$ coordinate points in 50 000--200 000 CPU hours, which is a factor ten slower than computations assuming complete redistribution. We illustrate the importance of the joint action of PRD and 3D effects for the Mg II h&k lines for disk-center intensities as well as the center-to-limb variation. The proposed method allows simulating PRD lines in time series of radiation-MHD models in order to interpret observations of chromospheric lines at high spatial resolution.

We report on a fan-shaped set of high-speed jets above a strongly magnetized light bridge (LB) of a sunspot observed in the H$\alpha$ line. We study the origin, dynamics and thermal properties of the jets using high-resolution imaging spectroscopy in H$\alpha$ from the Swedish 1-m Solar Telescope and data from the Solar Dynamics Observatory and Hinode. The H$\alpha$ jets have lengths of 7-38 Mm, are impulsively accelerated to a speed of $\sim 100$ km/s close to photospheric footpoints in the LB, and exhibit a constant deceleration consistent with solar effective gravity. They are predominantly launched from one edge of the light bridge, and their footpoints appear bright in the H$\alpha$ wings. AIA data indicates elongated brightenings that are nearly co-spatial with the H$\alpha$ jets. We interpret them as jets of at least transition region temperatures. The magnetic field in the light bridge has a strength of $0.8-2$ kG and it is nearly horizontal. All jet properties are consistent with magnetic reconnection as the driver.

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 OI 777 nm lines are among the most commonly used diagnostics for the oxygen abundances in the atmospheres of FGK-type stars. However, they form in conditions that are far from local thermodynamic equilibrium (LTE). We explore the departures from LTE of atomic oxygen, and their impact on OI lines, across the Stagger-grid of three-dimensional hydrodynamic model atmospheres. For the OI 777 nm triplet we find significant departures from LTE. These departures are larger in stars with larger effective temperatures, smaller surface gravities, and larger oxygen abundances. We present grids of predicted 3D non-LTE based equivalent widths for the OI616nm, [OI] 630 nm, [OI] 636 nm, and OI 777 nm lines, as well as abundance corrections to 1D LTE based results.

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.

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.

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.

We revisit the Galactic chemical evolution of oxygen, addressing the systematic errors inherent in classical determinations of the oxygen abundance that arise from the use of one dimensional hydrostatic (1D) model atmospheres and from the assumption of local thermodynamic equilibrium (LTE). We perform detailed 3D non-LTE radiative transfer calculations for atomic oxygen lines across a grid of 3D hydrodynamic stag- ger model atmospheres for dwarfs and subgiants. We apply our grid of predicted line strengths of the [OI] 630 nm and OI 777 nm lines using accurate stellar parameters from the literature. We infer a steep decay in [O/Fe] for [Fe/H] $\gtrsim$ -1.0, a plateau [O/Fe] $\approx$ 0.5 down to [Fe/H] $\approx$ -2.5 and an increasing trend for [Fe/H] $\lesssim$ -2.5. Our 3D non-LTE calculations yield overall concordant results from the two oxygen abundance diagnostics.

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.

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.

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.

Recently, new solar model atmospheres have been developed to replace classical 1D LTE hydrostatic models and used to for example derive the solar chemical composition. We aim to test various models against key observational constraints. In particular, a 3D model used to derive the solar abundances, a 3D MHD model (with an imposed 10 mT vertical magnetic field), 1D models from the PHOENIX project, the 1D MARCS model, and the 1D semi-empirical model of Holweger & M\"uller. We confront the models with observational diagnostics of the temperature profile: continuum centre-to-limb variations (CLV), absolute continuum fluxes, and the wings of hydrogen lines. We also test the 3D models for the intensity distribution of the granulation and spectral line shapes. The predictions from the 3D model are in excellent agreement with the continuum CLV observations, performing even better than the Holweger & M\"uller model (constructed largely to fulfil such observations). The predictions of the 1D theoretical models are worse, given their steeper temperature gradients. For the continuum fluxes, predictions for most models agree well with the observations. No model fits all hydrogen lines perfectly, but again the 3D model comes ahead. The 3D model also reproduces the observed continuum intensity fluctuations and spectral line shapes very well. The excellent agreement of the 3D model with the observables reinforces the view that its temperature structure is realistic. It outperforms the MHD simulation in all diagnostics, implying that recent claims for revised abundances based on MHD modelling are premature. Several weaknesses in the 1D models are exposed. The differences between the PHOENIX LTE and NLTE models are small. We conclude that the 3D hydrodynamical model is superior to any of the tested 1D models, which gives further confidence in the solar abundance analyses based on it.

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.

Radiative transfer modeling of spectral lines including partial redistribution (PRD) effects requires the evaluation of the ratio of the emission to the absorption profile. This quantity requires a large amount of computational work if one employs the angle-dependent redistribution function, which prohibits its use in 3D radiative transfer computations with model atmospheres containing velocity fields. We aim to provide a method to compute the emission to absorption profile ratio that requires less computational work but retains the effect of angle-dependent scattering in the resulting line profiles. We present a method to compute the profile ratio that employs the angle-averaged redistribution function and wavelength transforms to and from the rest frame of the scattering particles. We compare the emergent line profiles of the \MgII\,k and \Lyalpha\ lines computed with angle-dependent PRD, angle-averaged PRD and our new method in two representative test atmospheres. The new method yields a good approximation of true angle-dependent profile ratio and the resulting emergent line profiles while keeping the computational speed and simplicity of angle-averaged PRD theory.

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.

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.

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.

Numerical simulations of the solar chromosphere have become increasingly realistic over the past 5 years. However, many observed chromospheric structures and behavior are not reproduced. Current models do not show fibrils in Ca II 8542, and neither reproduce the Ca II 8542 bisector. The emergent H-alpha line core intensity computed from the models show granulation instead of chromospheric shocks or fibrils. I discuss these deficiencies and speculate about what physics should be included to alleviate these shortcomings.

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.

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.

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.

The continuum intensity at wavelengths around 1 mm provides an excellent way to probe the solar chromosphere. Future high-resolution millimetre arrays, such as the Atacama Large Millimeter Array (ALMA), will thus produce valuable input for the ongoing controversy on the thermal structure and the dynamics of this layer. Synthetic brightness temperature maps are calculated on basis of three-dimensional radiation (magneto-)hydrodynamic (MHD) simulations. While the millimetre continuum at 0.3mm originates mainly from the upper photosphere, the longer wavelengths considered here map the low and middle chromosphere. The effective formation height increases generally with wavelength and also from disk-centre towards the solar limb. The average intensity contribution functions are usually rather broad and in some cases they are even double-peaked as there are contributions from hot shock waves and cool post-shock regions in the model chromosphere. Taking into account the deviations from ionisation equilibrium for hydrogen gives a less strong variation of the electron density and with it of the optical depth. The result is a narrower formation height range. The average brightness temperature increases with wavelength and towards the limb. The relative contrast depends on wavelength in the same way as the average intensity but decreases towards the limb. The dependence of the brightness temperature distribution on wavelength and disk-position can be explained with the differences in formation height and the variation of temperature fluctuations with height in the model atmospheres.

An approximate method for solving the rate equations for the hydrogen populations was extended and implemented in the three-dimensional radiation (magneto-)hydrodynamics code CO5BOLD. The method is based on a model atom with six energy levels and fixed radiative rates. It has been tested extensively in one-dimensional simulations. The extended method has been used to create a three-dimensional model that extends from the upper convection zone to the chromosphere. The ionisation degree of hydrogen in our time-dependent simulation is comparable to the corresponding equilibrium value up to 500 km above optical depth unity. Above this height, the non-equilibrium ionisation degree is fairly constant over time and space, and tends to be at a value set by hot propagating shock waves. The hydrogen level populations and electron density are much more constant than the corresponding values for statistical equilibrium, too. In contrast, the equilibrium ionisation degree varies by more than 20 orders of magnitude between hot, shocked regions and cool, non-shocked regions. The simulation shows for the first time in 3D that the chromospheric hydrogen ionisation degree and electron density cannot be calculated in equilibrium. Our simulation can provide realistic values of those quantities for detailed radiative transfer computations.

We investigate encounters between giant molecular clouds (GMCs) and star clusters. We propose a single expression for the energy gain of a cluster due to an encounter with a GMC, valid for all encounter distances and GMC properties. This relation is verified with N-body simulations of cluster-GMC encounters and excellent agreement is found. The fractional mass loss from the cluster is 0.25 times the fractional energy gain. This is because 75% of the injected energy goes to the velocities of escaping stars, that are higher than the escape velocity. We derive an expression for the cluster disruption time (t_dis) based on the mass loss from the simulations, taking into account the effect of gravitational focusing by the GMC. The disruption time depends on the cluster mass (M_c) and half-mass radius (r_h) as t_dis=2.0 S (M_c/10^4 M_sun)(3.75 pc/r_h)^3 Gyr, with S=1 for the solar neighbourhood and inversely proportional to the GMC density. The observed shallow relation between cluster radius and mass gives t_dis a power-law dependence on the mass with index 0.7, similar to that found from observations and from simulations of clusters dissolving in tidal fields (0.62). The constant of 2.0 Gyr is about a factor of 3.5 shorter than found from earlier simulations of clusters dissolving under the combined effect of galactic tidal field and stellar evolution. It is somewhat higher than the observationally determined value of 1.3 Gyr. It suggests, however, that the combined effect of tidal field and encounters with GMCs can explain the lack of old open clusters in the solar neighbourhood. GMC encounters can also explain the (very) short disruption time that was observed for star clusters in the central region of M51, since there rho_n is an order of magnitude higher than in the solar neighbourhood.