Opacity for realistic 3D MHD simulations of cool stellar atmospheres

Perdomo García, A.; Vitas, N.; Khomenko, E.; Collados, M.; Allende Prieto, C.; Hubeny, I.; Osorio, Y.
Bibliographical reference

Astronomy and Astrophysics

Advertised on:
Number of authors
IAC number of authors
Refereed citations
Context. Realistic three-dimensional time-dependent simulations of stellar near-surface convection employ the opacity binning method for the efficient and accurate computation of the radiative energy exchange. The method provides several orders of magnitude of speedup, 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 were used to construct opacity distribution function (ODF) that was then verified through detailed comparison with the results of the ATLAS code. The opacity binning method was 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 the 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.
Related projects
Project Image
Magnetism, Polarization and Radiative Transfer in Astrophysics

Magnetic fields pervade all astrophysical plasmas and govern most of the variability in the Universe at intermediate time scales. They are present in stars across the whole Hertzsprung-Russell diagram, in galaxies, and even perhaps in the intergalactic medium. Polarized light provides the most reliable source of information at our disposal for the

Tanausú del
Pino Alemán
Project Image
Solar and Stellar Magnetism

Magnetic fields are at the base of star formation and stellar structure and evolution. When stars are born, magnetic fields brake the rotation during the collapse of the mollecular cloud. In the end of the life of a star, magnetic fields can play a key role in the form of the strong winds that lead to the last stages of stellar evolution. During

Felipe García
Solar Eruption
Numerical Simulation of Astrophysical Processes

Numerical simulation through complex computer codes has been a fundamental tool in physics and technology research for decades. The rapid growth of computing capabilities, coupled with significant advances in numerical mathematics, has made this branch of research accessible to medium-sized research centers, bridging the gap between theoretical and

Daniel Elías
Nóbrega Siverio
spectrum of mercury lamp
Chemical Abundances in Stars

Stellar spectroscopy allows us to determine the properties and chemical compositions of stars. From this information for stars of different ages in the Milky Way, it is possible to reconstruct the chemical evolution of the Galaxy, as well as the origin of the elements heavier than boron, created mainly in stellar interiors. It is also possible to

Allende Prieto