Context: Hot coronal jets are a basic observed feature of the solar atmosphere whose physical origin is still being actively debated. Aims: We study six recurrent jets occurring in the active region NOAA 12644 on April 04, 2017. They are observed in all the hot filters of AIA as well as cool surges in IRIS slit-jaw high spatial and temporal resolution images. Methods: The AIA filters allow us to study the temperature and the emission measure of the jets using the filter ratio method. We study the pre-jet phases by analyzing the intensity oscillations at the base of the jets with the wavelet technique. Results: A fine co-alignment of the AIA and IRIS data shows that the jets are initiated at the top of a canopy-like, double-chambered structure with cool emission on one side and hot emission in the other. The hot jets are collimated in the hot temperature filters, have high velocities (around 250 km/s) and accompanied by the cool surges and ejected kernels both moving at about 45 km/s. In the pre-phase of the jets, at their base we find quasi-periodic intensity oscillations in phase with small ejections; they have a period between 2 and 6 minutes and are reminiscent of acoustic or MHD waves. Conclusions: This series of jets and surges provides a good case-study to test the 2D and 3D magnetohydrodynamic (MHD) models that result from magnetic flux emergence. The double-chambered structure found in the observations corresponds to the cold and hot loop regions found in the models beneath the current sheet that contains the reconnection site. The cool surge with kernels is comparable with the cool ejection and plasmoids that naturally appear in the models.

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

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

The emergence of magnetic flux into the solar atmosphere produces, among other phenomena, Ellerman bombs (EBs), which are observed in H alpha and are due to magnetic reconnection in the photosphere below the chromospheric canopy. Signs of reconnection are also observed in other spectral lines, typical of the chromosphere or the transition region. An example are the UV bursts observed in the transition region lines of Si IV and the upper chromospheric lines of Mg II. In this work we analyze high cadence, high resolution coordinated observations between the Swedish 1-m Solar Telescope (SST) and the Interface Region Imaging Spectrograph (IRIS) spacecraft. H alpha images from the SST provide us with the positions, timings and trajectories of EBs in an emerging flux region. Simultaneous, co-aligned IRIS slit-jaw images at 133 (C II, transition region), 140 (Si IV, transition region) and 279.6 (Mg II k, core, upper chromosphere) nm, as well as spectroscopy in the far and near ultraviolet from the fast spectrograph raster, allow us to study the possible chromospheric/transition region counterparts of those EBs. Our main goal is to study the possible temporal and spatial relationship between several reconnection events at different layers in the atmosphere (namely EBs and UV bursts), the timing history between them, and the connection of these dynamical phenomena to the ejection of surges in the chromosphere. We also investigate the properties of an extended UV burst and their variations across the burst domain. Our results suggest a scenario where simultaneous and co-spatial EBs and UV bursts are part of the same reconnection system occurring sequentially along a vertical or nearly vertical current sheet. Heating and bidirectional jets trace the location where reconnection takes place. These results support and expand those obtained from recent numerical simulations of magnetic flux emergence.

Dynamical jets are generally found on Light bridges (LBs), which are key to studying sunspots decays. So far, their formation mechanism is not fully understood. In this paper, we used state-of-the-art observations from the Goode Solar Telescope, the Interface Region Imaging Spectrograph, the Spectro-Polarimeter on board Hinode and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory to analyze the fan shape jets on LBs in detail. Continuous upward motion of the jets in ascending phase is found from the H$\alpha$ velocity, which lasts for 12 minutes and is associated with the H$\alpha$ line wing enhancements. Two mini jets appear upon the bright fronts of the fan shape jets visible in the AIA 171 {\AA} and 193 {\AA} channels, with a time interval as short as 1 minute. Two kinds of small scale convective motions are identified in the photospheric images, along with the H$\alpha$ line wing enhancements. One seems to be associated with the formation of a new convection cell and the other manifests as the motion of a dark lane passing through the convection cell. The finding of three lobes Stokes V profiles and their inversion with NICOLE code indicates that there is magnetic field lines with opposite polarities in LBs. From the H$\alpha$ -0.8 {\AA} images, we found ribbon like brightenings propagating along the LBs, possibly indicating slipping reconnection. Our observation supports that the fan shape jets under study are caused by the magnetic reconnection and photospheric convective motions play an important role in triggering the magnetic reconnection.

