Bayesian Stokes inversion with normalizing flows

Díaz Baso, C. J.; Asensio Ramos, A.; de la Cruz Rodríguez, J.
Referencia bibliográfica

Astronomy and Astrophysics

Fecha de publicación:
3
2022
Número de autores
3
Número de autores del IAC
1
Número de citas
12
Número de citas referidas
9
Descripción
Stokes inversion techniques are very powerful methods for obtaining information on the thermodynamic and magnetic properties of solar and stellar atmospheres. In recent years, highly sophisticated inversion codes have been developed that are now routinely applied to spectro-polarimetric observations. Most of these inversion codes are designed to find an optimum solution to the nonlinear inverse problem. However, to obtain the location of potentially multimodal cases (ambiguities), the degeneracies and the uncertainties of each parameter inferred from the inversions algorithms - such as Markov chain Monte Carlo (MCMC) - require evaluation of the likelihood of the model thousand of times and are computationally costly. Variational methods are a quick alternative to Monte Carlo methods, and approximate the posterior distribution by a parametrized distribution. In this study, we introduce a highly flexible variational inference method to perform fast Bayesian inference, known as normalizing flows. Normalizing flows are a set of invertible, differentiable, and parametric transformations that convert a simple distribution into an approximation of any other complex distribution. If the transformations are conditioned on observations, the normalizing flows can be trained to return Bayesian posterior probability estimates for any observation. We illustrate the ability of the method using a simple Milne-Eddington model and a complex non-local thermodynamic equilibrium (NLTE) inversion. The method is extremely general and other more complex forward models can be applied. The training procedure need only be performed once for a given prior parameter space and the resulting network can then generate samples describing the posterior distribution several orders of magnitude faster than existing techniques.
Proyectos relacionados
Imagen del Proyecto
Magnestismo Solar y Estelar
Los campos magnéticos son uno de los ingredientes fundamentales en la formación de estrellas y su evolución. En el nacimiento de una estrella, los campos magnéticos llegan a frenar su rotación durante el colapso de la nube molecular, y en el fin de la vida de una estrella, el magnetismo puede ser clave en la forma en la que se pierden las capas
Tobías
Felipe García
Project Image
Magnetismo, Polarización y Transferencia Radiativa en Astrofísica
Los campos magnéticos están presentes en todos los plasmas astrofísicos y controlan la mayor parte de la variabilidad que se observa en el Universo a escalas temporales intermedias. Se encuentran en estrellas, a lo largo de todo el diagrama de Hertzsprung-Russell, en galaxias, e incluso quizás en el medio intergaláctico. La polarización de la luz
Tanausú del
Pino Alemán