Abstract
Coronal bright points (CBPs) are ubiquitous structures in the solar atmosphere composed of hot small-scale loops observed in extreme-ultraviolet (EUV) or X-rays in the quiet Sun and coronal ...holes. They are key elements to understanding the heating of the corona; nonetheless, basic questions regarding their heating mechanisms, the chromosphere underneath, or the effects of flux emergence in these structures remain open. We have used the Bifrost code to carry out a 2D experiment in which a coronal-hole magnetic null-point configuration evolves perturbed by realistic granulation. To compare with observations, synthetic SDO/AIA, Solar Orbiter EUI-HRI, and IRIS images have been computed. The experiment shows the self-consistent creation of a CBP through the action of stochastic granular motions alone, mediated by magnetic reconnection in the corona. The reconnection is intermittent and oscillatory, and it leads to coronal and transition-region temperature loops that are identifiable in our EUV/UV observables. During the CBP lifetime, convergence and cancellation at the surface of its underlying opposite polarities takes place. The chromosphere below the CBP shows a number of peculiar features concerning its density and the spicules in it. The final stage of the CBP is eruptive: Magnetic flux emergence at the granular scale disrupts the CBP topology, leading to different ejections, such as UV bursts, surges, and EUV coronal jets. Apart from explaining observed CBP features, our results pave the way for further studies combining simulations and coordinated observations in different atmospheric layers.
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 6563 and Ca ii 8542 . 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 surge and a Si iv burst that occurred on 2016 September 03 in active region AR 12585. To that end, we use coordinated observations from the Interface Region Imaging Spectrograph and the Swedish 1-m Solar Telescope. 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 nonequilibrium 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.
ABSTRACT Cool and dense ejections, typically H 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 accelerations close to . 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.
Abstract
The solar atmosphere is filled with clusters of hot small-scale loops commonly known as coronal bright points (CBPs). These ubiquitous structures stand out in the Sun by their strong X-ray ...and/or extreme-ultraviolet (EUV) emission for hours to days, which makes them a crucial piece when solving the solar coronal heating puzzle. In addition, they can be the source of coronal jets and small-scale filament eruptions. Here we present a novel 3D numerical model using the Bifrost code that explains the sustained CBP heating for several hours. We find that stochastic photospheric convective motions alone significantly stress the CBP magnetic field topology, leading to important Joule and viscous heating concentrated around the CBP’s inner spine at a few megameters above the solar surface. We also detect continuous upflows with faint EUV signals resembling observational dark coronal jets and small-scale eruptions when H
α
fibrils interact with the reconnection site. We validate our model by comparing simultaneous CBP observations from the Solar Dynamics Observatory (SDO) and the Swedish 1‐m Solar Telescope (SST) with observable diagnostics calculated from the numerical results for EUV wavelengths as well as for the H
α
line using the Multi3D synthesis code. Additionally, we provide synthetic observables to be compared with Hinode, Solar Orbiter, and the Interface Region Imaging Spectrograph (IRIS). Our results constitute a step forward in the understanding of the many different facets of the solar coronal heating problem.
Magnetic reconnection is thought to drive a wide variety of dynamic phenomena in the solar atmosphere. Yet, the detailed physical mechanisms driving reconnection are difficult to discern in the ...remote sensing observations that are used to study the solar atmosphere. In this Letter, we exploit the high-resolution instruments Interface Region Imaging Spectrograph and the new CHROMIS Fabry-Pérot instrument at the Swedish 1-m Solar Telescope (SST) to identify the intermittency of magnetic reconnection and its association with the formation of plasmoids in so-called UV bursts in the low solar atmosphere. The Si iv 1403 UV burst spectra from the transition region show evidence of highly broadened line profiles with often non-Gaussian and triangular shapes, in addition to signatures of bidirectional flows. Such profiles had previously been linked, in idealized numerical simulations, to magnetic reconnection driven by the plasmoid instability. Simultaneous CHROMIS images in the chromospheric Ca ii K 3934 line now provide compelling evidence for the presence of plasmoids by revealing highly dynamic and rapidly moving brightenings that are smaller than 0 2 and that evolve on timescales of the order of seconds. Our interpretation of the observations is supported by detailed comparisons with synthetic observables from advanced numerical simulations of magnetic reconnection and associated plasmoids in the chromosphere. Our results highlight how subarcsecond imaging spectroscopy sensitive to a wide range of temperatures combined with advanced numerical simulations that are realistic enough to compare with observations can directly reveal the small-scale physical processes that drive the wide range of phenomena in the solar atmosphere.
Ambipolar diffusion in the Bifrost code Nóbrega-Siverio, D.; Martínez-Sykora, J.; Moreno-Insertis, F. ...
Astronomy & astrophysics,
06/2020, Volume:
638
Journal Article
Peer reviewed
Open access
Context.
Ambipolar diffusion is a physical mechanism related to the drift between charged and neutral particles in a partially ionized plasma that is key to 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.
Aims.
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.
Methods.
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 the main feature of the module, we have implemented the super time stepping (STS) technique, which allows an important acceleration of the calculations. We have also introduced hyperdiffusion terms to guarantee the stability of the code.
Results.
We show that to have an accurate value for the ambipolar diffusion coefficient in the solar atmosphere it is necessary to include as atomic elements in the equation of state not only hydrogen and helium, but also the main electron donors like sodium, silicon, and potassium. In addition, we establish a range of criteria to set up an automatic selection of the free parameters of the STS method that guarantees the best performance, optimizing the stability and speed for the ambipolar diffusion calculations. We validate the STS implementation by comparison with a self-similar analytical solution.
Context.
Magnetic flux emergence from the solar interior has been shown to be a key mechanism for unleashing a wide variety of 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.
Aims.
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.
Methods.
Using the radiation-magnetohydrodynamic 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 H
2
molecules, and the ambipolar diffusion term of the generalized Ohm’s law.
Results.
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 H
2
molecules within the emerged region, thus mistakenly magnifying the exothermic contribution of the H
2
molecule formation to the thermal energy during the flux emergence process. We find that the ambipolar diffusion does not significantly affect the amount of total unsigned emerged magnetic flux, but it is important in the shocks that cross the emerged region, heating the plasma on characteristic times ranging from 0.1 to 100 s. We also briefly discuss the importance of including elements heavier than hydrogen in the equation of state so as not to overestimate the role of ambipolar diffusion in the atmosphere.
Context. Understanding the role of magnetic reconnection in the heating and dynamics of the solar atmosphere requires detailed observational data of any observable aspect of the reconnection process, ...including small-scale features such as plasmoids. Aims. Here, we examine the capability of active and upcoming instruments to detect plasmoids generated by reconnection in the corona including low-density regimes. Methods. We used the Bifrost code to perform simulations of plasmoid-mediated reconnection in the corona with a 2D idealized setup: a fan-spine topology with uniform density including thermal conduction. Through a forward-modeling of extreme-ultraviolet (EUV) observables, we checked whether our simulated plasmoids could be detected with the instruments of Solar Dynamics Observatory (SDO) and Solar Orbiter (SO), as well as the upcoming Multi-Slit Solar Explorer (MUSE) and Solar-C missions. Results. Short-lived (∼10 − 20 s) small-scale (∼0.2 − 0.5 Mm) coronal plasmoids are not resolvable with the Atmospheric Imaging Assembly (AIA) on board SDO. In contrast, they could be captured with the EUV High-Resolution Imager at the Extreme Ultraviolet Imager (EUI-HRI EUV ) of SO. The spatial and temporal high-resolution planned for the MUSE spectrograph (SG) is adequate to obtain full spectral information of these plasmoids. To achieve a sufficient signal-to-noise ratio (S/N) for ∼0.8 MK plasmoids in the MUSE/SG 171 Å channel, full-raster images are attainable for regions with electron densities above 10 9 cm −3 , while sit-and-stare observations are recommended for lower-density regions. The future Solar-C mission could also capture these coronal plasmoids using the EUV High-Throughput Spectroscopic Telescope (EUVST), considering rapid changes in Doppler shift and line widths in different EUV lines caused by plasmoid motions along the current sheet. Conclusions. With the combined spectra of MUSE/SG and Solar-C/EUVST in multiple emission lines, along with high-resolution images from SO/EUI-HRI EUV and MUSE/CI, it should be possible to gain new insights about plasmoid formation in the corona.
Advancements in instrumentation have revealed a multitude of small-scale extreme-ultraviolet (EUV) events in the solar atmosphere
and considerable effort is currently undergoing to unravel them. Our ...aim is to employ high-resolution and high-sensitivity magnetograms to gain a detailed understanding of the magnetic origin of such phenomena. We used coordinated observations from the Swedish 1-m Solar Telescope (SST), the Interface Region Imaging
Spectrograph ( and the Solar Dynamics Observatory ( to analyze an ephemeral magnetic flux emergence episode and the following chain of small-scale energetic events. These unique observations clearly link these phenomena together. The high-resolution (0 $) magnetograms obtained with SST/CRISP allowed us to reliably measure the
magnetic field at the photosphere and to detect the emerging bipole that caused the subsequent eruptive atmospheric events.
Notably, this small-scale emergence episode remains indiscernible in the lower
resolution SDO/HMI magnetograms (0 $).
We report the appearance of a dark bubble in related to the emerging bipole, a sign of the canonical expanding magnetic dome predicted in flux emergence simulations. Evidence of reconnection are also found, first through an Ellerman bomb and later by the launch of a surge next to a UV burst. The UV burst exhibits a weak EUV counterpart in the coronal SDO/AIA channels. By calculating the differential emission measure (DEM), its plasma is shown to reach a temperature beyond 1 MK and to have densities between the upper chromosphere and transition region. Our study showcases the importance of high-resolution magnetograms in revealing the mechanisms that trigger phenomena such as EBs, UV bursts, and surges.
This could hold implications for small-scale events similar to those recently reported in the EUV using Solar Orbiter.
The finding of temperatures beyond 1 MK in the UV burst plasma strongly suggests that we are examining analogous features.
Therefore, we recommend caution when drawing conclusions from full-disk magnetograms that lack the necessary resolution to reveal their true magnetic origin.
Context.
Magnetic reconnection is a fundamental mechanism in astrophysics. A common challenge in mimicking this process numerically in particular for the Sun is that the solar electrical resistivity ...is small compared to the diffusive effects caused by the discrete nature of codes.
Aims.
We aim to study different anomalous resistivity models and their respective effects on simulations related to magnetic reconnection in the Sun.
Methods.
We used the Bifrost code to perform a 2D numerical reconnection experiment in the corona that is driven by converging opposite polarities at the solar surface. This experiment was run with three different commonly used resistivity models: 1) the hyper-diffusion model originally implemented in Bifrost, 2) a resistivity proportional to the current density, and 3) a resistivity proportional to the square of the electron drift velocity. The study was complemented with a 1D experiment of a Harris current sheet with the same resistivity models.
Results.
The 2D experiment shows that the three resistivity models are capable of producing results in satisfactory agreement with each other in terms of the current sheet length, inflow velocity, and Poynting influx. Even though Petschek-like reconnection occurred with the current density-proportional resistivity while the other two cases mainly followed plasmoid-mediated reconnection, the large-scale evolution of thermodynamical quantities such as temperature and density are quite similar between the three cases. For the 1D experiment, some recalibration of the diffusion parameters is needed to obtain comparable results. Specifically the hyper-diffusion and the drift velocity-dependent resistivity model needed only minor adjustments, while the current density-proportional model needed a rescaling of several orders of magnitude.
Conclusions.
The Bifrost hyper-diffusion model is as suitable for simulations of magnetic reconnection as other common resistivity models and has the advantage of being applicable to any region in the solar atmosphere without the need for significant recalibration.