-
Hydrodynamic simulations of cool stellar atmospheres with \texttt{MANCHA}
Authors:
Andrea Perdomo García,
Nikola Vitas,
Elena Khomenko,
Manuel Collados
Abstract:
Aims. Our aim is to test how different binning strategies previously studied in one-dimensional models perform in three-dimensional radiative hydrodynamic simulations of stellar atmospheres. Methods. Realistic box-in-a-star simulations of the near-surface convection and photosphere of three spectral types (G2V, K0V, and M2V) were run with the MANCHA code with grey opacity. After reaching the stati…
▽ More
Aims. Our aim is to test how different binning strategies previously studied in one-dimensional models perform in three-dimensional radiative hydrodynamic simulations of stellar atmospheres. Methods. Realistic box-in-a-star simulations of the near-surface convection and photosphere of three spectral types (G2V, K0V, and M2V) were run with the MANCHA code with grey opacity. After reaching the stationary state, one snapshot of each of the three stellar simulations was used to compute the radiative energy exchange rate with grey opacity, opacity binned in four $τ$-bins, and opacity binned in 18 $\{ τ, λ\}$-bins. These rates were compared with the ones computed with opacity distribution functions. Then, stellar simulations were run with grey, four-bin, and 18-bin opacities to see the impact of the opacity setup on the mean stratification of the temperature and its gradient after time evolution. Results. The simulations of main sequence cool stars with the MANCHA code are consistent with those in the literature. For the three stars, the radiative energy exchange rates computed with 18 bins are remarkably close to the ones computed with the opacity distribution functions. The rates computed with four bins are similar to the rates computed with 18 bins, and present a significant improvement with respect to the rates computed with the Rosseland opacity, especially above the stellar surface. The Rosseland mean can reproduce the proper rates in sub-surface layers, but produces large errors for the atmospheric layers of the G2V and K0V stars. In the case of the M2V star, the Rosseland mean fails even in sub-surface layers, owing to the importance of the contribution from molecular lines in the opacity, underestimated by the harmonic mean. Similar conclusions are reached studying the mean stratification of the temperature and its gradient after time evolution.
△ Less
Submitted 4 June, 2024;
originally announced June 2024.
-
The influence of thermal pressure gradients and ionization (im)balance on the ambipolar diffusion and charge-neutral drifts
Authors:
M. M. Gómez-Míguez,
D. Martínez-Gómez,
E. Khomenko,
N. Vitas
Abstract:
Solar partially ionized plasma is frequently modeled using single-fluid (1F) or two-fluid (2F) approaches. In the 1F case, charge-neutral interactions are often described through ambipolar diffusion, while the 2F model fully considers charge-neutral drifts. Here, we expand the definition of the ambipolar diffusion coefficient to include inelastic collisions (ion/rec) in two cases: a VAL3C 1D model…
▽ More
Solar partially ionized plasma is frequently modeled using single-fluid (1F) or two-fluid (2F) approaches. In the 1F case, charge-neutral interactions are often described through ambipolar diffusion, while the 2F model fully considers charge-neutral drifts. Here, we expand the definition of the ambipolar diffusion coefficient to include inelastic collisions (ion/rec) in two cases: a VAL3C 1D model and a 2F simulations of the Rayleigh-Taylor instability (RTI) in a solar prominence thread based on \cite{PopLukKho2021aa, PopLukKho2021ab}. On one side, we evaluate the relative importance of the inelastic contribution, compared to elastic and charge-exchange collisions. On the other side, we compare the contributions of ion/rec, thermal pressure, viscosity, and magnetic forces to the charge-neutral drift velocity of the turbulent flow of the RTI. Our analysis reveals that the contribution of inelastic collisions to the ambipolar diffusion coefficient is negligible across the chromosphere, allowing the classical definition of this coefficient to be safely used in 1F modeling. However, in the transition region, the contribution of inelastic collisions can become as significant as that of elastic collisions. Furthermore, we ascertain that the thermal pressure force predominantly influences the charge-neutral drifts in the RTI model, surpassing the impact of the magnetic force.
△ Less
Submitted 21 February, 2024;
originally announced February 2024.
-
Mancha3D code: Multi-purpose Advanced Non-ideal MHD Code for High resolution simulations in Astrophysics
Authors:
M. Modestov,
E. Khomenko,
N. Vitas,
A. de Vicente,
A. Navarro,
P. A. Gonzalez-Morales,
M. Collados,
T. Felipe,
D. Martinez-Gomez,
P. Hunana,
M. Luna,
M. Koll Pistarini,
B. Popescu Braileanu,
A. Perdomo Garcia,
V. Liakh,
I. Santamaria,
M. M. Gomez Miguez
Abstract:
The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic processes in solar/stellar atmospheres. The code includes non-ideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high quality realistic simulations of magneto-convection, as well as idealized simulations of particular processes, suc…
▽ More
The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic processes in solar/stellar atmospheres. The code includes non-ideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high quality realistic simulations of magneto-convection, as well as idealized simulations of particular processes, such as wave propagation, instabilities or energetic events. The paper summarizes the equations and methods used in the Mancha3D code. It also describes its numerical stability and parallel performance and efficiency. The code is based on a finite difference discretization and memory-saving Runge-Kutta (RK) scheme. It handles non-ideal effects through super-time stepping and Hall diffusion schemes, and takes into account thermal conduction by solving an additional hyperbolic equation for the heat flux. The code is easily configurable to perform different kinds of simulations. Several examples of the code usage are given. It is demonstrated that splitting variables into equilibrium and perturbation parts is essential for simulations of wave propagation in a static background. A perfectly matched layer (PML) boundary condition built into the code greatly facilitates a non-reflective open boundary implementation. Spatial filtering is an important numerical remedy to eliminate grid-size perturbations enhancing the code stability. Parallel performance analysis reveals that the code is strongly memory bound, which is a natural consequence of the numerical techniques used, such as split variables and PML boundary conditions. Both strong and weak scalings show adequate performance up till several thousands of CPUs.
△ Less
Submitted 7 December, 2023;
originally announced December 2023.
-
Opacity for realistic 3D MHD simulations of cool stellar atmospheres
Authors:
A. Perdomo García,
N. Vitas,
E. Khomenko,
M. Collados,
C. Allende Prieto,
I. Hubeny,
Y. Osorio
Abstract:
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 o…
▽ More
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.
△ Less
Submitted 6 June, 2023;
originally announced June 2023.
-
A study of the capabilities for inferring atmospheric information from high-spatial-resolution simulations
Authors:
C. Quintero Noda,
E. Khomenko,
M. Collados,
B. Ruiz Cobo,
R. Gafeira,
N. Vitas,
M. Rempel,
R. J. Campbell,
A. Pastor Yabar,
H. Uitenbroek,
D. Orozco Suárez
Abstract:
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 obtai…
▽ More
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.
△ Less
Submitted 2 June, 2023;
originally announced June 2023.
-
Modeling the thermal conduction in the solar atmosphere with the code MANCHA3D
Authors:
Anamaría Navarro,
E. Khomenko,
M. Modestov,
N. Vitas
Abstract:
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 con…
▽ More
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.
△ Less
Submitted 20 May, 2022; v1 submitted 18 May, 2022;
originally announced May 2022.
-
Newly formed downflow lanes in exploding granules in the solar photosphere
Authors:
M. Ellwarth,
C. E. Fischer,
N. Vitas,
S. Schmiz,
W. Schmidt
Abstract:
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…
▽ More
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Å and the Fe I 5250.2Å 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.
△ Less
Submitted 1 July, 2021;
originally announced July 2021.
-
Influence of ambipolar and Hall effects on vorticity in 3D simulations of magneto-convection
Authors:
E. Khomenko,
M. Collados,
N. Vitas,
P. A. Gonzalez-Morales
Abstract:
This paper presents the results of the analysis of 3D simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation-runs are compared: one including both ambipolar diffusion and Hall effect; one including only ambipolar diffusion; and one without any of these two effects. The magnetic field is amplified from initial field to…
▽ More
This paper presents the results of the analysis of 3D simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation-runs are compared: one including both ambipolar diffusion and Hall effect; one including only ambipolar diffusion; and one without any of these two effects. The magnetic field is amplified from initial field to saturation level by the action of turbulent local dynamo. In each of these cases, we study 2 hours of simulated solar time after the local dynamo reaches the saturation regime. We analyze the power spectra of vorticity, of magnetic field fluctuations and of the different components of the magnetic Poynting flux responsible for the transport of vertical or horizontal perturbations. Our preliminary results show that the ambipolar diffusion produces a strong reduction of vorticity in the upper chromospheric layers and that it dissipates the vortical perturbations converting them into thermal energy. The Hall effect acts in the opposite way, strongly enhancing the vorticity. When the Hall effect is included, the magnetic field in the simulations becomes, on average, more vertical and long-lived flux tube-like structures are produced. We trace a single magnetic structure to study its evolution pattern and the magnetic field intensification, and their possible relation to the Hall effect.
△ Less
Submitted 21 September, 2020;
originally announced September 2020.
-
Joint action of Hall and ambipolar effects in 3D magneto-convection simulations of the quiet Sun. I. Dissipation and generation of waves
Authors:
P. A. Gonzalez-Morales,
E. Khomenko,
N. Vitas,
M. Collados
Abstract:
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 serie…
▽ More
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.
△ Less
Submitted 24 August, 2020;
originally announced August 2020.
-
Three-dimensional simulations of solar magneto-convection including effects of partial ionization
Authors:
E. Khomenko,
N. Vitas,
M. Collados,
A. de Vicente
Abstract:
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…
▽ More
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.
△ Less
Submitted 3 July, 2018;
originally announced July 2018.
-
Numerical simulations of quiet Sun magnetic fields seeded by Biermann battery
Authors:
E. Khomenko,
N. Vitas,
M. Collados,
A. de Vicente
Abstract:
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 sufficien…
▽ More
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.
△ Less
Submitted 19 June, 2017;
originally announced June 2017.
-
DeepVel: deep learning for the estimation of horizontal velocities at the solar surface
Authors:
A. Asensio Ramos,
I. S. Requerey,
N. Vitas
Abstract:
Many phenomena taking place in the solar photosphere are controlled by plasma motions. Although the line-of-sight component of the velocity can be estimated using the Doppler effect, we do not have direct spectroscopic access to the components that are perpendicular to the line-of-sight. These components are typically estimated using methods based on local correlation tracking. We have designed De…
▽ More
Many phenomena taking place in the solar photosphere are controlled by plasma motions. Although the line-of-sight component of the velocity can be estimated using the Doppler effect, we do not have direct spectroscopic access to the components that are perpendicular to the line-of-sight. These components are typically estimated using methods based on local correlation tracking. We have designed DeepVel, an end-to-end deep neural network that produces an estimation of the velocity at every single pixel and at every time step and at three different heights in the atmosphere from just two consecutive continuum images. We confront DeepVel with local correlation tracking, pointing out that they give very similar results in the time- and spatially-averaged cases. We use the network to study the evolution in height of the horizontal velocity field in fragmenting granules, supporting the buoyancy-braking mechanism for the formation of integranular lanes in these granules. We also show that DeepVel can capture very small vortices, so that we can potentially expand the scaling cascade of vortices to very small sizes and durations.
△ Less
Submitted 24 April, 2017; v1 submitted 15 March, 2017;
originally announced March 2017.
-
First detection of sign-reversed linear polarization from the forbidden [O I] 630.03 nm line
Authors:
A. G. de Wijn,
H. Socas-Navarro,
N. Vitas
Abstract:
We report on the detection of linear polarization of the forbidden [O i] 630.03 nm spectral line. The observations were carried out in the broader context of the determination of the solar oxygen abundance, an important problem in astrophysics that still remains unresolved. We obtained spectro-polarimetric data of the forbidden [O i] line at 630.03 nm as well as other neighboring permitted lines w…
▽ More
We report on the detection of linear polarization of the forbidden [O i] 630.03 nm spectral line. The observations were carried out in the broader context of the determination of the solar oxygen abundance, an important problem in astrophysics that still remains unresolved. We obtained spectro-polarimetric data of the forbidden [O i] line at 630.03 nm as well as other neighboring permitted lines with the Solar Optical Telescope of the Hinode satellite. A novel averaging technique was used, yielding very high signal-to-noise ratios in excess of $10^5$. We confirm that the linear polarization is sign-reversed compared to permitted lines as a result of the line being dominated by a magnetic dipole transition. Our observations open a new window for solar oxygen abundance studies, offering an alternative method to disentangle the Ni i blend from the [O i] line at 630.03 nm that has the advantage of simple LTE formation physics.
△ Less
Submitted 30 January, 2017;
originally announced January 2017.
-
How different are the Liège and Hamburg atlases of the solar spectrum?
Authors:
H. -P. Doerr,
N. Vitas,
D. Fabbian
Abstract:
Context: The high-fidelity solar spectral atlas prepared by Delbouille and co-workers (Liège atlas), and the one by Neckel and co-workers (Hamburg atlas), are widely recognised as the most important reference spectra of the Sun at disc-centre in the visible wavelength range. Both datasets serve as fundamental resources for many researchers, in particular for chemical abundance analysis. But despit…
▽ More
Context: The high-fidelity solar spectral atlas prepared by Delbouille and co-workers (Liège atlas), and the one by Neckel and co-workers (Hamburg atlas), are widely recognised as the most important reference spectra of the Sun at disc-centre in the visible wavelength range. Both datasets serve as fundamental resources for many researchers, in particular for chemical abundance analysis. But despite their similar published specifications (spectral resolution, noise level), the shapes of spectral lines in the two atlases differ significantly and systematically.
Aims: Knowledge of any instrumental degradations is imperative to fully exploit the information content of spectroscopic data. We seek to investigate the magnitude---and to explain the possible sources---of these differences. We provide the wavelength-dependent correction parameters that need to be taken into account when the spectra are to be compared with e.g. synthetic data.
Methods: A parametrically degraded version of the Hamburg spectrum is fitted to the Liège spectrum. The parameters of the model (wavelength shift, broadening, intensity scaling, intensity offset) represent the different characteristics of the respective instruments, observational strategies and data processing.
Results: The wavelength scales of the Liège and Hamburg atlases differ on average by 0.5 mÅ with a standard deviation of $\pm$2 mÅ, except for a peculiar region around 5500 Å. The continuum levels are offset by up to 18% below 5000 Å but stably at a 0.8% difference towards the red. We find no evidence for spectral straylight in the Liège spectrum. Its resolving power is almost independent of wavelength but limited to about 216000, between two to six times less than specified. When accounting for the degradations determined in this work, the spectra of the two atlases agree within a few parts in 10$^3$.
△ Less
Submitted 13 April, 2016;
originally announced April 2016.
-
Fluid description of multi-component solar partially ionized plasma
Authors:
Elena Khomenko,
Manuel Collados,
Antonio Diaz,
Nikola Vitas
Abstract:
We derive self-consistent formalism for the description of multi-component partially ionized solar plasma, by means of the coupled equations for the charged and neutral components for an arbitrary number of chemical species, and the radiation field. All approximations and assumptions are carefully considered. Generalized Ohm's law is derived for the single-fluid and two-fluid formalism. Our approa…
▽ More
We derive self-consistent formalism for the description of multi-component partially ionized solar plasma, by means of the coupled equations for the charged and neutral components for an arbitrary number of chemical species, and the radiation field. All approximations and assumptions are carefully considered. Generalized Ohm's law is derived for the single-fluid and two-fluid formalism. Our approach is analytical with some order-of-magnitude support calculations. After general equations are developed we particularize to some frequently considered cases as for the interaction of matter and radiation.
△ Less
Submitted 8 August, 2014;
originally announced August 2014.
-
Ellerman bombs: fallacies, fads, usage
Authors:
Robert J. Rutten,
Gregal J. M. Vissers,
Luc H. M. Rouppe van der Voort,
Peter Sütterlin,
Nikola Vitas
Abstract:
Ellerman bombs are short-lived brightenings of the outer wings of Halpha that occur in active regions with much flux emergence. We point out fads and fallacies in the extensive Ellerman bomb literature, discuss their appearance in various spectral diagnostics, and advocate their use as indicators of field reconfiguration in active-region topography using AIA 1700 A images.
Ellerman bombs are short-lived brightenings of the outer wings of Halpha that occur in active regions with much flux emergence. We point out fads and fallacies in the extensive Ellerman bomb literature, discuss their appearance in various spectral diagnostics, and advocate their use as indicators of field reconfiguration in active-region topography using AIA 1700 A images.
△ Less
Submitted 5 April, 2013; v1 submitted 4 April, 2013;
originally announced April 2013.
-
Line shape effects on intensity measurements of solar features: Brightness correction to SoHO MDI continuum images
Authors:
S. Criscuoli,
I. Ermolli,
D. Del Moro,
F. Giorgi,
A. Tritschler,
H. Uitenbroek,
N. Vitas
Abstract:
Continuum intensity observations obtained with the Michelson Doppler Imager (MDI) on-board the SoHO mission provide long time series of filtergrams that are ideal for studying the evolution of large-scale phenomena in the solar atmosphere and their dependence on solar activity. These filtergrams, however, are not taken in a pure continuum spectral band, but are constructed from a proxy, namely a c…
▽ More
Continuum intensity observations obtained with the Michelson Doppler Imager (MDI) on-board the SoHO mission provide long time series of filtergrams that are ideal for studying the evolution of large-scale phenomena in the solar atmosphere and their dependence on solar activity. These filtergrams, however, are not taken in a pure continuum spectral band, but are constructed from a proxy, namely a combination of filtergrams sampling the Ni I 676.8 nm line. We studied the sensitivity of this continuum proxy to the shape of the nickel line and to the degradation in the instrumental transmission profiles. We compared continuum intensity measurements in the nearby of nickel line with MDI proxy values in three sets of high resolution spectro-polarimetric data obtained with the Interferometric Bidimensional Spectrometer (IBIS), and in synthetic data, obtained from multi-dimensional simulations of magneto-convection and one-dimensional atmosphere models. We found that MDI continuum measurements require brightness corrections which depend on magnetic field strength, temperature and, to a smaller extent, plasma velocity. The correction ranges from 2% to 25% in sunspots, and is, on average, less than 2% for other features. The brightness correction also varies with position on the disk, with larger variations obtained for sunspots, and smaller variations obtained for quiet sun, faculae and micropores. Correction factors derived from observations agree with those deduced from the numerical simulations when observational effects are taken into account. Finally, we found that the investigated potential uncertainties in the transmission characteristics of MDI filters only slightly affect the brightness correction to proxy measurements.
△ Less
Submitted 12 December, 2010;
originally announced December 2010.
-
Explanation of the activity sensitivity of Mn I 5394.7 Å
Authors:
N. Vitas,
B. Viticchiè,
R. J. Rutten,
A. Vögler
Abstract:
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-di…
▽ More
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.
△ Less
Submitted 21 November, 2008;
originally announced November 2008.
-
On the solar abundance of indium
Authors:
N. Vitas,
I. Vince,
M. Lugaro,
O. Andriyenko,
M. Gosic,
R. J. Rutten
Abstract:
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 de…
▽ More
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.
△ Less
Submitted 14 November, 2007;
originally announced November 2007.