Three-dimensional simulations of solar magneto-convection including effects of partial ionization

Over the last decades, realistic 3D radiative-MHD simulations have become the dominant theoretical tool for understanding the complex interactions between the plasma and the magnetic field on the Sun. Most of such simulations are based on approximations of magnetohydrodynamics, without directly considering the consequences of the very low degree of ionization of the solar plasma in the photosphere and bottom chromosphere. The presence of large amount of neutrals leads to a partial decoupling of the plasma and the magnetic field. As a consequence of that, a series of non-ideal effects (ambipolar diffusion, Hall effect and battery effect) arises. The ambipolar effect is the dominant one in the solar chromosphere. Here we report on the first three-dimensional realistic simulations of magneto-convection including ambipolar diffusion and battery effects. The simulations are done using the newly developed Mancha3D code. Our results reveal that ambipolar diffusion causes measurable effects on the amplitudes of waves excited by convection in the simulations, on the absorption of Poynting flux and heating and on the formation of chromospheric structures. We provide a low limit on the chromospheric temperature increase due to the ambipolar effect using the simulations with battery-excited dynamo fields.Comment: To appear in Astronomy & Astrophysic

Numerical simulations of quiet Sun magnetic fields seeded by Biermann battery

The magnetic fields of the quiet Sun cover at any time more than 90\% of its surface and their magnetic energy budget is crucial to explain the thermal structure of the solar atmosphere. One of the possible origins of these fields is due to the action of local dynamo in the upper convection zone of the Sun. Existing simulations of the local solar dynamo require an initial seed field, and sufficiently high spatial resolution, in order to achieve the amplification of the seed field to the observed values in the quiet Sun. Here we report an alternative model of seeding based on the action of the Bierman battery effect. This effect generates a magnetic field due to the local imbalances in electron pressure in the partially ionized solar plasma. We show that the battery effect self-consistently creates from zero an initial seed field of a strength of the order of micro G, and together with dynamo amplification, allows the generation of quiet Sun magnetic fields of a similar strength to those from solar observations.Comment: To appear in Astronomy & Astrophysic

Modeling the thermal conduction in the solar atmosphere with the code MANCHA3D

Thermal conductivity is one of the important mechanisms of heat transfer in the solar corona. In the limit of strongly magnetized plasma, it is typically modeled by Spitzer's expression where the heat flux is aligned with the magnetic field. This paper describes the implementation of the heat conduction into the code MANCHA3D with an aim of extending single-fluid MHD simulations from the upper convection zone into the solar corona. Two different schemes to model heat conduction are implemented: (1) a standard scheme where a parabolic term is added to the energy equation, and (2) a scheme where the hyperbolic heat flux equation is solved. The first scheme limits the time step due to the explicit integration of a parabolic term, which makes the simulations computationally expensive. The second scheme solves the limitations on the time step by artificially limiting the heat conduction speed to computationally manageable values. The validation of both schemes is carried out with standard tests in one, two, and three spatial dimensions. Furthermore, we implement the model for heat flux derived by Braginskii (1965) in its most general form, when the expression for the heat flux depends on the ratio of the collisional to cyclotron frequencies of the plasma, and, therefore on the magnetic field strength. Additionally, our implementation takes into account the heat conduction in parallel, perpendicular, and transverse directions, and provides the contributions from ions and electrons separately. The model also transitions smoothly between field-aligned conductivity and isotropic conductivity for regions with a low or null magnetic field. Finally, we present a two-dimensional test for heat conduction using realistic values of the solar atmosphere where we prove the robustness of the two schemes implemented.Comment: 11 pages, 8 figure

Joint action of Hall and ambipolar effects in 3D magneto-convection simulations of the quiet Sun. I. Dissipation and generation of waves

