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.

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.

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.

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.