Surges are ubiquitous cool ejections in the solar atmosphere that often appear associated with transient phenomena like UV bursts or coronal jets. Recent observations from the Interface Region Imaging Spectrograph (IRIS) show that surges, although traditionally related to chromospheric lines, can exhibit enhanced emission in Si IV with brighter spectral profiles than for the average transition region (TR). In this paper, we explain why surges are natural sites to show enhanced emissivity in TR lines. We performed 2.5D radiative-MHD numerical experiments using the Bifrost code including the nonequilibrium ionization of silicon and oxygen. A surge is obtained as a by-product of magnetic flux emergence; the TR enveloping the emerged domain is strongly affected by nonequilibrium effects: assuming statistical equilibrium would produce an absence of Si IV and O IV ions in most of the region. Studying the properties of the surge plasma emitting in the Si IV 1402.77 {\AA} and O IV 1401.16 {\AA} lines, we find that a) the timescales for the optically-thin losses and heat conduction are very short, leading to departures from statistical equilibrium, and b) the surge emits in Si IV more and has an emissivity ratio of Si IV to O IV larger than a standard TR. Using synthetic spectra, we conclude the importance of line-of-sight effects: given the involved geometry of the surge, the line-of-sight can cut the emitting layer at small angles and/or cross it multiple times, causing prominent, spatially intermittent brightenings both in Si IV and O IV.

Surges often appear as a result of the emergence of magnetized plasma from the solar interior. Traditionally, they are observed in chromospheric lines such as H$\alpha$ 6563 \AA and Ca II 8542 \AA. However, whether there is a response to the surge appearance and evolution in the Si IV lines or, in fact, in many other transition region lines has not been studied. In this paper we analyze a simultaneous episode of an H$\alpha$ surge and a Si IV burst that occurred on 2016 September 03 in active region AR12585. To that end, we use coordinated observations from the Interface Region Imaging Spectrograph (IRIS) and the Swedish 1-m Solar Telescope (SST). For the first time, we report emission of Si IV within the surge, finding profiles that are brighter and broader than the average. Furthermore, the brightest Si IV patches within the domain of the surge are located mainly near its footpoints. To understand the relation between the surges and the emission in transition region lines like Si IV, we have carried out 2.5D radiative MHD (RMHD) experiments of magnetic flux emergence episodes using the Bifrost code and including the non-equilibrium ionization of silicon. Through spectral synthesis we explain several features of the observations. We show that the presence of Si IV emission patches within the surge, their location near the surge footpoints and various observed spectral features are a natural consequence of the emergence of magnetized plasma from the interior to the atmosphere and the ensuing reconnection processes.

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

Cool and dense ejections, typically H$_{\alpha}$ surges, often appear alongside EUV or X-Ray coronal jets as a result of the emergence of magnetized plasma from the solar interior. Idealized numerical experiments explain those ejections as being indirectly associated with the magnetic reconnection taking place between the emerging and preexisting systems. However, those experiments miss basic elements that can importantly affect the surge phenomenon. In this paper we study the cool surges using a realistic treatment of the radiation transfer and material plasma properties. To that end, the Bifrost code is used, which has advanced modules for the equation of state of the plasma, photospheric and chromospheric radiation transfer, heat conduction and optically thin radiative cooling. We carry out a 2.5D experiment of the emergence of magnetized plasma through (meso)granular convection cells and the low atmosphere to the corona. Through detailed Lagrange tracing, we study the formation and evolution of the cool ejection and, in particular, the role of the entropy sources: this allows us to discern families of evolutionary patterns for the plasma elements. In the launch phase many elements suffer accelerations well in excess of gravity; when nearing the apex of their individual trajectories, instead, the plasma elements follow quasi-parabolic trajectories with acceleration close to g$_{\odot^{}}$. We show how the formation of the cool ejection is mediated by a wedge-like structure composed of two shocks, one of which leads to the detachment of the surge from the original emerged plasma dome.