The partial ionization of the solar plasma causes several nonideal effects such as the ambipolar diffusion, the Hall effect, and the Biermann battery effect. Here we report on the first three-dimensional realistic simulations of solar local dynamo where all three effects were taken into account. The simulations started with a snapshot of already saturated battery-seeded dynamo, where two new series were developed: one with solely ambipolar diffusion and another one also taking into account the Hall term in the generalized Ohm's law. The simulations were then run for about 4 hours of solar time to reach the stationary regime and improve the statistics. In parallel, a purely MHD dynamo simulation was also run for the same amount of time. The simulations are compared in a statistical way. The results show that, with the inclusion of the ambipolar diffusion, the amplitudes of the incompressible perturbations related to Alfven waves are reduced, and the Poynting flux is absorbed, with a frequency dependence. The Hall effect causes the opposite action: significant excess of incompressible perturbations is generated and an excess of the Poynting flux is observed in the chromospheric layers. The model with ambipolar diffusion shows, on average, sharper current sheets and slightly more abundant fast magneto-acoustic shocks in the chromosphere. The model with the Hall effect has higher temperatures at the lower chromosphere and stronger and more vertical magnetic field concentrations all over the chromosphere. The study of high-frequency waves reveals that significant power of incompressible perturbations is associated with areas with intense and more vertical magnetic fields and larger temperatures. We find a positive correlation between the magnitude of the ambipolar heating and the temperature increase at the same location after a characteristic time of 10^2 sec.Comment: To appear in Astronomy & Astrophysic

Newly formed downflow lanes in exploding granules in the solar photosphere

Exploding granules have drawn renewed interest because of their interaction with the magnetic field. Especially the newly forming downflow lanes developing in their centre seem to be eligible candidates for the intensification of magnetic fields. We analyse spectroscopic data from two different instruments in order to study the intricate velocity pattern within the newly forming downflow lanes in detail. We aim to examine general properties of a number of exploding granules. To gain a better understanding of the formation process of the developing intergranular lane in exploding granules, we study the temporal evolution and height dependence of the line-of-sight velocities at their formation location. Additionally, we search for evidence that exploding granules act as acoustic sources. We investigated the evolution of several exploding granules using data taken with the Interferometric Bidimensional Spectrometer and the Imaging Magnetograph eXperiment. Velocities for different heights of the solar atmosphere were determined by computing bisectors of the Fe I 6173.0{\AA} and the Fe I 5250.2{\AA} lines. We performed a wavelet analysis to study the intensity and velocity oscillations within and around exploding granules. We also compared our findings with predictions of numerical simulations. We found that exploding granules have significantly longer lifetimes than regular granules. Exploding granules larger than 3.8 arcsec form an independent intergranular lane during their decay phase, while smaller granules usually fade away or disappear into the intergranular area. For all exploding granules that form a new intergranular downflow lane, we find a temporal height-dependent shift with respect to the maximum of the downflow velocity. Our suggestion that this results from a complex atmospheric structure within the newly forming downflow lane is supported by the simulations.Comment: 13 pages; accepted for publication in A&

On the solar abundance of indium

The generally adopted value for the solar abundance of indium is over six times higher than the meteoritic value. We address this discrepancy through numerical synthesis of the 451.13 nm line on which all indium abundance studies are based, both for the quiet-sun and the sunspot umbra spectrum, employing standard atmosphere models and accounting for hyperfine structure and Zeeman splitting in detail. The results, as well as a re-appraisal of indium nucleosynthesis, suggest that the solar indium abundance is close to the meteoritic value, and that some unidentified ion line causes the 451.13 nm feature in the quiet-sun spectrum.Comment: 7 pages, 9 figures, accepted for publication in MNRAS Main Journa

Opacity for realistic 3D MHD simulations of cool stellar atmospheres

Context. Realistic 3D time-dependent simulations of stellar near-surface convection employ the opacity binning method for efficient and accurate computation of the radiative energy exchange. The method provides several orders of magnitude of speed-up, but its implementation includes a number of free parameters. Aims. Our aim is to evaluate the accuracy of the opacity binning method as a function of the choice of these free parameters. Methods. The monochromatic opacities computed with the SYNSPEC code are used to construct opacity distribution function (ODF) that is then verified through detailed comparison with the results of the ATLAS code. The opacity binning method is implemented with the SYNSPEC opacities for four representative cool main-sequence stellar spectral types (F3V, G2V, K0V, and M2V). Results. The ODFs from SYNSPEC and ATLAS show consistent results for the opacity and bolometric radiative energy exchange rate Q in case of the F, G, and K -- type stars. Significant differences, coming mainly from the molecular line lists, are found for the M -- type star. It is possible to optimise a small number of bins to reduce the deviation of the results coming from the opacity grouping with respect to the ODF for the F, G, and K -- type stars. In the case of the M -- type star, the inclusion of splitting in wavelength is needed in the grouping to get similar results, with a subsequent increase in computing time. In the limit of a large number of bins, the deviation for all the binning configurations tested saturates and the results do not converge to the ODF solution. Due to this saturation, the Q rate cannot be improved by increasing the number of bins to more than about 20 bins. The more effective strategy is to select the optimal location of fewer bins.Comment: 26 pages, 29 figures. Accepted for publication in section 15 (Numerical methods and codes) of Astronomy and Astrophysic

Explanation of the activity sensitivity of Mn I 5394.7 \AA

There is a long-standing controversy concerning the reason why the Mn I 5394.7 A line in the solar irradiance spectrum brightens more at larger activity than most other photospheric lines. The claim that this activity sensitivity is caused by spectral interlocking to chromospheric emission in Mg II h & k is disputed. Classical one-dimensional modeling is used for demonstration; modern three-dimensional MHD simulation for verification and analysis. The Mn I 5394.7 A line thanks its unusual sensitivity to solar activity to its hyperfine structure. This overrides the thermal and granular Doppler smearing through which the other, narrower, photospheric lines lose such sensitivity. We take the nearby Fe I 5395.2 A line as example of the latter and analyze the formation of both lines in detail to demonstrate and explain granular Doppler brightening. We show that this affects all narrow lines. Neither the chromosphere nor Mg II h & k play a role, nor is it correct to describe the activity sensitivity of Mn I 5394.7 A through plage models with outward increasing temperature contrast. The Mn I 5394.7 A line represents a proxy diagnostic of strong-field magnetic concentrations in the deep solar photosphere comparable to the G band and the blue wing of H-alpha, but not a better one than these. The Mn I lines are more promising as diagnostic of weak fields in high-resolution Stokes polarimetry.Comment: 12 pages, 8 figures, accepted by A&

A study of the capabilities for inferring atmospheric information from high-spatial-resolution simulations

In this work, we study the accuracy that can be achieved when inferring the atmospheric information from realistic numerical magneto-hydrodynamic simulations that reproduce the spatial resolution we will obtain with future observations made by the 4m class telescopes DKIST and EST. We first study multiple inversion configurations using the SIR code and the Fe I transitions at 630 nm until we obtain minor differences between the input and the inferred atmosphere in a wide range of heights. Also, we examine how the inversion accuracy depends on the noise level of the Stokes profiles. The results indicate that when the majority of the inverted pixels come from strongly magnetised areas, there are almost no restrictions in terms of the noise, obtaining good results for noise amplitudes up to 1$\times10^{-3}$ of $I_c$. At the same time, the situation is different for observations where the dominant magnetic structures are weak, and noise restraints are more demanding. Moreover, we find that the accuracy of the fits is almost the same as that obtained without noise when the noise levels are on the order of 1$\times10^{-4}$of $I_c$. We, therefore, advise aiming for noise values on the order of or lower than 5$\times10^{-4}$ of $I_c$ if observers seek reliable interpretations of the results for the magnetic field vector reliably. We expect those noise levels to be achievable by next-generation 4m class telescopes thanks to an optimised polarisation calibration and the large collecting area of the primary mirror.Comment: 14 pages, 13 figure